Command line: [exonerate -m p2g --showtargetgff -q selenol_zf.aa.fa.canviat -t genomeselenol_JXSJ01000513.1.fa] Hostname: [ws116064U] C4 Alignment: ------------ Query: SPP00000637_2.0 # Protein # Selenoprotein L (SELENOL) # Zebrafish Target: JXSJ01000513.1:subseq(100740,100567) Miichthys miiuy scaffold513, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 428 Query range: 64 -> 294 Target range: 55712 -> 50002 65 : AsnGlnLysGluAlaGluSerLeuTrpGluArgPhePheGlnAsnThrGluValArgGlu : 84 .!..!.!!.:!!!.!|||.!!:!!!! ! ||| ! !|||.!. !!|||!.! !!||| GluAsnAsnGlnGlyGluGlyVal***AlaArgHisHisGlnGluProGluAla***Glu 55712 : GAAAACAATCAAGGAGAAGGAGTTTGAGCTCGTCACCATCAAGAACCAGAAGCTTGAGAA : 55655 85 : SerValGluGlyLeuLeuGlnLeuGlu >>>> Target Intron 1 >>>> ValG : 95 ||||||:!!|||||| ::: ||| 4147 bp |||: SerValGlnGlyLeuAlaArgGlyGlu+- ++ValG 55654 : TCTGTGCAGGGCCTtgcaagaggagagga.........................agGTGC : 51475 96 : luTrpAspArgPheLeuLysArgLeuAspValGlnMetGlnMetSerAspThrValLeuS : 115 !!|||||| !||||||:!!|||||||||! !:!!:!:|||:!!|||! !! !:!!|||. lnTrpAspAlaPheLeuGlnArgLeuAspAspGluLeuGlnLeuSerAlaLysIleLeuA 51474 : AGTGGGATGCCTTTCTGCAGCGTTTGGATGATGAGTTACAACTGAGTGCTAAAATACTGG : 51415 116 : erGlnHisProAlaAlaGlnThrLeuSerSerAspMetArgPheLeuHisValGlnThrL : 135 .! ! !!! ! !:!!:!! !:!!||| !!|||! !.!!:!!:!!!.!::!||| spGlyAspHisGlnSerLysArgValSerProAspThrSerLeuIleAsnAlaArgThrG 51414 : ATGGGGACCATCAGTCAAAACGTGTTTCACCTGATACATCGCTCATCAACGCCAGAACAG : 51355 136 : ys{Gl} >>>> Target Intron 2 >>>> {u}AsnValSerLeuGlyGlnTyrL : 144 !{||} 82 bp {|}!..|||!!!!!!! !:!!|||| ly{Gl}++ ++{u}ThrValThrPheGluLysTyrL 51354 : GA{GA}gt.........................ag{A}ACAGTGACATTCGAGAAATATC : 51246 145 : euGlyLysGlyGluAsnLeuLeuLeuValLeuLeuArgHisPheGlyUnkLeuProUnkA : 164 |||||||||||:!!!!.:!!||||||||||||:!!|||!!.|||..! |||! ! | euGlyLysGlyLysLysIleLeuLeuValLeuIleArgGlnPheSerCysLeuLeuCysA 51245 : TCGGGAAAGGTAAAAAGATCCTGCTGGTGTTGATCCGTCAGTTCTCTTGTCTACTCTGCC : 51186 165 : rgAspHisLeuThrGluLeuLysAsnSerLys >>>> Target Intron 3 >>>> : 175 || !||||||! !!!:|||:!!!!.!.!:!! 580 bp rgLeuHisLeuLysAspLeuGluLysLysGln++ + 51185 : GCCTCCACCTTAAAGATCTGGAGAAGAAACAGgt.........................a : 50575 176 : MetLeuLeuGluAlaGlnSerLeuArgValLeuValValSerTyrGlySerLeuGluGl : 194 ! ! !|||!!:|||!!.|||:!!..!|||:!!|||||||||!:!|||!!!! !||||| +SerSerLeuAspAlaHisSerValGluValValValValSerPheGlyCysGlnGluGl 50574 : gAGCTCCCTCGATGCGCACTCAGTTGAAGTGGTGGTCGTTTCGTTCGGCTGTCAGGAAGG : 50516 195 : yAlaThrPheTrpLeuGluGlnThrGlyPheGluPheAspMetLeuLeuAspThrGluAr : 214 ||||:!! !||||||:!!.!.|||! !! !:!!!:!||||||||||||||| !!!!:|| yAlaSerHisTrpLeuGlnAspThrAspCysGlnTyrAspMetLeuLeuAspProAspAr 50515 : AGCTTCACACTGGCTGCAGGACACCGACTGTCAATACGACATGCTGTTGGACCCTGATAG : 50456 215 : gThr >>>> Target Intron 4 >>>> ValTyrLysMetPheGlyLeuGlySer : 224 |! ! 81 bp :!!||| ! ||||||||||||:!! gLys++ ++MetTyrAlaAlaPheGlyLeuGlyAla 50455 : GAAGgt.........................agATGTACGCAGCATTCGGCCTCGGAGCG : 50345 225 : SerMetSerLysValMetLysPheLysLeuMetPheHisTyrSerGluIleMetAlaMet : 244 |||:!!!.!||||||:!:!!.|||!:! !|||.!.! !|||:!!||| !:!!||| SerLeuLysLysValLeuAsnPheArgAsnMetLeuLeuTyrAlaGluTyrValAlaAsp 50344 : TCCCTGAAGAAAGTGTTAAACTTCAGAAACATGCTGCTCTACGCTGAATACGTGGCCGAC : 50285 245 : AsnArgThrLeuProGluMetProProGlnPheIleGluAspLeuPheGln >>>> Ta : 262 |||! ! !!!!|||..! !! !|||..!.!! !!!:|||:!!|||||| AsnMetGluPheProArgGlyLeuProSerIleGlnAspAspMetPheGln++ 50284 : AATATGGAGTTTCCACGAGGGCTTCCATCGATTCAAGACGACATGTTTCAGgt....... : 50229 263 : rget Intron 5 >>>> MetGlyGlyAspPheValLeuGluGlnAspGlyLysValI : 275 130 bp :::||||||||||||||||||!!::!! !!|||!:!|||: ++LeuGlyGlyAspPheValLeuAspGluHisGlyArgValL 50228 : ..................agctggGCGGCGACTTTGTTCTGGATGAACATGGGAGGGTGC : 50062 276 : lePheSerTyrArgCysLysSerProValAspArgProSerAlaThrGlnIleLeuAla : 294 !:||||||:!! !!|||:!!||||||:!!||||||||||||!.! !.!.||||||:!! euPheSerHisCysCysGlnSerProIleAspArgProSerValGluAspIleLeuSer 50061 : TGTTCTCTCACTGCTGTCAGAGTCCCATAGACAGACCCTCAGTGGAGGACATTTTATCT : 50003 vulgar: SPP00000637_2.0 64 294 . JXSJ01000513.1:subseq(100740,100567) 55712 50002 - 428 M 29 87 5 0 2 I 0 4143 3 0 2 M 42 126 S 0 2 5 0 2 I 0 78 3 0 2 S 1 1 M 38 114 5 0 2 I 0 576 3 0 2 M 41 123 5 0 2 I 0 77 3 0 2 M 46 138 5 0 2 I 0 126 3 0 2 M 33 99 # --- 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 # JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local gene 50003 55712 428 - . gene_id 1 ; sequence SPP00000637_2.0 ; gene_orientation + JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local cds 55626 55712 . - . JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local exon 55626 55712 . - . insertions 0 ; deletions 0 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice5 55624 55625 . - . intron_id 1 ; splice_site "ga" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local intron 51479 55625 . - . intron_id 1 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice3 51479 51480 . - . intron_id 0 ; splice_site "AG" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local cds 51351 51478 . - . JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local exon 51351 51478 . - . insertions 0 ; deletions 0 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice5 51349 51350 . - . intron_id 2 ; splice_site "GT" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local intron 51269 51350 . - . intron_id 2 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice3 51269 51270 . - . intron_id 1 ; splice_site "ag" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local cds 51154 51268 . - . JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local exon 51154 51268 . - . insertions 0 ; deletions 0 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice5 51152 51153 . - . intron_id 3 ; splice_site "GT" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local intron 50574 51153 . - . intron_id 3 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice3 50574 50575 . - . intron_id 2 ; splice_site "AG" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local cds 50451 50573 . - . JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local exon 50451 50573 . - . insertions 0 ; deletions 0 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice5 50449 50450 . - . intron_id 4 ; splice_site "GT" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local intron 50370 50450 . - . intron_id 4 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice3 50370 50371 . - . intron_id 3 ; splice_site "AG" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local cds 50232 50369 . - . JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local exon 50232 50369 . - . insertions 0 ; deletions 0 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice5 50230 50231 . - . intron_id 5 ; splice_site "GT" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local intron 50102 50231 . - . intron_id 5 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local splice3 50102 50103 . - . intron_id 4 ; splice_site "ag" JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local cds 50003 50101 . - . JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local exon 50003 50101 . - . insertions 0 ; deletions 0 JXSJ01000513.1:subseq(100740,100567) exonerate:protein2genome:local similarity 50003 55712 428 - . alignment_id 1 ; Query SPP00000637_2.0 ; Align 55713 65 87 ; Align 51479 94 126 ; Align 51268 137 114 ; Align 50574 175 123 ; Align 50370 216 138 ; Align 50102 262 99 # --- END OF GFF DUMP --- # -- completed exonerate analysis