Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q ../SelH-human.fa -t SelH-NIVO01034600.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000014_1.0 # Protein # Selenoprotein H (SelH) # Homo sapiens # Complete Target: NIVO01034600.1:subseq(4978439,100353) Ammotragus lervia isolate NAZ scaffold37_len5366835_cov71, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 474 Query range: 0 -> 122 Target range: 50574 -> 49999 1 : MetAlaProArgGlyArgLysArgLysAlaGluAlaAlaValValAlaValAlaGluLys : 20 |||||| !!||||||||||||||||||||||||||||||:!!!.!|||!.!||||||||| MetAlaSerArgGlyArgLysArgLysAlaGluAlaAlaLeuAlaAlaAlaAlaGluLys 50574 : ATGGCGTCCCGCGGGAGGAAGCGGAAGGCCGAGGCGGCGCTGGCCGCTGCGGCCGAGAAG : 50517 21 : ArgGluLysLeuAlaAsnGlyGlyGluGly<->MetGluGluAlaThrValValIleGlu : 39 |||||||||! !|||..!||| !|||||| :!!! !! ! !!!:!|||||||||||| ArgGluLysProAlaGlyGlyGlnGluGlyGluValAlaGlyProSerValValIleGlu 50516 : CGGGAGAAGCCGGCCGGCGGCCAGGAGGGGGAAGTGGCGGGGCCGAGCGTCGTCATCGAG : 50457 40 : His{Cy} >>>> Target Intron 1 >>>> {s}ThrSerUnkArgValTyrGly : 48 |||{||} 83 bp {|}|||||| |||||||||||| His{Cy}++ ++{s}ThrSer***ArgValTyrGly 50456 : CAC{TG}gt.........................ag{C}ACGAGCTGACGAGTCTACGGG : 50347 49 : ArgAsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGluAlaProGluLeuProValLys : 68 ||||||||||||||||||||||||||||||||||||:!!|||||||||||| !!|||||| ArgAsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGlnAlaProGluLeuAlaValLys 50346 : CGCAACGCCGCGGCCCTGAGCCAGGCGCTGCGCCTGCAGGCCCCCGAGCTGGCGGTGAAG : 50287 69 : ValAsnProThrLysProArgArgGlySerPheGluValThrLeuLeuArgProAspGly : 88 |||||||||.!!!:!|||||||||||||||||||||||||||||||||||| !!|||||| ValAsnProAlaArgProArgArgGlySerPheGluValThrLeuLeuArgAlaAspGly 50286 : GTGAACCCCGCCAGGCCGCGGAGGGGCAGCTTCGAGGTGACGCTGCTGCGCGCCGACGGC : 50227 89 : Ser{S} >>>> Target Intron 2 >>>> {er}AlaGluLeuTrpThrGlyIle : 97 |||{|} 123 bp {||}||||||||||||||||||:!! Ser{S}++ ++{er}AlaGluLeuTrpThrGlyLeu 50226 : AGC{A}gt.........................ag{GC}GCGGAGCTCTGGACGGGTCTT : 50077 98 : LysLysGlyProProArgLysLeuLysPheProGluProGlnGluValValGluGluLeu : 117 |||||||||||||||||||||||||||||||||||||||!!.! !|||:!!||||||||| LysLysGlyProProArgLysLeuLysPheProGluProHisValValLeuGluGluLeu 50076 : AAGAAGGGGCCCCCACGCAAACTCAAGTTTCCGGAGCCTCACGTGGTGCTGGAGGAGCTG : 50017 118 : LysLysTyrLeuSer : 122 ||||||||||||||| LysLysTyrLeuSer 50016 : AAGAAGTACCTTTCG : 50000 vulgar: SPP00000014_1.0 0 122 . NIVO01034600.1:subseq(4978439,100353) 50574 49999 - 474 M 30 90 G 0 3 M 10 30 S 0 2 5 0 2 I 0 79 3 0 2 S 1 1 M 48 144 S 0 1 5 0 2 I 0 119 3 0 2 S 1 2 M 32 96 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local gene 50000 50574 474 - . gene_id 1 ; sequence SPP00000014_1.0 ; gene_orientation + NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local cds 50450 50574 . - . NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local exon 50450 50574 . - . insertions 3 ; deletions 0 NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local splice5 50448 50449 . - . intron_id 1 ; splice_site "GT" NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local intron 50367 50449 . - . intron_id 1 NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local splice3 50367 50368 . - . intron_id 0 ; splice_site "AG" NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local cds 50221 50366 . - . NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local exon 50221 50366 . - . insertions 0 ; deletions 0 NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local splice5 50219 50220 . - . intron_id 2 ; splice_site "GT" NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local intron 50098 50220 . - . intron_id 2 NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local splice3 50098 50099 . - . intron_id 1 ; splice_site "AG" NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local cds 50000 50097 . - . NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local exon 50000 50097 . - . insertions 0 ; deletions 0 NIVO01034600.1:subseq(4978439,100353) exonerate:protein2genome:local similarity 50000 50574 474 - . alignment_id 1 ; Query SPP00000014_1.0 ; Align 50575 1 90 ; Align 50482 31 30 ; Align 50366 42 144 ; Align 50096 91 96 # --- END OF GFF DUMP --- # -- completed exonerate analysis