Command line: [exonerate -m p2g --showtargetgff -q DIO1.x.fa -t DIO1.extraction.KV884771.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: DIO1 # Protein # Iodothyronine deiodinase 1 (DIO1) # Zebrafish Target: KV884771.1:subseq(62576,100000) Monopterus albus unplaced genomic scaffold scaffold83.1, whole genome shotgun sequence Model: protein2genome:local Raw score: 718 Query range: 0 -> 226 Target range: 49999 -> 51307 1 : LeuArgLysLeuPheValTyrIleSerAlaValLeuMetValCysAlaAlaIleLeuArg : 20 |||!:!||||||.!.:!!! !:!:|||.!!|||! !! |||! ! !!.!||| !! ! LeuGlnLysLeuMetLeuCysLeuSerThrValCysAsnValPheTyrValIleAlaLeu 50000 : TTGCAAAAACTTATGCTCTGTTTGTCGACAGTATGCAATGTCTTCTACGTTATAGCACTC : 50057 21 : MetSerMetLeuLysLeuLeuSerPheIleSerProGlyArgMetArgLysIleHisMet : 40 ! !! !:!!!:!:!!:!! !||||||||||||.!!! !! !:!||| ! !!:: AsnLeuThrValArgIleMetHisPheIleSerProSerLeuThrLysLysHisIleIle 50058 : AATCTCACGGTCAGGATTATGCATTTCATTTCACCAAGCCTCACCAAGAAGCACATCAtc : 50117 41 : LysMetGlyGluArgThrThrMetThrGlnAsnProLysPheArgTyrGluAspTrpGly : 60 :::|||...|||:::|||||||||||||||||||||||||||!:!||||||||||||||| ArgMetAsnGluLysThrThrMetThrGlnAsnProLysPheLysTyrGluAspTrpGly 50118 : agaatgaatgaaaagaCAACTATGACCCAGAATCCCAAATTTAAGTATGAAGATTGGGGT : 50177 61 : ProAlaPhePheSerLeuAlaPheIleLysThrLeuPhePheValAsnTrpCysSerLeu : 80 |||.!!||| !|||!!! !|||:!!|||.!! !! !!:: !! ||| |||||| ProThrPheAsnSerPheAsnPheValLysAlaAlaSerTrpHisMetTrpLeuSerLeu 50178 : CCGACGTTTAACTCCTTTAATTTCGTCAAAGCCGCCTCCTGGCACATGTGGCTGTCTCTC : 50237 81 : GlyLeuGluAlaPheGluGlyHisSerAlaProAspSerAlaLeuIleThrLeuAspArg : 100 |||! !||||||||| !! ! ! ..!|||||||||||| !!:!!:!!|||:!!||| !! GlyGlnGluAlaPheLeuGluLysGluAlaProAspSerProValValThrMetAspGly 50238 : GGACAAGAGGCTTTTCTGGAAAAAGAAGCTCCGGACTCACCTGTGGTCACCATGGACGGA : 50297 101 : GlnLysThrSerValHisArgPheLeuLys{G} >>>> Target Intron 1 >>>> : 111 :!!|||||||||:!! !!:!!:!|||!:!{!} 91 bp LysLysThrSerIlePheGlnTyrLeuArg{A}++ 50298 : AAGAAAACGAGTATCTTCCAGTACTTGAGA{G}gt......................... : 50333 112 : {ly}AsnArgProLeuValLeuSerPheGlySerCysThrUnkProProPheLeuTyr : 129 { !}|||||||||||||||||||||||||||||||||||| |||||||||:!!||| ++{sp}AsnArgProLeuValLeuSerPheGlySerCysThr***ProProPheMetTyr 50334 : ag{AC}AACAGGCCGCTGGTGTTGAGTTTTGGAAGCTGCACCTGACCCCCGTTTATGTAT : 50475 130 : LysLeuAspGluPheLysGlnLeuValLysAspPheSerAsnValAlaAspPheLeuIle : 149 !:!|||!!:||||||||||||||||||||||||||||||:!!|||||||||||||||:!: ArgLeuGluGluPheLysGlnLeuValLysAspPheSerAspValAlaAspPheLeuVal 50476 : AGACTTGAGGAGTTCAAACAACTCGTCAAGGACTTCAGCGATGTGGCGGACTTTCTTGTG : 50535 150 : ValTyrLeuAlaGluAlaHisAlaThr{A} >>>> Target Intron 2 >>>> { : 159 ||||||:!!||||||||||||:!!|||{|} 539 bp { ValTyrIleAlaGluAlaHisSerThr{A}++ ++{ 50536 : GTTTACATCGCCGAGGCACATTCAACA{G}gt.........................ag{ : 51105 160 : sp}AlaTrpAlaPheLysAsnAsnValAspIleSerValHisLysAsnLeuGluGluArg : 178 ||}!.!|||||||||!!.|||||| !||||||!.! !|||:::!:!||||||!!:||| sp}GlyTrpAlaPheAsnAsnAsnTyrAspIleLysHisHisArgSerLeuGluAspArg 51106 : AT}GGTTGGGCCTTTAACAACAACTATGACATCAAGCACCACCGTAGCCTGGAGGACAGG : 51161 179 : LeuAlaAlaAlaArgThrLeuLeuLysGluAspProProCysProValValValAspGlu : 198 |||:!!|||||||||! !|||:!!:!!:!!!!:|||! !||||||||||||||||||||| LeuSerAlaAlaArgIleLeuValGlnLysGluProLeuCysProValValValAspGlu 51162 : CTGTCCGCAGCACGGATCCTGGTTCAGAAGGAGCCCCTGTGTCCGGTGGTTGTGGATGAA : 51221 199 : MetAsnAsnIleThrAlaSerLysTyrGlyAlaLeuProGluArgLeuTyrValIleGln : 218 !!:|||:!!:!!||||||! !|||||||||||||||||||||||||||||||||:!:||| IleAsnAspValThrAlaIleLysTyrGlyAlaLeuProGluArgLeuTyrValLeuGln 51222 : ATAAATGATGTCACTGCCATAAAGTATGGTGCTCTACCTGAAAGGCTTTATGTGCTGCAG : 51281 219 : SerGlyLysValIleTyrGlnVal : 226 :!!|||||||||:!:|||:!!||| AlaGlyLysValValTyrLysVal 51282 : GCTGGGAAAGTCGTGTACAAGGTG : 51307 vulgar: DIO1 0 226 . KV884771.1:subseq(62576,100000) 49999 51307 + 718 M 110 330 S 0 1 5 0 2 I 0 87 3 0 2 S 1 2 M 47 141 S 0 1 5 0 2 I 0 535 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 2017-11-23 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884771.1:subseq(62576,100000) exonerate:protein2genome:local gene 50000 51307 718 + . gene_id 1 ; sequence DIO1 ; gene_orientation + KV884771.1:subseq(62576,100000) exonerate:protein2genome:local cds 50000 50330 . + . KV884771.1:subseq(62576,100000) exonerate:protein2genome:local exon 50000 50330 . + . insertions 0 ; deletions 0 KV884771.1:subseq(62576,100000) exonerate:protein2genome:local splice5 50331 50332 . + . intron_id 1 ; splice_site "GT" KV884771.1:subseq(62576,100000) exonerate:protein2genome:local intron 50331 50421 . + . intron_id 1 KV884771.1:subseq(62576,100000) exonerate:protein2genome:local splice3 50420 50421 . + . intron_id 0 ; splice_site "AG" KV884771.1:subseq(62576,100000) exonerate:protein2genome:local cds 50422 50565 . + . KV884771.1:subseq(62576,100000) exonerate:protein2genome:local exon 50422 50565 . + . insertions 0 ; deletions 0 KV884771.1:subseq(62576,100000) exonerate:protein2genome:local splice5 50566 50567 . + . intron_id 2 ; splice_site "GT" KV884771.1:subseq(62576,100000) exonerate:protein2genome:local intron 50566 51104 . + . intron_id 2 KV884771.1:subseq(62576,100000) exonerate:protein2genome:local splice3 51103 51104 . + . intron_id 1 ; splice_site "AG" KV884771.1:subseq(62576,100000) exonerate:protein2genome:local cds 51105 51307 . + . KV884771.1:subseq(62576,100000) exonerate:protein2genome:local exon 51105 51307 . + . insertions 0 ; deletions 0 KV884771.1:subseq(62576,100000) exonerate:protein2genome:local similarity 50000 51307 718 + . alignment_id 1 ; Query DIO1 ; Align 50000 1 330 ; Align 50424 112 141 ; Align 51107 160 201 # --- END OF GFF DUMP --- # -- completed exonerate analysis