Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q ../SelR2-human.fa -t SelR2-NIVO01003597.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000022_1.0 # Protein # Methionine-R-sufoxide reductase 2 (SelR2) # Homo sapiens # Complete Target: NIVO01003597.1:subseq(1758717,113857) Ammotragus lervia isolate NAZ scaffold3385_len2400928_cov71, whole genome shotgun sequence Model: protein2genome:local Raw score: 738 Query range: 12 -> 198 Target range: 42432 -> 63869 13 : AlaAlaProGluArgAlaAlaArgLeuLeuTrpLeuLeuArgGlyLeuThrLeuGlyThr : 32 ||||||! !!!: ! !|||||||||||| !! !||||||||||||.!!||| !! ! AlaAlaLeuAspValMetAlaArgLeuLeuArgAlaLeuArgGlyLeuAlaLeuArgGlu 42433 : GCCGCGCTAGACGTCATGGCGCGCCTCCTGCGGGCTCTGCGGGGCCTGGCCCTCCGCGAG : 42490 33 : AlaProArgArgAlaValArgGlyGlnAlaGlyGlyGlyGlyProGlyThrGlyProGly : 52 |||||| !! !!.!!!.!||||||!:!|||! ! !!||||||! ! !!.!!||| !!||| AlaProGlyTrpThrAlaArgGlyArgAlaAspCysGlyGlyLeuArgAlaGlyAlaGly 42491 : GCGCCCGGGTGGACTGCGCGGGGCCGAGCGGACTGCGGCGGCCTCCGCGCAGGTGCTGGA : 42550 53 : LeuGlyGluAla{G} >>>> Target Intron 1 >>>> {ly}SerLeuAlaThr : 61 ! ! !:!!|||{|} 7432 bp {||}||||||.!!! ! SerProGlnAla{G}++ ++{ly}SerLeuThrLys 42551 : TCCCCGCAAGCA{G}gt.........................ag{GG}TCTCTTACAAAA : 50009 62 : CysGluLeuProLeuAlaLysSerGluTrpGlnLysLysLeuThrProGluGlnPheTyr : 81 !:!!! ! !!|||:!!|||||||||||||||||||||||||||||||||||||||:!! HisGlnProAlaLeuSerLysSerGluTrpGlnLysLysLeuThrProGluGlnPheHis 50010 : CATCAGCCGGCTCTCTCCAAGAGTGAGTGGCAGAAGAAGCTGACACCAGAGCAGTTCCAT : 50069 82 : ValThrArgGluLysGlyThrGluPro >>>> Target Intron 2 >>>> ProP : 92 ||||||||||||||||||||||||||| 5947 bp |||| ValThrArgGluLysGlyThrGluPro++ ++ProP 50070 : GTCACGAGAGAAAAAGGGACAGAGCCGgt.........................agCCTT : 56049 93 : heSerGlyIleTyrLeuAsnAsnLysGluAlaGlyMetTyrHisCysValCysCysAspS : 112 ||||||||||||||||||||||||||||| !!|||||||||||||||||||||||||||| heSerGlyIleTyrLeuAsnAsnLysGluProGlyMetTyrHisCysValCysCysAspS 56050 : TCAGTGGGATCTACCTGAATAACAAGGAACCAGGGATGTATCACTGTGTGTGTTGCGACA : 56109 113 : erProLeuPhe{Se} >>>> Target Intron 3 >>>> {r}SerGluLysLysT : 121 |||||||||||{||} 6504 bp {|}||||||||||||! erProLeuPhe{Se}++ ++{r}SerGluLysLysP 56110 : GCCCGCTCTTC{AG}gt.........................ag{C}TCTGAGAAGAAGT : 62640 122 : yrCysSerGlyThrGlyTrpProSerPheSerGluAlaHisGlyThrSerGlySerAspG : 141 :!|||||||||||||||||||||||||||||||||||||||||||||||||||||||||| heCysSerGlyThrGlyTrpProSerPheSerGluAlaHisGlyThrSerGlySerAspG 62641 : TCTGCTCTGGCACGGGATGGCCTTCGTTCTCTGAGGCTCACGGGACCTCTGGCTCTGATG : 62700 142 : luSerHisThrGlyIleLeuArgArgLeuAspThrSerLeuGlySerAlaArgThrGluV : 161 |||||:!!|||||||||||||||||| !||||||||||||||| !||||||||||||| luSerAsnThrGlyIleLeuArgArgAlaAspThrSerLeuGlyLeuAlaArgThrGluV 62701 : AGAGTAACACGGGGATCCTGAGACGTGCTGACACCTCGTTGGGACTAGCTCGCACAGAGG : 62760 162 : alValCysLysGln >>>> Target Intron 4 >>>> CysGluAlaHisLeuGl : 171 |||||||||||||| 996 bp ||||||||||||||||| alValCysLysGln++ ++CysGluAlaHisLeuGl 62761 : TCGTCTGTAAGCAGgt.........................agTGCGAAGCTCACCTTGG : 63786 172 : yHisValPheProAspGlyProGlyProAsnGlyGlnArgPheCysIleAsnSerValAl : 191 |||||||||||||||||||||||||||| !||||||||||||||||||||||||||||| yHisValPheProAspGlyProGlyProAlaGlyGlnArgPheCysIleAsnSerValAl 63787 : CCATGTGTTCCCCGATGGACCTGGGCCTGCTGGGCAGAGGTTCTGCATCAACAGTGTGGC : 63846 192 : aLeuLysPheLysProArgLys : 198 ||||!:!|||||||||!:!||| aLeuArgPheLysProLysLys 63847 : GCTCAGGTTCAAGCCAAAGAAA : 63869 vulgar: SPP00000022_1.0 12 198 . NIVO01003597.1:subseq(1758717,113857) 42432 63869 + 738 M 44 132 S 0 1 5 0 2 I 0 7428 3 0 2 S 1 2 M 33 99 5 0 2 I 0 5943 3 0 2 M 25 75 S 0 2 5 0 2 I 0 6500 3 0 2 S 1 1 M 49 147 5 0 2 I 0 992 3 0 2 M 33 99 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local gene 42433 63869 738 + . gene_id 1 ; sequence SPP00000022_1.0 ; gene_orientation + NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local cds 42433 42565 . + . NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local exon 42433 42565 . + . insertions 0 ; deletions 0 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice5 42566 42567 . + . intron_id 1 ; splice_site "GT" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local intron 42566 49997 . + . intron_id 1 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice3 49996 49997 . + . intron_id 0 ; splice_site "AG" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local cds 49998 50098 . + . NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local exon 49998 50098 . + . insertions 0 ; deletions 0 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice5 50099 50100 . + . intron_id 2 ; splice_site "GT" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local intron 50099 56045 . + . intron_id 2 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice3 56044 56045 . + . intron_id 1 ; splice_site "AG" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local cds 56046 56122 . + . NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local exon 56046 56122 . + . insertions 0 ; deletions 0 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice5 56123 56124 . + . intron_id 3 ; splice_site "GT" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local intron 56123 62626 . + . intron_id 3 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice3 62625 62626 . + . intron_id 2 ; splice_site "AG" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local cds 62627 62774 . + . NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local exon 62627 62774 . + . insertions 0 ; deletions 0 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice5 62775 62776 . + . intron_id 4 ; splice_site "GT" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local intron 62775 63770 . + . intron_id 4 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local splice3 63769 63770 . + . intron_id 3 ; splice_site "AG" NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local cds 63771 63869 . + . NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local exon 63771 63869 . + . insertions 0 ; deletions 0 NIVO01003597.1:subseq(1758717,113857) exonerate:protein2genome:local similarity 42433 63869 738 + . alignment_id 1 ; Query SPP00000022_1.0 ; Align 42433 13 132 ; Align 50000 58 99 ; Align 56046 91 75 ; Align 62628 117 147 ; Align 63771 166 99 # --- END OF GFF DUMP --- # -- completed exonerate analysis