Command line: [exonerate -m p2g --showtargetgff -q SelU1n/SelU1n_U.fa -t SelU1n/SUBSEQ/KB229058.1_0_subseq.txt --exhaustive yes] Hostname: [luke] C4 Alignment: ------------ Query: SPP00000112_2.0 # Protein # Selenoprotein X1 (SelX1) # Human (Homo sapiens) Target: KB229058.1:subseq(3421380,25179) Odobenus rosmarus divergens isolate Ivan unplaced genomic scaffold Scaffold181, whole genome shotgun sequence Model: protein2genome:local Raw score: 916 Query range: 0 -> 192 Target range: 18953 -> 25179 1 : MetSerPheLeuGlnAspProSerPhePheThrMetGlyMetTrpSerIleGlyAlaGly : 20 ||||||||||||||||||! !|||||||||!:!|||||||||||||||:!!||||||||| MetSerPheLeuGlnAspLeuSerPhePheSerMetGlyMetTrpSerLeuGlyAlaGly 18954 : ATGTCTTTCCTCCAGGATCTGAGTTTCTTCAGCATGGGGATGTGGTCCCTTGGCGCGGGG : 19011 21 : AlaLeuGlyAlaAlaAlaLeuAlaLeuLeuLeuAlaAsnThrAspValPheLeuSerLys : 40 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AlaLeuGlyAlaAlaAlaLeuAlaLeuLeuLeuAlaAsnThrAspValPheLeuSerLys 19012 : GCCCTGGGGGCGGCTGCCTTGGCCCTGCTCCTGGCCAACACAGACGTGTTTTTGTCCAAG : 19071 41 : ProGlnLysAlaAlaLeuGluTyrLeuGluAspIleAspLeuLysThrLeuGluLys{G} : 60 !!|||||||||.!!||||||||||||||||||||||||||||||||||||! !|||{|} SerGlnLysAlaThrLeuGluTyrLeuGluAspIleAspLeuLysThrLeuGlyLys{G} 19072 : TCCCAGAAGGCGACGCTGGAATATTTGGAGGACATAGACCTGAAAACCCTGGGAAAG{G} : 19131 61 : >>>> Target Intron 1 >>>> {lu}ProArgThrPheLysAlaLysGluLeu : 69 1993 bp {||}||||||||||||||||||||||||||| ++ ++{lu}ProArgThrPheLysAlaLysGluLeu 19132 : gt.........................ag{AA}CCAAGGACTTTCAAAGCAAAGGAGCTC : 21151 70 : TrpGluLysAsnGlyAlaValIleMetAlaValArgArgProGlyCysPheLeuCysArg : 89 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| TrpGluLysAsnGlyAlaValIleMetAlaValArgArgProGlyCysPheLeuCysArg 21152 : TGGGAAAAGAATGGAGCTGTGATTATGGCTGTGCGCAGGCCAGGCTGTTTCCTCTGCCGA : 21211 90 : Glu >>>> Target Intron 2 >>>> GluAlaAlaAspLeuSerSerLeuLysS : 100 ||| 2982 bp ||||||||||||||||||||||||||| Glu++ ++GluAlaAlaAspLeuSerSerLeuLysP 21212 : GAGgt.........................agGAGGCTGCGGACCTGTCTTCCCTGAAGC : 24226 101 : erMetLeuAspGlnLeuGlyValProLeuTyrAlaValValLysGluHisIleArgThrG : 120 !!! !||||||:!!|||||||||||||||||||||||||||||||||!!.|||||||||| roLysLeuAspGluLeuGlyValProLeuTyrAlaValValLysGluGlnIleArgThrG 24227 : CCAAGCTGGACGAGCTGGGAGTCCCCCTGTATGCCGTGGTGAAGGAACAAATCAGGACCG : 24286 121 : luValLysAspPheGlnProTyrPheLysGlyGluIlePheLeuAspGluLys >>>> : 138 |||||:!!||||||||||||||||||||||||||||||||||||||||||||| luValGlnAspPheGlnProTyrPheLysGlyGluIlePheLeuAspGluLys++ 24287 : AAGTGCAGGACTTCCAGCCTTATTTCAAAGGAGAAATCTTCCTGGATGAAAAGgt..... : 24342 139 : Target Intron 3 >>>> LysLysPheTyrGlyProGlnArgArgLysMetMetPh : 150 675 bp |||||||||||||||||||||||||||||||||||||| ++LysLysPheTyrGlyProGlnArgArgLysMetMetPh 24343 : ....................agAAAAAGTTCTACGGTCCCCAAAGACGGAAGATGATGTT : 25051 151 : eMetGlyPheIleArgLeuGlyValTrpTyrAsnPhePheArgAlaTrpAsnGlyGlyPh : 170 |||||||||||||||||||||||||||||||||||||||||||||| !!||||||||||| eMetGlyPheIleArgLeuGlyValTrpTyrAsnPhePheArgAlaArgAsnGlyGlyPh 25052 : TATGGGATTTATCCGTCTGGGTGTCTGGTACAACTTCTTCCGAGCCCGGAACGGAGGCTT : 25111 171 : eSerGlyAsnLeuGluGlyGluGlyPheIleLeuGlyGlyValPheValValGlySerGl : 190 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| eSerGlyAsnLeuGluGlyGluGlyPheIleLeuGlyGlyValPheValValGlySerGl 25112 : TTCTGGAAACCTGGAAGGCGAAGGCTTCATCCTTGGGGGAGTTTTTGTGGTGGGATCAGG : 25171 191 : yLysGln : 192 ||||||| yLysGln 25172 : AAAGCAG : 25179 vulgar: SPP00000112_2.0 0 192 . KB229058.1:subseq(3421380,25179) 18953 25179 + 916 M 59 177 S 0 1 5 0 2 I 0 1989 3 0 2 S 1 2 M 30 90 5 0 2 I 0 2978 3 0 2 M 47 141 5 0 2 I 0 671 3 0 2 M 55 165 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2015-11-25 ##type DNA # # # seqname source feature start end score strand frame attributes # KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local gene 18954 25179 916 + . gene_id 1 ; sequence SPP00000112_2.0 ; gene_orientation + KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local cds 18954 19131 . + . KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local exon 18954 19131 . + . insertions 0 ; deletions 0 KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local splice5 19132 19133 . + . intron_id 1 ; splice_site "GT" KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local intron 19132 21124 . + . intron_id 1 KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local splice3 21123 21124 . + . intron_id 0 ; splice_site "AG" KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local cds 21125 21216 . + . KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local exon 21125 21216 . + . insertions 0 ; deletions 0 KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local splice5 21217 21218 . + . intron_id 2 ; splice_site "GT" KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local intron 21217 24198 . + . intron_id 2 KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local splice3 24197 24198 . + . intron_id 1 ; splice_site "AG" KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local cds 24199 24339 . + . KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local exon 24199 24339 . + . insertions 0 ; deletions 0 KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local splice5 24340 24341 . + . intron_id 3 ; splice_site "GT" KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local intron 24340 25014 . + . intron_id 3 KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local splice3 25013 25014 . + . intron_id 2 ; splice_site "AG" KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local cds 25015 25179 . + . KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local exon 25015 25179 . + . insertions 0 ; deletions 0 KB229058.1:subseq(3421380,25179) exonerate:protein2genome:local similarity 18954 25179 916 + . alignment_id 1 ; Query SPP00000112_2.0 ; Align 18954 1 177 ; Align 21127 61 90 ; Align 24199 91 141 ; Align 25015 138 165 # --- END OF GFF DUMP --- # -- completed exonerate analysis