Command line: [exonerate -m p2g --showtargetgff -q DIO1.fa -t genomic.fa] Hostname: [ws101087U] C4 Alignment: ------------ Query: tr|D3Y4D3|D3Y4D3_DANRE Iodothyronine deiodinase OS=Danio rerio PE=2 SV=1 Target: CM004438.1:subseq(17080000,150000) Ictalurus punctatus breed USDA103 chromosome 25, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 933 Query range: 11 -> 264 Target range: 4805 -> 4052 12 : SerAlaLeuLysAsnAlaAlaValCysValLeuLeuLeuProArgPheLeuLeuAlaAlaLe : 32 !.!||||||||||||||| !|||! !:!!:!!|||||||||||||||||||||.!!|||:! LysAlaLeuLysAsnAlaLeuValSerLeuMetLeuLeuProArgPheLeuLeuThrAlaVa 4805 : AAAGCGCTCAAAAACGCGCTGGTGTCCCTGATGCTGCTGCCGCGCTTCCTGCTGACCGCCGT : 4745 33 : uMetLeuCysLeuLeuAspPheLeuCysIleArgArgLysLeuLeuLeuLysMetGlnGluG : 53 !:!!|||!! !!!|||||||||:!!|||||||||||||||:!!|||:!!||||||!::||| lLeuLeuTrpPheLeuAspPheValCysIleArgArgLysValLeuValLysMetArgGluA 4744 : GCTGCTGTGGTTTCTCGACTTCGTCTGCATCCGACGGAAGGTGCTGGTGAAGATGCGCGAAC : 4682 54 : lyAlaPheSerSerProAspAspProProLeuArgValSerAspSerAsnLysMetPheThr : 73 !!! ! !|||||||||||||||||| !!:!!! !|||||||||||||||::!||||||||| rgGluGlySerSerProAspAspProSerValThrValSerAspSerAsnArgMetPheThr 4681 : GCGAGGGCTCGAGTCCCGACGACCCGTCGGTGACGGTGTCAGACTCCAACCGGATGTTCACG : 4622 74 : LeuGluSerLeuArgAlaValTrpTyrGlyGlnLysLeuAspPhePheLysSerAlaArgLe : 94 :!!||||||||||||||||||||||||.!!|||||||||||||||! !|||:!!|||!.!|| ValGluSerLeuArgAlaValTrpTyrSerGlnLysLeuAspPheCysLysAlaAlaHisLe 4621 : GTGGAGTCACTGCGCGCCGTGTGGTACAGCCAGAAGTTGGACTTCTGCAAGGCGGCGCACCT : 4559 95 : uGlyGlyAlaAlaProAsnThrGluValPheProLeuAspGlyAspAlaArgAlaAlaGluA : 115 |||| !.!!|||||||||!:!|||||| ! |||||| |||||| !!:! !||| | uGlyLeuThrAlaProAsnSerGluValValProLeu---GlyAspArgLysArgAla---A 4558 : GGGGCTGACGGCGCCCAACAGCGAGGTGGTACCGCTC---GGCGATCGGAAACGCGCG---C : 4502 116 : rgIleLeuAspTyrAlaArgGlyArgArgProLeuIleLeuAsnPheGlySerCysSerUnk : 135 ||||||||||||||||||||||| !||||||||||||||||||||||||||||||||| rgIleLeuAspTyrAlaArgGlySerArgProLeuIleLeuAsnPheGlySerCysSer*** 4501 : GCATACTGGACTACGCCCGCGGCTCGAGACCGCTCATCCTTAACTTCGGCAGCTGCTCATGA : 4442 136 : ProProPheMetThrArgLeuSerAlaPheGlnArgValAlaArgGlnTyrAlaAspIleAl : 156 |||||||||||||||||||||:!!||||||!:!||||||||| !||||||||||||||||| ProProPheMetThrArgLeuAlaAlaPheArgArgValAlaAspGlnTyrAlaAspIleAl 4441 : CCACCGTTCATGACGCGCCTGGCCGCGTTCCGGCGCGTTGCTGACCAGTACGCCGACATCGC : 4379 157 : aAspSerLeuLeuValTyrIleGluGluAlaHisProProAspGlyTrpValSerSerAspA : 177 |||||||||||||:!!||||||||||||||||||||| !!|||||||||||||||!!!|||| aAspSerLeuLeuIleTyrIleGluGluAlaHisProSerAspGlyTrpValSerThrAspA 4378 : CGACTCACTGCTTATATACATCGAGGAGGCGCATCCGTCCGACGGCTGGGTGAGCACGGACG : 4316 178 : laProValGlnIleProArgHisArgCysLeuGluAspArgLeuArgAlaAlaGlnMetLeu : 197 ||||| !|||||||||||||||||||||:!!:!!||||||||||||||||||!:::!::!! laProTyrGlnIleProArgHisArgCysIleLysAspArgLeuArgAlaAlaArgLeuMet 4315 : CCCCGTACCAGATCCCACGCCACCGCTGTATAAAGGACCGTCTCCGAGCCGCGCGCCTCATG : 4256 198 : HisArgGluAlaAlaGlySerAlaValValValAspSerMetGlnAsnSerCysAsnAlaAl : 218 ! ! !!.! !!|||||| !||||||||||||!!!|||.!.||||||! !|||||||| ThrAlaThrValProGlySerIleValValValAspThrMetAspAsnSerSerAsnAlaAl 4255 : ACCGCCACGGTACCAGGCAGCATCGTGGTCGTGGACACCATGGACAACTCGTCCAACGCGGC : 4193 219 : aTyrGlyAlaTyrPheGluArgLeuTyrIleValLysAspAlaAlaValValTyrGlnGlyG : 239 ||||||||||||||||||||||||||||:!:|||||||||! ! !|||||||||||||||| aTyrGlyAlaTyrPheGluArgLeuTyrValValLysAspGluLysValValTyrGlnGlyG 4192 : GTACGGAGCCTACTTCGAGCGGCTCTACGTGGTCAAGGACGAGAAGGTGGTGTATCAGGGAG : 4130 240 : lyArgGlyProGluGlyTyrArgIleAlaGluLeuArgAspTrpLeuGluArgTyrArgSer : 259 ||||||||||||||||||||||||||:!!|||||||||:!!|||||||||!:!||||||..! lyArgGlyProGluGlyTyrArgIleSerGluLeuArgAsnTrpLeuGluGlnTyrArgGln 4129 : GCCGGGGACCCGAGGGCTACCGCATTTCCGAACTGAGGAACTGGCTCGAGCAATATCGCCAA : 4070 260 : GlyLeuGlnGluSer : 264 !|||:!!! ||| LysLeuGluGlySer 4069 : AAGCTCGAAGGCTCC : 4053 vulgar: tr|D3Y4D3|D3Y4D3_DANRE 11 264 . CM004438.1:subseq(17080000,150000) 4805 4052 - 933 M 95 285 G 1 0 M 6 18 G 1 0 M 150 450 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-23 ##type DNA # # # seqname source feature start end score strand frame attributes # CM004438.1:subseq(17080000,150000) exonerate:protein2genome:local gene 4053 4805 933 - . gene_id 1 ; sequence tr|D3Y4D3|D3Y4D3_DANRE ; gene_orientation . CM004438.1:subseq(17080000,150000) exonerate:protein2genome:local cds 4053 4805 . - . CM004438.1:subseq(17080000,150000) exonerate:protein2genome:local exon 4053 4805 . - . insertions 0 ; deletions 2 CM004438.1:subseq(17080000,150000) exonerate:protein2genome:local similarity 4053 4805 933 - . alignment_id 1 ; Query tr|D3Y4D3|D3Y4D3_DANRE ; Align 4806 12 285 ; Align 4521 108 18 ; Align 4503 115 450 # --- END OF GFF DUMP --- # -- completed exonerate analysis