Command line: [exonerate -m p2g --showtargetgff -q selenou1a_zf.aa.fa.canviat -t genomeselenou1a_JXSJ01000036.1.fa] Hostname: [ws101073u] C4 Alignment: ------------ Query: SPP00000652_2.0 # Protein # Selenoprotein U1a (SELENOU1a) # Zebrafish Target: JXSJ01000036.1:subseq(1809825,53445) Miichthys miiuy scaffold36, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 825 Query range: 7 -> 222 Target range: 51945 -> 49999 8 : PheGluGlyGluThrLeuAlaMetGlyMetTrpSerLeuGlyLeuGlyAlaValGlyAla : 27 .!.||||||!!: !|||!.!||||||||||||||||||||||||||||||||||||||| LeuGluGlyAspLeuLeuGlyMetGlyMetTrpSerLeuGlyLeuGlyAlaValGlyAla 51945 : CTGGAGGGTGACCTGTTGGGGATGGGGATGTGGTCGTTGGGTCTGGGCGCTGTTGGCGCT : 51888 28 : AlaIleAlaGlyLeuIleLeuAlaAsnThrAspPheLeuLeuThrLysSerAlaProAla : 47 |||:!!||||||:!!.!!|||||||||||||||.!.! !|||||||||:!!||| !||| AlaLeuAlaGlyIlePheLeuAlaAsnThrAspLeuCysLeuThrLysAlaAlaAsnAla 51887 : GCGCTAGCCGGGATCTTCCTGGCCAACACCGATCTGTGTCTGACTAAAGCTGCCAACGCA : 51828 48 : ThrValAspTyrLeuAlaAsnAlaAspLeuLysThrIleAspGly{A} >>>> Targe : 63 |||:!!!!::!!|||! !:!!.!!|||||| :!!! ! ! !{|} ThrLeuGluHisLeuGluAspThrAspLeuCysSerThrThrAsp{A}++ 51827 : ACGCTGGAGCACCTTGAAGATACTGACCTGTGCTCCACCACAGAT{G}gt.......... : 51777 64 : t Intron 1 >>>> {sp}GluArgSerLeuLysAlaLysAlaLeuTrpGluLysSer : 76 564 bp {||}:!!!:! !:!!|||||||||::!|||||||||! !!:! ++{sp}GlnLysValIleLysAlaLysSerLeuTrpGluThrAsn 51776 : ...............ag{at}CAGAAGGTCATCAAAGCAAAGAGCCTATGGGAGACGAAT : 51177 77 : GlyAlaValIleMetAlaValArgArgProGlyUnkPheLeuCysArgGlu >>>> Ta : 94 |||||||||:!!||||||||||||||||||||| ||||||||||||||| GlyAlaValValMetAlaValArgArgProGly***PheLeuCysArgGlu++ 51176 : GGGGCTGTGGTCATGGCCGTACGACGGCCTGGATGATTTTTGTGCAGAGAGgt....... : 51121 95 : rget Intron 2 >>>> GluAlaSerGluLeuSerSerLeuLysProGlnLeuAspG : 107 89 bp ||||||||||||||||||||||||||||||||||||!!:| ++GluAlaSerGluLeuSerSerLeuLysProGlnLeuGluG 51120 : ..................agGAGGCCTCTGAGCTGTCCTCTCTGAAGCCCCAGCTGGAAG : 50995 108 : luLeuGlyValProLeuTyrAlaValValLysGluAsnValGlyThrGluIleGlnAspP : 127 ||||||||||||||||| !|||||||||||||||||||||||||||||||||||||||| luLeuGlyValProLeuValAlaValValLysGluAsnValGlyThrGluIleGlnAspP 50994 : AGCTCGGGGTCCCTCTGGTCGCCGTGGTGAAGGAGAACGTCGGCACAGAGATCCAGGACT : 50935 128 : heArgProHisPheAlaGlyGluIlePheLeuAspGluLys >>>> Target Intro : 141 |||||! !||||||||||||!!:|||!:!:!!||||||:!! 569 bp heArgLeuHisPheAlaGlyAspIleTyrValAspGluGln++ 50934 : TCAGACTGCACTTCGCTGGGGACATCTACGTAGATGAACAGgt................. : 50891 142 : n 3 >>>> GlnAlaPheTyrGlyProGlnGlnArgLysMetGlyGlyLeuGlyPheIl : 157 :!! !||||||||||||! !:!!|||||||||||||||||||||||||| ++LysHisPheTyrGlyProLeuGluArgLysMetGlyGlyLeuGlyPheIl 50890 : ........agAAACACTTCTACGGCCCACTGGAGAGGAAGATGGGCGGTCTGGGCTTCAT : 50276 158 : eArgLeuGlyValTrpGlnAsnPheValArgAlaTrpArgAlaGlyTyrGlnGlyAsnMe : 177 |||||||||||||||||||||||||:!!||||||||||||:!!||||||||||||||||| eArgLeuGlyValTrpGlnAsnPheMetArgAlaTrpArgSerGlyTyrGlnGlyAsnMe 50275 : TCGCCTCGGAGTGTGGCAGAACTTCATGCGGGCTTGGAGGTCCGGTTACCAGGGCAACAT : 50216 178 : tAsnGlyGluGlyPheIleLeuGlyGlyValPheValMetGlySerGlyGlyGln >>> : 196 |||||||||||||||||||||||||||||||||||||!!:|||:!!|||! !||| tAsnGlyGluGlyPheIleLeuGlyGlyValPheValIleGlyAlaGlyGluGln++ 50215 : GAACGGCGAGGGCTTCATCCTCGGAGGAGTGTTTGTCATCGGAGCAGGAGAGCAGgt... : 50157 197 : > Target Intron 4 >>>> GlyValLeuLeuGluHisArgGluLysGluPheGly : 207 79 bp |||:!!|||||||||||||||||||||||||||||| ++GlyIleLeuLeuGluHisArgGluLysGluPheGly 50156 : ......................agggGATTCTCCTGGAACATCGTGAGAAGGAGTTCGGC : 50047 208 : AspLysValSerLeuGluSerValLeuGluAlaAlaLysLysVal : 222 :!!||||||..! !! !..!||||||||||||!.!||||||:!! AsnLysValAspThrAlaAspValLeuGluAlaValLysLysIle 50046 : AATAAGGTGGATACTGCAGATGTTTTAGAGGCTGTTAAGAAAATT : 50000 vulgar: SPP00000652_2.0 7 222 . JXSJ01000036.1:subseq(1809825,53445) 51945 49999 - 825 M 55 165 S 0 1 5 0 2 I 0 560 3 0 2 S 1 2 M 30 90 5 0 2 I 0 85 3 0 2 M 47 141 5 0 2 I 0 565 3 0 2 M 55 165 5 0 2 I 0 75 3 0 2 M 27 81 # --- 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 # JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local gene 50000 51945 825 - . gene_id 1 ; sequence SPP00000652_2.0 ; gene_orientation + JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local cds 51780 51945 . - . JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local exon 51780 51945 . - . insertions 0 ; deletions 0 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice5 51778 51779 . - . intron_id 1 ; splice_site "GT" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local intron 51216 51779 . - . intron_id 1 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice3 51216 51217 . - . intron_id 0 ; splice_site "ag" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local cds 51124 51215 . - . JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local exon 51124 51215 . - . insertions 0 ; deletions 0 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice5 51122 51123 . - . intron_id 2 ; splice_site "GT" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local intron 51035 51123 . - . intron_id 2 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice3 51035 51036 . - . intron_id 1 ; splice_site "AG" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local cds 50894 51034 . - . JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local exon 50894 51034 . - . insertions 0 ; deletions 0 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice5 50892 50893 . - . intron_id 3 ; splice_site "GT" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local intron 50325 50893 . - . intron_id 3 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice3 50325 50326 . - . intron_id 2 ; splice_site "AG" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local cds 50160 50324 . - . JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local exon 50160 50324 . - . insertions 0 ; deletions 0 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice5 50158 50159 . - . intron_id 4 ; splice_site "GT" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local intron 50081 50159 . - . intron_id 4 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local splice3 50081 50082 . - . intron_id 3 ; splice_site "ag" JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local cds 50000 50080 . - . JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local exon 50000 50080 . - . insertions 0 ; deletions 0 JXSJ01000036.1:subseq(1809825,53445) exonerate:protein2genome:local similarity 50000 51945 825 - . alignment_id 1 ; Query SPP00000652_2.0 ; Align 51946 8 165 ; Align 51214 64 90 ; Align 51035 94 141 ; Align 50325 141 165 ; Align 50081 196 81 # --- END OF GFF DUMP --- # -- completed exonerate analysis