Command line: [exonerate -m p2g --showtargetgff -q sel15_zf.aa.fa.canviat -t genomesel15_JXSJ01000006.1.fa] Hostname: [ws116064U] C4 Alignment: ------------ Query: SPP00000632_2.0 # Protein # 15 kDa selenoprotein (Sel15) # Zebrafish Target: JXSJ01000006.1:subseq(2230275,103608) Miichthys miiuy scaffold6, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 607 Query range: 8 -> 146 Target range: 53617 -> 49999 9 : GlnLeuAlaSerTyrGlyAlaGluLeuSerSerGluAlaCysArgGluLeuGlyPheSer : 28 ||||||:!!:!!|||||||||!!:|||||||||||||||||||||||||||||||||||| GlnLeuSerAlaTyrGlyAlaAspLeuSerSerGluAlaCysArgGluLeuGlyPheSer 53617 : cagCTGTCGGCATACGGAGCTGATCTGTCGTCTGAGGCATGCAGGGAGCTGGGCTTCTCC : 53560 29 : SerAsnLeuLeuCysSerSerCysGluLeuLeuGlyGlnPheSerLeuAsnGlnLeuAsp : 48 ||||||||||||||||||||||||!!:|||||||||:!!|||||||||!.!:!!|||.!. SerAsnLeuLeuCysSerSerCysAspLeuLeuGlyGluPheSerLeuThrLysLeuGln 53559 : AGCAACCTGCTGTGTAGCTCCTGTGACCTGCTCGGGGAGTTCAGCCTCACCAAGCTCCAG : 53500 49 : LeuProCysArgGlnCysCysGlnGluGluAlaGlnLeuGluAsnArgLys >>>> Ta : 66 ! ! !||||||||||||||||||:!!|||||||||:!!||| |||||| ProGlyCysArgGlnCysCysGlnGlnGluAlaGlnMetGluAlaArgLys++ 53499 : CCCGGCTGCCGGCAGTGCTGTCAGCAGGAGGCCCAGATGGAAGCGCGCAAGgt....... : 53444 67 : rget Intron 1 >>>> LeuTyrProGlyAlaIleLeuGluValCysGlyUnkLysL : 79 2530 bp |||||| !!|||||||||||||||||||||||| |||| ++LeuTyrAlaGlyAlaIleLeuGluValCysGly***LysL 53443 : ..................agCTCTATGCTGGGGCCATCCTCGAGGTGTGTGGATGAAAAC : 50877 80 : euGlyArgPheProGlnValGln{A} >>>> Target Intron 2 >>>> {la}P : 88 |||||||||||||||||||||||{|} 496 bp {||}| euGlyArgPheProGlnValGln{A}++ ++{la}P 50876 : TGGGGAGGTTCCCCCAAGTCCAA{G}gt.........................ag{ct}t : 50354 89 : heValArgSerAspLysProLysLeuPheArgGlyLeuGlnIleLys >>>> Target : 104 |||||||||||||||||||||||:!!|||!:!||||||||||||||| 1 heValArgSerAspLysProLysMetPheLysGlyLeuGlnIleLys++ 50353 : TTGTCAGGAGTGACAAGCCAAAAATGTTCAAGGGTCTTCAGATCAAGgt........... : 50304 105 : Intron 3 >>>> TyrValArgGlySerAspProValLeuLysLeuLeuAspAspAs : 118 78 bp |||||||||||||||||||||||||||||||||||||||||||| ++TyrValArgGlySerAspProValLeuLysLeuLeuAspAspAs 50303 : ..............agTACGTTAGAGGCTCAGATCCTGTATTAAAGCTTTTGGACGATAA : 50086 119 : nGlyAsnIleAlaGluGluLeuSerIleLeuLysTrpAsnThrAspSerValGluGluPh : 138 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| nGlyAsnIleAlaGluGluLeuSerIleLeuLysTrpAsnThrAspSerValGluGluPh 50085 : CGGGAACATCGCTGAAGAGCTCAGCATCCTCAAGTGGAACACAGACAGCGTAGAGGAGTT : 50026 139 : eLeuSerGluLysLeuGluArgIle : 146 ||||||||||||||||!!:|||||| eLeuSerGluLysLeuAspArgIle 50025 : CCTGAGTGAGAAATTAGACCGGATA : 50000 vulgar: SPP00000632_2.0 8 146 . JXSJ01000006.1:subseq(2230275,103608) 53617 49999 - 607 M 57 171 5 0 2 I 0 2526 3 0 2 M 21 63 S 0 1 5 0 2 I 0 492 3 0 2 S 1 2 M 16 48 5 0 2 I 0 174 3 0 2 M 43 129 # --- 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 # JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local gene 50000 53617 607 - . gene_id 1 ; sequence SPP00000632_2.0 ; gene_orientation + JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local cds 53447 53617 . - . JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local exon 53447 53617 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local splice5 53445 53446 . - . intron_id 1 ; splice_site "GT" JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local intron 50917 53446 . - . intron_id 1 JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local splice3 50917 50918 . - . intron_id 0 ; splice_site "AG" JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local cds 50853 50916 . - . JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local exon 50853 50916 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local splice5 50851 50852 . - . intron_id 2 ; splice_site "GT" JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local intron 50357 50852 . - . intron_id 2 JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local splice3 50357 50358 . - . intron_id 1 ; splice_site "ag" JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local cds 50307 50356 . - . JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local exon 50307 50356 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local splice5 50305 50306 . - . intron_id 3 ; splice_site "GT" JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local intron 50129 50306 . - . intron_id 3 JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local splice3 50129 50130 . - . intron_id 2 ; splice_site "AG" JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local cds 50000 50128 . - . JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local exon 50000 50128 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(2230275,103608) exonerate:protein2genome:local similarity 50000 53617 607 - . alignment_id 1 ; Query SPP00000632_2.0 ; Align 53618 9 171 ; Align 50917 66 63 ; Align 50355 88 48 ; Align 50129 104 129 # --- END OF GFF DUMP --- # -- completed exonerate analysis