Command line: [exonerate -m p2g --showtargetgff -q Sel15.x.fa -t Sel15.extraction.KV884860.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: Sel15 # Protein # 15 kDa selenoprotein (Sel15) # Zebrafish Target: KV884860.1:subseq(2942292,100000) Monopterus albus unplaced genomic scaffold scaffold172.1, whole genome shotgun sequence Model: protein2genome:local Raw score: 603 Query range: 0 -> 146 Target range: 46844 -> 50137 1 : LeuPheLeuLeuValPhePheAlaGlnLeuAlaSerTyrGlyAlaGluLeuSerSerGlu : 20 |||.!. |||::: |||||||||:!!||||||!.!||||||||||||||| LeuLeuProLeuLeuSerSerLeuGlnLeuAlaAlaTyrGlyGlyGluLeuSerSerGlu 46845 : CTTCTacctcttctttcctctctccagcTGGCGGCATATGGTGGTGAGCTGTCGTCAGAG : 46902 21 : AlaCysArgGluLeuGlyPheSerSerAsnLeuLeuCysSerSerCysGluLeuLeuGly : 40 ||||||||||||||||||||||||||||||||||||||||||||||||!!:||||||||| AlaCysArgGluLeuGlyPheSerSerAsnLeuLeuCysSerSerCysAspLeuLeuGly 46903 : GCGTGCAGGGAGCTGGGCTTCTCCAGCAACCTGCTGTGCAGCTCCTGTGACCTGCTCGGG : 46962 41 : GlnPheSerLeuAsnGlnLeuAspLeuProCysArgGlnCysCysGlnGluGluAlaGln : 60 :!!|||||||||!.!:!!|||.!.! ! !||||||||||||||||||:!!||||||||| GluPheSerLeuThrLysLeuGlnProAsnCysArgGlnCysCysGlnGlnGluAlaGln 46963 : GAGTTCAGCCTCACCAAGCTCCAGCCCAACTGCCGGCAGTGCTGCCAGCAGGAGGCCCAG : 47022 61 : LeuGluAsnArgLys >>>> Target Intron 1 >>>> LeuTyrProGlyAlaI : 71 :!!||| |||||| 2189 bp |||||| !!||||||| MetGluAlaArgLys++ ++LeuTyrAlaGlyAlaI 47023 : ATGGAAGCACGCAAGgt.........................agCTCTATGCTGGAGCCA : 49244 72 : leLeuGluValCysGlyUnkLysLeuGlyArgPheProGlnValGln{A} >>>> Tar : 87 ||||||||||||||||| |||||||||||||||||||||||||||{|} leLeuGluValCysGly***LysLeuGlyArgPheProGlnValGln{A}++ 49245 : TCCTCGAGGTATGTGGATGAAAACTGGGGAGGTTCCCTCAAGTCCAA{G}gt........ : 49295 88 : get Intron 2 >>>> {la}PheValArgSerAspLysProLysLeuPheArgGlyL : 100 484 bp {||}||||||||||||||||||||||||:!!!:!!:!|||| ++{la}PheValArgSerAspLysProLysMetTyrLysGlyL 49296 : .................ag{CT}TTTGTCAGGAGTGACAAGCCGAAGATGTACAAGGGTC : 49815 101 : euGlnIleLys >>>> Target Intron 3 >>>> TyrValArgGlySerAspPr : 110 ||||||||||| 182 bp |||||||||||||||||||| euGlnIleLys++ ++TyrValArgGlySerAspPr 49816 : TTCAGATCAAGgt.........................agTACGTGAGAGGCTCAGATCC : 50027 111 : oValLeuLysLeuLeuAspAspAsnGlyAsnIleAlaGluGluLeuSerIleLeuLysTr : 130 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| oValLeuLysLeuLeuAspAspAsnGlyAsnIleAlaGluGluLeuSerIleLeuLysTr 50028 : TGTACTAAAGCTTTTGGATGATAACGGGAACATCGCTGAAGAGCTCAGCATCCTCAAGTG : 50087 131 : pAsnThrAspSerValGluGluPheLeuSerGluLysLeuGluArgIle : 146 ||||||||||||||||||||||||||||||||||||||||!!:|||||| pAsnThrAspSerValGluGluPheLeuSerGluLysLeuAspArgIle 50088 : GAACACAGACAGTGTGGAGGAGTTCCTGAGTGAGAAATTAGACCGAATA : 50137 vulgar: Sel15 0 146 . KV884860.1:subseq(2942292,100000) 46844 50137 + 603 M 65 195 5 0 2 I 0 2185 3 0 2 M 21 63 S 0 1 5 0 2 I 0 480 3 0 2 S 1 2 M 16 48 5 0 2 I 0 178 3 0 2 M 43 129 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-22 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local gene 46845 50137 603 + . gene_id 1 ; sequence Sel15 ; gene_orientation + KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local cds 46845 47039 . + . KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local exon 46845 47039 . + . insertions 0 ; deletions 0 KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local splice5 47040 47041 . + . intron_id 1 ; splice_site "GT" KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local intron 47040 49228 . + . intron_id 1 KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local splice3 49227 49228 . + . intron_id 0 ; splice_site "AG" KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local cds 49229 49292 . + . KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local exon 49229 49292 . + . insertions 0 ; deletions 0 KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local splice5 49293 49294 . + . intron_id 2 ; splice_site "GT" KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local intron 49293 49776 . + . intron_id 2 KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local splice3 49775 49776 . + . intron_id 1 ; splice_site "AG" KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local cds 49777 49826 . + . KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local exon 49777 49826 . + . insertions 0 ; deletions 0 KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local splice5 49827 49828 . + . intron_id 3 ; splice_site "GT" KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local intron 49827 50008 . + . intron_id 3 KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local splice3 50007 50008 . + . intron_id 2 ; splice_site "AG" KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local cds 50009 50137 . + . KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local exon 50009 50137 . + . insertions 0 ; deletions 0 KV884860.1:subseq(2942292,100000) exonerate:protein2genome:local similarity 46845 50137 603 + . alignment_id 1 ; Query Sel15 ; Align 46845 1 195 ; Align 49229 66 63 ; Align 49779 88 48 ; Align 50009 104 129 # --- END OF GFF DUMP --- # -- completed exonerate analysis