Command line: [exonerate -m p2g --showtargetgff -q gpx5_hs.aa.fa.canviat -t genomegpx5_JXSJ01000010.1.fa] Hostname: [ws101074U] C4 Alignment: ------------ Query: SPP00000036_2.0 # Protein # Glutathione peroxidase 5 (GPx5) # Human Target: JXSJ01000010.1:subseq(2769482,101618) Miichthys miiuy scaffold10, whole genome shotgun sequence Model: protein2genome:local Raw score: 585 Query range: 29 -> 221 Target range: 49993 -> 51618 30 : MetAspCysHisLysAspGluLysGlyThrIleTyrAspTyrGluAlaIleAlaLeuAsn : 49 ..!.!.|||.!...!! ! ! ! ||||||||||||.!.|||...|||! .!!|||||| GlnGlnCysGluSerAlaThrAspGlyThrIleTyrGlnTyrSerAlaLysThrLeuAsn 49994 : CAGCAGTGTGAGTCAGCCACAGATGGCACCATTTACCAATACAGTGCCAAGACTCTAAAC : 50051 50 : LysAsnGluTyrValSerPheLysGlnTyrValGlyLysHisIleLeuPheValAsnVal : 69 !!:!... !|||!:!|||!...!.|||..!|||:!! !:!!||||||||||||||| GlySerArgThrValAsnPheSerAspTyrThrGlyGluSerValLeuPheValAsnVal 50052 : GGGAGTCGCACTGTGAACTTCAGTGACTACACGGGCGAGAGTGTCCTCTTCGTCAACGTG : 50111 70 : AlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G} >>>> Target Intron 1 > : 81 |||||||||!! |||! !||| !|||||| !{|} 92 bp AlaThrTyr***GlyTyrThrPheGlnTyrVal{G}++ 50112 : GCCACATACTGAGGATACACCTTCCAGTATGTG{G}gt...................... : 50150 82 : >>> {lu}LeuAsnAlaLeuGlnGluGluLeuLysProTyrGlyLeuValValLeuGly : 98 {||}||||||||||||!!.||||||:!!|||||| !||||||..!:!!|||||| ++{lu}LeuAsnAlaLeuHisGluGluMetLysProLeuGlyLeuThrIleLeuGly 50151 : ...ag{AA}CTGAATGCACTACATGAGGAAATGAAACCTCTTGGACTCACCATCCTCGGC : 50290 99 : PheProCysAsnGlnPheGlyLysGlnGluProGlyAspAsnLysGluIleLeuProGly : 118 ||||||||||||||||||||||||||||||||||||.!.!!. ! ||||||||||||||| PheProCysAsnGlnPheGlyLysGlnGluProGlyGlnLysHisGluIleLeuProGly 50291 : TTTCCCTGCAACCAATTTGGGAAACAGGAACCAGGGCAGAAACACGAAATTCTGCCAGGT : 50350 119 : Leu{Ly} >>>> Target Intron 2 >>>> {s}TyrValArgProGlyGlyGly : 127 |||{||} 547 bp {|}:!!||||||||||||..!||| Leu{Ly}++ ++{s}HisValArgProGlyAsnGly 50351 : TTA{AA}gt.........................ag{G}CATGTCAGACCAGGCAATGGA : 50924 128 : PheValProSerPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLys : 147 |||||||||!:!|||! !||||||!!:||||||||||||||||||:!! ! ||||||:!! PheValProAsnPheLeuLeuPheAspLysGlyAspValAsnGlyLysAspGluGlnGlu 50925 : TTTGTTCCAAACTTCCTGCTGTTTGACAAAGGTGATGTGAATGGAAAAGACGAGCAAGAG : 50984 148 : ValPheSerPheLeuLys >>>> Target Intron 3 >>>> HisSerCysProH : 158 ||||||!!!||||||||| 407 bp :!!|||||||||! ValPheThrPheLeuLys++ ++AsnSerCysProP 50985 : GTTTTCACTTTTCTCAAGgt.........................agAACTCCTGCCCTC : 51424 159 : isProSerGluIleLeuGlyThr<->PheLysSerIleSerTrpAspProValLysValH : 177 !.!!!!:! !!!!|||!.! !!!.!!!:!:! !|||.!.|||:!!|||:!! roValGlyAspAsnPheGlyAsnProThrAsnArgLeuPheTrpGlnProLeuLysMetS 51425 : CAGTCGGAGACAACTTCGGCAACCCCACCAACAGATTGTTCTGGCAGCCTCTGAAGATGA : 51484 178 : isAspIleArgTrpAsnPheGluLysPheLeuValGlyProAspGlyIleProValMetA : 197 !||||||!:!||||||||||||||||||||||||||||||||||||! |||||||||| erAspIleLysTrpAsnPheGluLysPheLeuValGlyProAspGlyLysProValMetA 51485 : GCGACATCAAGTGGAACTTCGAAAAGTTTCTGGTGGGCCCAGACGGGAAGCCAGTGATGA : 51544 198 : rgTrpSerHisArgAlaThrValSerSerValLysThrAspIleLeuAlaTyrLeuLysG : 217 ||||| !! !!!!.!!.!:!!|||..!|||::!:!!||||||! ! !|||||| : rgTrpHisProSerValAsnIleSerGluValArgSerAspIleArgLysTyrLeuLeuL 51545 : GGTGGCACCCAAGTGTCAACATTTCTGAAGTCCGATCTGACATCCGCAAATACCTGCTCA : 51604 218 : lnPheLysThrLys : 221 !!.!. ! |||||| ysLeuTyrThrLys 51605 : AACTATACACAAAG : 51618 vulgar: SPP00000036_2.0 29 221 . JXSJ01000010.1:subseq(2769482,101618) 49993 51618 + 585 M 51 153 S 0 1 5 0 2 I 0 88 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 543 3 0 2 S 1 1 M 33 99 5 0 2 I 0 403 3 0 2 M 12 36 G 0 3 M 56 168 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-12-01 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local gene 49994 51618 585 + . gene_id 1 ; sequence SPP00000036_2.0 ; gene_orientation + JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local cds 49994 50147 . + . JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local exon 49994 50147 . + . insertions 0 ; deletions 0 JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local splice5 50148 50149 . + . intron_id 1 ; splice_site "GT" JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local intron 50148 50239 . + . intron_id 1 JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local splice3 50238 50239 . + . intron_id 0 ; splice_site "AG" JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local cds 50240 50357 . + . JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local exon 50240 50357 . + . insertions 0 ; deletions 0 JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local splice5 50358 50359 . + . intron_id 2 ; splice_site "GT" JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local intron 50358 50904 . + . intron_id 2 JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local splice3 50903 50904 . + . intron_id 1 ; splice_site "AG" JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local cds 50905 51004 . + . JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local exon 50905 51004 . + . insertions 0 ; deletions 0 JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local splice5 51005 51006 . + . intron_id 3 ; splice_site "GT" JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local intron 51005 51411 . + . intron_id 3 JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local splice3 51410 51411 . + . intron_id 2 ; splice_site "AG" JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local cds 51412 51618 . + . JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local exon 51412 51618 . + . insertions 3 ; deletions 0 JXSJ01000010.1:subseq(2769482,101618) exonerate:protein2genome:local similarity 49994 51618 585 + . alignment_id 1 ; Query SPP00000036_2.0 ; Align 49994 30 153 ; Align 50242 82 114 ; Align 50906 121 99 ; Align 51412 154 36 ; Align 51451 166 168 # --- END OF GFF DUMP --- # -- completed exonerate analysis