Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q ../DI2-human.fa -t DI2-NIVO01066248.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000002_1.0 # Protein # Iodothyronine deiodinase 2 (DI2) # Homo sapiens # Complete Target: NIVO01066248.1:subseq(6673153,109011) Ammotragus lervia isolate NAZ scaffold348_len6808047_cov72, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 1271 Query range: 0 -> 265 Target range: 59011 -> 49999 1 : MetGlyIleLeuSerValAspLeuLeuIleThrLeuGlnIleLeuProValPhePheSer : 20 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| MetGlyIleLeuSerValAspLeuLeuIleThrLeuGlnIleLeuProValPhePheSer 59011 : ATGGGCATCCTCAGCGTAGACTTGCTGATCACACTGCAGATTCTGCCAGTTTTTTTCTCC : 58954 21 : AsnCysLeuPheLeuAlaLeuTyrAspSerValIleLeuLeuLysHisValValLeuLeu : 40 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AsnCysLeuPheLeuAlaLeuTyrAspSerValIleLeuLeuLysHisValValLeuLeu 58953 : AACTGCCTCTTCCTGGCGCTCTATGACTCGGTCATTCTCCTCAAGCACGTGGTGCTGCTG : 58894 41 : LeuSerArgSerLysSerThrArgGlyGluTrpArgArgMetLeuThrSerGluGlyLeu : 60 |||||||||||||||||||||||||||:!!|||||||||||||||||||||||||||:!! LeuSerArgSerLysSerThrArgGlyGlnTrpArgArgMetLeuThrSerGluGlyMet 58893 : CTGAGCCGCTCCAAGTCCACTCGCGGGCAGTGGAGGCGCATGCTGACCTCAGAAGGAATG : 58834 61 : ArgCysValTrpLysSerPheLeuLeuAspAlaTyrLysGln >>>> Target Intr : 75 ||||||:!!||||||||||||||||||||||||||||||||| 8217 bp ArgCysIleTrpLysSerPheLeuLeuAspAlaTyrLysGln++ 58833 : CGCTGCATCTGGAAAAGCTTCCTCCTGGACGCCTACAAACAGgt................ : 58787 76 : on 1 >>>> ValLysLeuGlyGluAspAlaProAsnSerSerValValHisValSerS : 91 |||||||||||||||||||||||||||||| !!||||||! !||||||| ++ValLysLeuGlyGluAspAlaProAsnSerArgValValProValSerS 58786 : .........agGTGAAACTGGGTGAAGATGCCCCCAATTCCCGCGTGGTGCCTGTCTCCA : 50524 92 : erThrGluGlyGlyAspAsnSerGlyAsnGlyThrGlnGluLysIleAlaGluGlyAlaT : 111 || !!||||||||||||!.!||||||||||||.!!|||||||||! !|||!!:||| !! erProGluGlyGlyAspThrSerGlyAsnGlyAlaGlnGluLysThrAlaAspGlyProG 50523 : GTCCGGAAGGAGGTGACACCAGTGGAAATGGTGCCCAGGAGAAAACAGCGGATGGACCTG : 50464 112 : hrCysHisLeuLeuAspPheAlaSerProGluArgProLeuValValAsnPheGlySerA : 131 !|||||||||||||||||||||||||||||||||||||||||||||||||||||||||| luCysHisLeuLeuAspPheAlaSerProGluArgProLeuValValAsnPheGlySerA 50463 : AATGCCACCTTCTGGACTTTGCCAGCCCTGAGCGCCCGCTGGTGGTCAACTTCGGCTCGG : 50404 132 : laThrUnkProProPheThrSerGlnLeuProAlaPheArgLysLeuValGluGluPheS : 151 ||||| ||||||||||||!:!|||||||||||||||!!!||||||||||||||||||| laThr***ProProPheThrAsnGlnLeuProAlaPheSerLysLeuValGluGluPheS 50403 : CCACTTGACCTCCTTTCACTAACCAGCTGCCAGCCTTCAGCAAACTGGTGGAAGAGTTCT : 50344 152 : erSerValAlaAspPheLeuLeuValTyrIleAspGluAlaHisProSerAspGlyTrpA : 171 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| erSerValAlaAspPheLeuLeuValTyrIleAspGluAlaHisProSerAspGlyTrpA 50343 : CATCCGTGGCTGACTTCCTGTTGGTCTACATTGATGAGGCTCATCCTTCAGATGGTTGGG : 50284 172 : laIleProGlyAspSerSerLeuSerPheGluValLysLysHisGlnAsnGlnGluAspA : 191 ||:!:|||||||||||||||||| ||||||||||||||||||!:!||||||||||||| laValProGlyAspSerSerLeuPhePheGluValLysLysHisArgAsnGlnGluAspA 50283 : CGGTGCCTGGggattcttctttgtttttcgaGGTGAAGAAGCACCGGAACCAGGAAGACC : 50224 192 : rgCysAlaAlaAlaGlnGlnLeuLeuGluArgPheSerLeuProProGlnCysArgValV : 211 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| rgCysAlaAlaAlaGlnGlnLeuLeuGluArgPheSerLeuProProGlnCysArgValV 50223 : GATGCGCAGCAGCCCAGCAGCTCCTGGAGCGTTTCTCCTTGCCGCCCCAGTGCCGAGTTG : 50164 212 : alAlaAspArgMetAspAsnAsnAlaAsnIleAlaTyrGlyValAlaPheGluArgValC : 231 |||||||||||||||||||||||||||||:!!|||||||||||||||||||||||||||| alAlaAspArgMetAspAsnAsnAlaAsnValAlaTyrGlyValAlaPheGluArgValC 50163 : TGGCTGACCGCATGGACAATAATGCCAACGTCGCTTATGGGGTAGCCTTTGAACGCGTGT : 50104 232 : ysIleValGlnArgGlnLysIleAlaTyrLeuGlyGlyLysGlyProPheSerTyrAsnL : 251 ||||||||||||||||||||||||||||||||||||||||||||||||||! !||||||| ysIleValGlnArgGlnLysIleAlaTyrLeuGlyGlyLysGlyProPhePheTyrAsnL 50103 : GCATTGTGCAGAGACAGAAAATTGCTTATCTGGGAGGAAAGGGTCCCTTCTTTTACAACC : 50044 252 : euGlnGluValArgHisTrpLeuGluLysAsnPheSerLysArg : 265 ||||||||||||||!.!||||||||||||||||||||||||||| euGlnGluValArgArgTrpLeuGluLysAsnPheSerLysArg 50043 : TGCAAGAAGTCCGGCGTTGGCTGGAGAAGAATTTCAGCAAGAGA : 50000 vulgar: SPP00000002_1.0 0 265 . NIVO01066248.1:subseq(6673153,109011) 59011 49999 - 1271 M 74 222 5 0 2 I 0 8213 3 0 2 M 191 573 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-16 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local gene 50000 59011 1271 - . gene_id 1 ; sequence SPP00000002_1.0 ; gene_orientation + NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local cds 58790 59011 . - . NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local exon 58790 59011 . - . insertions 0 ; deletions 0 NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local splice5 58788 58789 . - . intron_id 1 ; splice_site "GT" NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local intron 50573 58789 . - . intron_id 1 NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local splice3 50573 50574 . - . intron_id 0 ; splice_site "AG" NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local cds 50000 50572 . - . NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local exon 50000 50572 . - . insertions 0 ; deletions 0 NIVO01066248.1:subseq(6673153,109011) exonerate:protein2genome:local similarity 50000 59011 1271 - . alignment_id 1 ; Query SPP00000002_1.0 ; Align 59012 1 222 ; Align 50573 75 573 # --- END OF GFF DUMP --- # -- completed exonerate analysis