Command line: [exonerate -m p2g --showtargetgff -q gpx3a.fa -t genomic.fa] Hostname: [ws101183u] C4 Alignment: ------------ Query: SPP00000617_2.0 # Protein # Glutathione peroxidase 3a (GPx3a) # Zebrafish Target: CM004427.1:subseq(4100000,1500000) Ictalurus punctatus breed USDA103 chromosome 14, whole genome shotgun sequence Model: protein2genome:local Raw score: 734 Query range: 25 -> 211 Target range: 37061 -> 39362 26 : AsnThrGlnAlaCysAsnSerAlaAlaGlyAspSerPheHisAsnTyrGlyAlaLysThr : 45 ::! !||| !||||||..!|||!.!! !! !|||.!.:!!|||||||||!.!|||||| SerGlnGlnLeuCysAsnGlnAlaValAspGlySerIleTyrAsnTyrGlyGlyLysThr 37062 : TCTCAACAGTTATGTAACCAGGCTGTGGATGGCTCCATATACAACTATGGAGGGAAGACG : 37119 46 : IleAsnGlyThrGlnPheIleProPheSerHisTyrAlaGlyLysHisValLeuValVal : 65 :!:!.!! ! ! ! !||||||||| !!!!.|||:!!||||||:!!|||||||||||| LeuThrGluGlyLeuHisIleProPheArgGlnTyrSerGlyLysTyrValLeuValVal 37120 : TTGACTGAAGGTTTGCACATTCCCTTCCGTCAGTATTCCGGGAAGTACGTCCTCGTTGTC : 37179 66 : AsnValAlaThrTyrUnkGlyLeuThrPheGlnTyrVal{G} >>>> Target Intr : 79 |||||||||:!!!:! ||||||||||||||||||:!!{|} 198 bp AsnValAlaSerPhe***GlyLeuThrPheGlnTyrLeu{G}++ 37180 : AACGTGGCCTCCTTCTGAGGACTTACCTTCCAGTATTTG{G}gt................ : 37224 80 : on 1 >>>> {lu}LeuAsnAlaLeuHisGluGluLeuArgHisLeuGlyPheThrIle : 94 {||}||||||||||||!!.|||||||||!:! !!||||||||||||||| ++{lu}LeuAsnAlaLeuGlnGluGluLeuLysAspLeuGlyPheThrIle 37225 : .........ag{AA}CTGAATGCACTACAGGAGGAGCTCAAAGACCTGGGCTTCACCATT : 37464 95 : LeuGlyPheProCysAspGlnPheGlyLysGlnGluProGlyGluAsnAsnGluIleLeu : 114 |||||||||||||||:!!|||||||||||||||||||||||| !|||.!.||||||||| LeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyLeuAsnGluGluIleLeu 37465 : CTCGGATTTCCGTGTAACCAGTTTGGGAAACAGGAACCCGGGTTAAATGAAGAGATCTTG : 37524 115 : SerAlaLeu{Ly} >>>> Target Intron 2 >>>> {s}TyrValArgProGly : 123 !! !|||{||} 1251 bp {|}||||||||||||||| ProLeuLeu{Ly}++ ++{s}TyrValArgProGly 37525 : CCACTGCTA{AA}gt.........................ag{G}TATGTCCGCCCAGGC : 38802 124 : AsnGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyAspGlyGlu : 143 |||||||||||||||||||||||||||||||||!..|||||||||||||||!!:..!||| AsnGlyPheValProAsnPheGlnLeuPheGluSerGlyAspValAsnGlyGluAsnGlu 38803 : AACGGATTTGTTCCAAACTTCCAGCTCTTTGAAAGCGGCGATGTGAACGGTGAAAACGAG : 38862 144 : GlnAlaLeuPheThrPheLeuLys >>>> Target Intron 3 >>>> AsnAlaC : 154 |||::!:!!||||||||||||||| 294 bp ||||||| GlnSerValPheThrPheLeuLys++ ++AsnAlaC 38863 : CAAAGTGTTTTCACTTTCCTGAAGgt.........................agAACGCCT : 39189 155 : ysProProValGlyGluSerPheGlyAlaThrSerAsnArgLeuPheTrpGluProLeuL : 174 ||||||||||||||.!.|||||||||::! !!!!!..!|||||||||||||||||||||| ysProProValGlyAsnSerPheGlySerProThrGlyArgLeuPheTrpGluProLeuL 39190 : GCCCACCTGTCGGCAACAGCTTCGGAAGCCCCACCGGCAGACTGTTTTGGGAGCCGCTGA : 39249 175 : ysValAsnAspIleLysTrpAsnPheGluLysPheLeuLeuAspProAspGlyArgProV : 194 |||||!.!||||||||||||||||||||||||||||||:!!! !|||||||||!:!|||: ysValThrAspIleLysTrpAsnPheGluLysPheLeuValGlyProAspGlyLysProL 39250 : AAGTCACCGACATCAAGTGGAACTTCGAGAAGTTTCTGGTGGGCCCTGATGGAAAGCCGC : 39309 195 : alMetArgTrpPheProArgValAsnValSerGluValArgAlaAspIleLeu : 211 !!||||||||||||||||||||||||||||||||||||||||||||||||:!! euMetArgTrpPheProArgValAsnValSerGluValArgAlaAspIleIle 39310 : TGATGAGGTGGTTTCCAAGGGTGAACGTTTCCGAAGTCCGGGCCGACATCATA : 39362 vulgar: SPP00000617_2.0 25 211 . CM004427.1:subseq(4100000,1500000) 37061 39362 + 734 M 53 159 S 0 1 5 0 2 I 0 194 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 1247 3 0 2 S 1 1 M 33 99 5 0 2 I 0 290 3 0 2 M 60 180 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-14 ##type DNA # # # seqname source feature start end score strand frame attributes # CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local gene 37062 39362 734 + . gene_id 1 ; sequence SPP00000617_2.0 ; gene_orientation + CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local cds 37062 37221 . + . CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local exon 37062 37221 . + . insertions 0 ; deletions 0 CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local splice5 37222 37223 . + . intron_id 1 ; splice_site "GT" CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local intron 37222 37419 . + . intron_id 1 CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local splice3 37418 37419 . + . intron_id 0 ; splice_site "AG" CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local cds 37420 37537 . + . CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local exon 37420 37537 . + . insertions 0 ; deletions 0 CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local splice5 37538 37539 . + . intron_id 2 ; splice_site "GT" CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local intron 37538 38788 . + . intron_id 2 CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local splice3 38787 38788 . + . intron_id 1 ; splice_site "AG" CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local cds 38789 38888 . + . CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local exon 38789 38888 . + . insertions 0 ; deletions 0 CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local splice5 38889 38890 . + . intron_id 3 ; splice_site "GT" CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local intron 38889 39182 . + . intron_id 3 CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local splice3 39181 39182 . + . intron_id 2 ; splice_site "AG" CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local cds 39183 39362 . + . CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local exon 39183 39362 . + . insertions 0 ; deletions 0 CM004427.1:subseq(4100000,1500000) exonerate:protein2genome:local similarity 37062 39362 734 + . alignment_id 1 ; Query SPP00000617_2.0 ; Align 37062 26 159 ; Align 37422 80 114 ; Align 38790 119 99 ; Align 39183 152 180 # --- END OF GFF DUMP --- # -- completed exonerate analysis