Command line: [exonerate -m p2g --showtargetgff -q selenot2_zf.aa.fa.canviat -t genomeselenot2_JXSJ01000199.1.fa] Hostname: [ws116066U] C4 Alignment: ------------ Query: SPP00000650_2.0 # Protein # Selenoprotein T2 (SELENOT2) # Zebrafish Target: JXSJ01000199.1:subseq(170221,101612) Miichthys miiuy scaffold199, whole genome shotgun sequence Model: protein2genome:local Raw score: 754 Query range: 0 -> 209 Target range: 49999 -> 51612 1 : MetAlaGluTyrSerGlnThrGlyIleLeuThrAlaLeuLeuLeuPheThrValValThr : 20 ||||||||||||||||||.!!|||:!!|||.!!||||||||||||||||||!.!:!!||| MetAlaGluTyrSerGlnAlaGlyLeuLeuAlaAlaLeuLeuLeuPheThrAlaLeuThr 50000 : ATGGCGGAGTACAGCCAGGCGGGCCTCCTGGCGGCTCTGCTCCTCTTCACGGCTCTGACG : 50057 21 : ValLysAspIleTyrValGlyArgAsnSer{Va} >>>> Target Intron 1 >>> : 31 |||:::|||:!!|||:!!||||||:::|||{:!} 33 bp ValArgAspValTyrLeuGlyArgSerSer{Le}+- 50058 : GTCCGTGACGTCTACCTCGGCAGATCGAGC{CT}gc........................ : 50094 32 : > {l}ThrGlnGlnGluAsnThrGlyProAspIleAsnThrGlnArgGlnAsnLysHis : 49 {!} !...! !!!:!.!..! !||||||! !.!. !! !!:!! !.!||| ! ++{u}LeuSerProAspThrValLeuProAspThrGluProGlyLysProThrLysSer 50095 : .ag{C}CTCTCCCCGGATACCGTACTACCGGACACAGAACCCGGAAAACCCACCAAGTCC : 50177 50 : ThrPheTyrThrGlyProValLeuLysPheGlnTyr{Cy} >>>> Target Intron : 62 :!!.!.||||||||||||||||||||||||||||||{||} 253 bp SerLeuTyrThrGlyProValLeuLysPheGlnTyr{Cy}++ 50178 : TCCCTGTACACCGGGCCTGTGCTCAAGTTCCAGTAC{TG}gt.................. : 50220 63 : 2 >>>> {s}IleSerUnkGlyTyrSerLysValPheGlnGluTyrSerArgSerIle : 78 {|}|||||| |||||||||||||||||||||||||||||||||:!!||| ++{s}IleSer***GlyTyrSerLysValPheGlnGluTyrSerArgAlaIle 50221 : .......ag{t}atCTCCTGAGGGTACAGCAAAGTGTTCCAGGAGTACTCACGAGCCATC : 50517 79 : SerGlnLeuTyrProAspIleArgIleGluGlyAspAsnTyrProProLysProIleAsn : 98 |||||||||||||||||||||||||||||||||!!:||||||||||||! |||.!!||| SerGlnLeuTyrProAspIleArgIleGluGlyGluAsnTyrProProThrProPheAsn 50518 : AGCCAGCTGTACCCGGACATCCGCATTGAGGGAGAAAACTACCCTCCAACCCCCTTCAAC : 50577 99 : {Ly} >>>> Target Intron 3 >>>> {s}TyrLeuGlyAsnPheLeuSerTyr : 107 {!:} 294 bp {!} :::||||||.!.:!!|||||| {Ar}++ ++{g}CysValGlyAsnLeuValSerTyr 50578 : {AG}gt.........................ag{g}tgtgttgGGAACCTGGTCTCCTAC : 50898 108 : PheLysLeuLeuAlaIleAlaLeuIleValThrGlyGlnAsnProPheGlnMetPheGly : 127 .!.|||||||||:::||| |||||||||!:!||||||||||||||| :!!.!.||| LeuLysLeuLeuSerIleLeuLeuIleValSerGlyGlnAsnProPheIleLeuLeuGly 50899 : CTGAAGctgctctccatcctcctcatcGTCAGCGGTCAGAACCCGTTCATCCTGCTGGGC : 50958 128 : MetAsnThrProArgIleTrpAlaTrpGlyGlnGluAsnLys >>>> Target Intr : 142 :!::!!||||||||| !|||.!!|||.!!|||!!:|||||| 73 bp LeuAspThrProArgAlaTrpThrTrpSerGlnAspAsnLys++ 50959 : CTCGACACTCCGCGGGCCTGGACCTGGAGTCAGGACAACAAGgt................ : 51005 143 : on 4 >>>> IlePheSerCysLeuMetAlaPhePheIleSerAsnMetLeuGluThrH : 158 |||||||||||||||||||||||||||:!!!!!||||||:!!||||||| ++IlePheSerCysLeuMetAlaPhePheLeuCysAsnMetMetGluThrH 51006 : .........agATCTTCTCCTGTCTGATGGCCTTCTTCCTCTGCAACATGATGGAGACTC : 51124 159 : isPheLeuSerThrGlyAlaPheGluIleThrLeuAsn{A} >>>> Target Intro : 171 ||||||||||||||||||||||||||:!!|||||||||{|} 333 bp isPheLeuSerThrGlyAlaPheGluValThrLeuAsn{A}++ 51125 : ACTTCCTGTCTACGGGAGCCTTCGAGGTCACTCTGAAC{G}gt................. : 51166 172 : n 5 >>>> {sp}IleProIleTrpSerLysLeuGlnSerGlyTyrValProAsnIleG : 187 {||}:!!|||:!!||||||||||||||||||||||||||||||||||||| ++{sp}ValProLeuTrpSerLysLeuGlnSerGlyTyrValProAsnIleG 51167 : ........ag{ac}gTCCCTCTATGGTCGAAGCTTCAGTCGGGTTACGTCCCGAACATCC : 51544 188 : lnGluLeuPheGlnIleLeuAspAsnHisLeuLysMetAsnGlnAlaAspLysMetAsnP : 207 |||||:!!||||||||||||||||||||||||||||||||||||!.!||| |||!:!| lnGluIlePheGlnIleLeuAspAsnHisLeuLysMetAsnGlnValAspProMetSerP 51545 : AGGAGATCTTCCAGATCCTCGACAACCACCTGAAGATGAACCAGGTGGACCCCATGAGCT : 51604 208 : heProSer : 209 || !!||| heThrSer 51605 : TCACCTCC : 51612 vulgar: SPP00000650_2.0 0 209 . JXSJ01000199.1:subseq(170221,101612) 49999 51612 + 754 M 30 90 S 0 2 5 0 2 I 0 29 3 0 2 S 1 1 M 30 90 S 0 2 5 0 2 I 0 249 3 0 2 S 1 1 M 36 108 S 0 2 5 0 2 I 0 290 3 0 2 S 1 1 M 42 126 5 0 2 I 0 69 3 0 2 M 29 87 S 0 1 5 0 2 I 0 329 3 0 2 S 1 2 M 38 114 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local gene 50000 51612 754 + . gene_id 1 ; sequence SPP00000650_2.0 ; gene_orientation + JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local cds 50000 50091 . + . JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local exon 50000 50091 . + . insertions 0 ; deletions 0 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice5 50092 50093 . + . intron_id 1 ; splice_site "GC" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local intron 50092 50124 . + . intron_id 1 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice3 50123 50124 . + . intron_id 0 ; splice_site "AG" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local cds 50125 50217 . + . JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local exon 50125 50217 . + . insertions 0 ; deletions 0 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice5 50218 50219 . + . intron_id 2 ; splice_site "GT" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local intron 50218 50470 . + . intron_id 2 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice3 50469 50470 . + . intron_id 1 ; splice_site "ag" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local cds 50471 50581 . + . JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local exon 50471 50581 . + . insertions 0 ; deletions 0 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice5 50582 50583 . + . intron_id 3 ; splice_site "GT" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local intron 50582 50875 . + . intron_id 3 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice3 50874 50875 . + . intron_id 2 ; splice_site "ag" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local cds 50876 51002 . + . JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local exon 50876 51002 . + . insertions 0 ; deletions 0 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice5 51003 51004 . + . intron_id 4 ; splice_site "GT" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local intron 51003 51075 . + . intron_id 4 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice3 51074 51075 . + . intron_id 3 ; splice_site "AG" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local cds 51076 51163 . + . JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local exon 51076 51163 . + . insertions 0 ; deletions 0 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice5 51164 51165 . + . intron_id 5 ; splice_site "GT" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local intron 51164 51496 . + . intron_id 5 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local splice3 51495 51496 . + . intron_id 4 ; splice_site "ag" JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local cds 51497 51612 . + . JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local exon 51497 51612 . + . insertions 0 ; deletions 0 JXSJ01000199.1:subseq(170221,101612) exonerate:protein2genome:local similarity 50000 51612 754 + . alignment_id 1 ; Query SPP00000650_2.0 ; Align 50000 1 90 ; Align 50126 32 90 ; Align 50472 63 108 ; Align 50877 100 126 ; Align 51076 142 87 ; Align 51499 172 114 # --- END OF GFF DUMP --- # -- completed exonerate analysis