Command line: [exonerate -m p2g --showtargetgff -q selenou3_zf.aa.fa.canviat -t genomeselenou3_JXSJ01000043.1.fa] Hostname: [ws101073u] C4 Alignment: ------------ Query: SPP00000656_2.0 # Protein # Selenoprotein U3 (SELENOU3) # Zebrafish Target: JXSJ01000043.1:subseq(24604,103059) Miichthys miiuy scaffold43, whole genome shotgun sequence Model: protein2genome:local Raw score: 705 Query range: 0 -> 179 Target range: 49999 -> 58267 1 : ValGluIleGlySerLeuTrpArgGluGlnAlaValValLeuPhePheLeuArgArgPhe : 20 ||||||:!! !|||||||||!:!!!:||| !!||||||||||||||||||||||||||| ValGluLeuGlnSerLeuTrpGlnAspGlnProValValLeuPhePheLeuArgArgPhe 50000 : GTAGAGCTTCAGTCTCTGTGGCAGGACCAGCCGGTGGTATTGTTCTTCCTGCGCAGATTT : 50057 21 : GlyCysGlnValCysArgTrpMetAlaAlaGluValSerLysLeuGluLysAspLeuLys : 40 |||||||||:!!|||||||||||||||:!!|||:!!|||||||||||| !||||||!:! GlyCysGlnIleCysArgTrpMetAlaSerGluIleSerLysLeuGluProAspLeuArg 50058 : GGCTGTCAGATATGTCGGTGGATGGCATCAGAGATCAGCAAACTAGAGCCGGACCTGAGA : 50117 41 : AlaHisGlyIleAlaLeuValGlyIleGlyProGluGluThrGlyValLysGluPheLys : 60 ||| ! !!:!!.!!|||||||||:!!|||||||||||| !|||:!!|||||||||||| AlaSerArgValThrLeuValGlyValGlyProGluGluPheGlyLeuLysGluPheLys 50118 : GCCAGCCGTGTAACGCTGGTGGGCGTCGGACCAGAAGAGTTTGGGCTGAAGGAGTTCAAG : 50177 61 : AspGlyGlyPhePheLysGly{A} >>>> Target Intron 1 >>>> {sp}Ile : 69 !!:||||||||||||||||||{.} 685 bp {.!}||| GluGlyGlyPhePheLysGly{S}++ ++{er}Ile 50178 : GAAGGAGGGTTTTTCAAAGGA{T}gt.........................ag{CC}ATT : 50889 70 : TyrIleAspGluMetLysGlnCysTyrLysAspLeuGlyPheLys{Ar} >>>> Targ : 85 ||||||||||||...|||:!! !|||||||||||||||||||||{||} TyrIleAspGluGlnLysLysThrTyrLysAspLeuGlyPheLys{Ar}++ 50890 : TACATTGACGAACAAAAGAAAACCTACAAGGACTTGGGCTTCAAA{AG}gt......... : 50941 86 : et Intron 2 >>>> {g}TyrAsnAlaIleAsnValValProAlaAlaMetGlyLys : 98 365 bp {|}|||!..|||:!!!:!|||||||||||||||:!!|||||| ++{g}TyrThrAlaValSerValValProAlaAlaLeuGlyLys 50942 : ................ag{A}TACACAGCTGTAAGTGTGGTTCCTGCTGCTTTGGGGAAG : 51341 99 : LysValArgGluIleAlaSerLys >>>> Target Intron 3 >>>> AlaSerA : 109 ||||||||||||||||||:!!||| 140 bp |||!.!| LysValArgGluIleAlaAlaLys++ ++AlaLysA 51342 : AAAGTACGAGAGATAGCTGCAAAGgt.........................agGCCAAAG : 51514 110 : laGluGlyIleGlnGlyAsnPheSerGlyAspLeuLeuGlnSerGlyGlyMetLeuIleV : 129 ||!!:||||||||||||||||||||||||||||||||||||||||||||||||||||||| laAspGlyIleGlnGlyAsnPheSerGlyAspLeuLeuGlnSerGlyGlyMetLeuIleV 51515 : CTGATGGCATCCAGGGAAACTTCTCAGGAGATCTGCTGCAGAGTGGAGGCATGCTTATTG : 51574 130 : alAlaLys{G} >>>> Target Intron 4 >>>> {ly}GlyGluLysValLeuL : 138 ||||||||{|} 1360 bp {||}|||!!:|||||||||| alAlaLys{G}++ ++{ly}GlyAspLysValLeuL 51575 : TGGCCAAA{G}gt.........................ag{GT}GGAGATAAGGTCCTAC : 52961 139 : euHisPheIleGlnLysSerProAlaAspAsnProProLeuGluGluIleThrLysAlaL : 158 |||||||||||||| ! ||||||!.!|||:!! !|||||||||!!:|||:!!||||||| euHisPheIleGlnAspSerProGlyAspHisLeuProLeuGluAspIleSerLysAlaL 52962 : TGCACTTTATCCAGGATTCACCGGGAGACCACTTACCTCTGGAGGATATTTCCAAGGCTT : 53021 159 : euGlyIleSerAlaAsnValGlnAlaGlyGluLysProGln >>>> Target Intro : 172 ||||||||||||||!..|||:!!||||||:!!!:!||| ! 5181 bp euGlyIleSerAlaThrValLysAlaGlyGlnArgProVal++ 53022 : TGGGTATCTCAGCCACAGTAAAAGCAGGACAGAGACCAGTGgt................. : 53065 173 : n 5 >>>> CysAsnGluAspValCysThrArg : 179 ||||||!!:||||||||||||||| ++CysAsnAspAspValCysThrArg 53066 : ........agTGCAACGATGATGTTTGTACTCGA : 58267 vulgar: SPP00000656_2.0 0 179 . JXSJ01000043.1:subseq(24604,103059) 49999 58267 + 705 M 67 201 S 0 1 5 0 2 I 0 681 3 0 2 S 1 2 M 16 48 S 0 2 5 0 2 I 0 361 3 0 2 S 1 1 M 21 63 5 0 2 I 0 136 3 0 2 M 25 75 S 0 1 5 0 2 I 0 1356 3 0 2 S 1 2 M 39 117 5 0 2 I 0 5177 3 0 2 M 8 24 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-21 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local gene 50000 58267 705 + . gene_id 1 ; sequence SPP00000656_2.0 ; gene_orientation + JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local cds 50000 50201 . + . JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local exon 50000 50201 . + . insertions 0 ; deletions 0 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice5 50202 50203 . + . intron_id 1 ; splice_site "GT" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local intron 50202 50886 . + . intron_id 1 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice3 50885 50886 . + . intron_id 0 ; splice_site "AG" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local cds 50887 50938 . + . JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local exon 50887 50938 . + . insertions 0 ; deletions 0 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice5 50939 50940 . + . intron_id 2 ; splice_site "GT" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local intron 50939 51303 . + . intron_id 2 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice3 51302 51303 . + . intron_id 1 ; splice_site "AG" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local cds 51304 51367 . + . JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local exon 51304 51367 . + . insertions 0 ; deletions 0 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice5 51368 51369 . + . intron_id 3 ; splice_site "GT" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local intron 51368 51507 . + . intron_id 3 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice3 51506 51507 . + . intron_id 2 ; splice_site "aG" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local cds 51508 51583 . + . JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local exon 51508 51583 . + . insertions 0 ; deletions 0 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice5 51584 51585 . + . intron_id 4 ; splice_site "GT" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local intron 51584 52943 . + . intron_id 4 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice3 52942 52943 . + . intron_id 3 ; splice_site "AG" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local cds 52944 53062 . + . JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local exon 52944 53062 . + . insertions 0 ; deletions 0 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice5 53063 53064 . + . intron_id 5 ; splice_site "GT" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local intron 53063 58243 . + . intron_id 5 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local splice3 58242 58243 . + . intron_id 4 ; splice_site "AG" JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local cds 58244 58267 . + . JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local exon 58244 58267 . + . insertions 0 ; deletions 0 JXSJ01000043.1:subseq(24604,103059) exonerate:protein2genome:local similarity 50000 58267 705 + . alignment_id 1 ; Query SPP00000656_2.0 ; Align 50000 1 201 ; Align 50889 69 48 ; Align 51305 86 63 ; Align 51508 107 75 ; Align 52946 133 117 ; Align 58244 172 24 # --- END OF GFF DUMP --- # -- completed exonerate analysis