Command line: [exonerate -m p2g --showtargetgff -q dio1_zf.aa.fa.canviat -t genomedio1_JXSJ01000014.1.fa] Hostname: [ws116069U] C4 Alignment: ------------ Query: SPP00000611_2.0 # Protein # Iodothyronine deiodinase 1 (DIO1) # Zebrafish Target: JXSJ01000014.1:subseq(2346940,101317) Miichthys miiuy scaffold14, whole genome shotgun sequence Model: protein2genome:local Raw score: 699 Query range: 0 -> 226 Target range: 49999 -> 51317 1 : LeuArgLysLeuPheValTyrIleSerAlaValLeuMetValCysAlaAlaIleLeuArg : 20 |||!:!!:!|||.!!:!!|||!!:|||.!!!.!! !! :!!||| ! !:!!! !! ! LeuGlnArgLeuLeuLeuTyrMetSerThrAlaTrpAsnLeuCysPheMetLeuSerLeu 50000 : TTGCAAAGGCTGCTTTTGTACATGTCGACAGCATGGAATCTTTGTTTTATGCTATCACTG : 50057 21 : MetSerMetLeuLysLeuLeuSerPheIleSerProGlyArgMetArgLysIleHisMet : 40 ! !!!:! !!:!:!!||| ||||||||||||...! ! !:!|||:!! !:!: AsnLeuIleCysArgIleLeuHisPheIleSerProSerLeuAlaLysLysValIleLeu 50058 : AATCTCATATGCAGGATCTTgcactttatttcaccaagtCTTGCCAAGAAGGTAATCCTC : 50117 41 : LysMetGlyGluArgThrThrMetThrGlnAsnProLysPheArgTyrGluAspTrpGly : 60 |||||||||||||||! !||||||||||||||||||...|||:::||||||||||||||| LysMetGlyGluArgIleThrMetThrGlnAsnProAsnPheLysTyrGluAspTrpGly 50118 : AAAATGGGTGAAAGGATCACCATGACCCAGAACCCCaattttaaatatgaagattGGGGT : 50177 61 : ProAlaPhePheSerLeuAlaPheIleLysThrLeuPhePheValAsnTrpCysSerLeu : 80 ! ||| !.!!||||||||||||! !!!. !! ||| |||||| ---LeuThrThrSerThrThrPheIleLysThrSerLeuLysHisMetTrpValSerLeu 50178 : ---CTGACAACGTCTACAACGTTCATCAAAACCTCGTTGAAGCACATGTGGGTGTCTCTT : 50234 81 : GlyLeuGluAlaPheGluGlyHisSerAlaProAspSerAlaLeuIleThrLeuAspArg : 100 |||! !|||||||||! !||| !..!|||||||||||| !!:!!:!!!:::::::: GlyGlnGluAlaPheValGlyGlyGluAlaProAspSerProValValSerMetGluGly 50235 : GGACAAGAGGCTTTTGTGGGAGGTGAAGCTCCAGACTCACCTGTGGTCAgcatggaggga : 50294 101 : GlnLysThrSerValHisArgPheLeuLys{G} >>>> Target Intron 1 >>>> : 111 :::|||:::|||:!! !!:!!:!||||||{|} 88 bp GluLysSerSerIleCysLysTyrLeuLys{G}++ 50295 : gagaaatcaAGTATCTGCAAGTACCTGAAA{G}gt......................... : 50330 112 : {ly}AsnArgProLeuValLeuSerPheGlySerCysThrUnkProProPheLeuTyr : 129 {||}|||||||||||||||||||||||||||||||||||| |||||||||:!!!:! ++{ly}AsnArgProLeuValLeuSerPheGlySerCysThr***ProProPheMetPhe 50331 : ag{gc}aaCAGACCGCTGGTGTTGAGTTTTGGGAGTTGCACCTGACCCCCGTTTATGTTC : 50469 130 : LysLeuAspGluPheLysGlnLeuValLysAspPheSerAsnValAlaAspPheLeuIle : 149 ||||||!!:||||||||||||||||||||||||||| !!:!!|||||||||||||||:!: LysLeuGluGluPheLysGlnLeuValLysAspPheArgAspValAlaAspPheLeuVal 50470 : AAACTTGAGGAGTTCAAGCAACTCGTCAAGGACTTCCGTGATGTAGCTGACTTTCTTGTG : 50529 150 : ValTyrLeuAlaGluAlaHisAlaThr{A} >>>> Target Intron 2 >>>> { : 159 :!!|||:!!||||||||||||:!!|||{|} 555 bp { IleTyrIleAlaGluAlaHisSerThr{A}++ ++{ 50530 : ATCTACATCGCTGAGGCACACTCAACA{G}gt.........................ag{ : 51115 160 : sp}AlaTrpAlaPheLysAsnAsnValAspIleSerValHisLysAsnLeuGluGluArg : 178 ||}!.!|||||||||! |||||| !!||||||!:! !|||::!||||||||||||||| sp}GlyTrpAlaPheThrAsnAsnPheAspIleAsnGlnHisArgAsnLeuGluGluArg 51116 : AT}GGTTGGGCCTTTACCAACAACTTTGACATCAATCAGCACCGGAACCTGGAGGAAAGG : 51171 179 : LeuAlaAlaAlaArgThrLeuLeuLysGluAspProProCysProValValValAspGlu : 198 |||:!!||||||!:!! !|||:!!:!!:!!!!:|||! !||||||||||||||||||||| LeuSerAlaAlaGlnIleLeuValGlnLysGluProLeuCysProValValValAspGlu 51172 : TTGTCCGCTGCACAGATCCTGGTACAGAAAGAGCCTCTGTGTCCAGTGGTTgtggatgaa : 51231 199 : MetAsnAsnIleThrAlaSerLysTyrGlyAlaLeuProGluArgLeuTyrValIleGln : 218 |||...||| ..!|||! !|||||||||||||||||||||||||||||||||:!:||| MetLysAsnAlaAlaAlaIleLysTyrGlyAlaLeuProGluArgLeuTyrValLeuGln 51232 : atgaaaaatgcagcTGCCATAAAGTATGGTGCTCTGCCAGAGAGGCTTTATGTGCTGCAG : 51291 219 : SerGlyLysValIleTyrGlnVal : 226 :!!|||||||||:!:|||:!!||| AlaGlyLysValValTyrLysVal 51292 : GCTGGGAAAGTTGTGTACAAGGTG : 51317 vulgar: SPP00000611_2.0 0 226 . JXSJ01000014.1:subseq(2346940,101317) 49999 51317 + 699 M 60 180 G 1 0 M 49 147 S 0 1 5 0 2 I 0 84 3 0 2 S 1 2 M 47 141 S 0 1 5 0 2 I 0 551 3 0 2 S 1 2 M 67 201 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-18 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local gene 50000 51317 699 + . gene_id 1 ; sequence SPP00000611_2.0 ; gene_orientation + JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local cds 50000 50327 . + . JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local exon 50000 50327 . + . insertions 0 ; deletions 1 JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local splice5 50328 50329 . + . intron_id 1 ; splice_site "GT" JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local intron 50328 50415 . + . intron_id 1 JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local splice3 50414 50415 . + . intron_id 0 ; splice_site "ag" JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local cds 50416 50559 . + . JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local exon 50416 50559 . + . insertions 0 ; deletions 0 JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local splice5 50560 50561 . + . intron_id 2 ; splice_site "GT" JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local intron 50560 51114 . + . intron_id 2 JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local splice3 51113 51114 . + . intron_id 1 ; splice_site "AG" JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local cds 51115 51317 . + . JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local exon 51115 51317 . + . insertions 0 ; deletions 0 JXSJ01000014.1:subseq(2346940,101317) exonerate:protein2genome:local similarity 50000 51317 699 + . alignment_id 1 ; Query SPP00000611_2.0 ; Align 50000 1 180 ; Align 50180 62 147 ; Align 50418 112 141 ; Align 51117 160 201 # --- END OF GFF DUMP --- # -- completed exonerate analysis