Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 166,037:1 | +T | 100% | intergenic (+4/‑27) | trnI‑Gly → / → trnI‑Arg | tRNA‑Gly/tRNA‑Arg |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 166,037 | 1 | . | T | 96.2% | 134.5 / ‑0.5 | 52 | intergenic (+4/‑27) | trnI‑Gly/trnI‑Arg | tRNA‑Gly/tRNA‑Arg |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base T (22/28); total (23/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑A‑CCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTAATAGAATAGATAGGAAA > BS168_reference_genome/165900‑166166 | gCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑at < 5:66152/139‑2 (MQ=255) tAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTACGCTCCATTT‑atcc > 6:247927/1‑135 (MQ=255) aaaTGTATTATTACCAATGGTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑at < 2:274135/135‑2 (MQ=255) aaTGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑atccta < 8:259031/138‑6 (MQ=255) aaTGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTaa < 6:99489/139‑1 (MQ=255) tattatTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACggg < 5:303545/139‑1 (MQ=255) tattatTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACggg < 7:58364/139‑1 (MQ=255) ttattaCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACggg > 1:207372/1‑137 (MQ=255) ttattaCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACggg < 2:207372/137‑1 (MQ=255) ttattaCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAAACCCGTATTGCGCTCCATTT‑ATACTAACGGGa > 6:141057/1‑138 (MQ=21) ttattaCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGat > 3:90256/1‑139 (MQ=255) tattaCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTGTTCCGCTCCATTT‑ATCCTAACGGGata > 7:231171/1‑139 (MQ=255) taCCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGTATGCCCTGAGCCGCTCCATTTAA‑CATAACGCGATTTTg > 6:28603/1‑139 (MQ=21) tctcATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGACCACCACCTTGCCAAGGTGGGGGTCGCGGGGTCGAATCCCCTCTGCCGCCCCAATT‑ATCCTAAAGGGATATTGTAaattgaat > 3:311984/1‑139 (MQ=21) ctcATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTACGCTCCATTT‑ATGCTAACGGGATATTTTAaattgaatt > 7:216284/1‑139 (MQ=21) aTAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTgc < 7:58190/138‑1 (MQ=255) ggAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGGCTTCCGGTCCATTT‑ATACTGAACGGTTTTTGTTAATTTAATTGCGCCCg > 5:197939/1‑139 (MQ=17) gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 1:24681/139‑1 (MQ=255) gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 5:121326/139‑1 (MQ=255) gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 5:210800/139‑1 (MQ=255) gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 6:299030/139‑1 (MQ=255) gaAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAg < 3:49297/139‑1 (MQ=255) aaGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCt < 8:21669/139‑1 (MQ=255) aaGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCt < 6:311906/139‑1 (MQ=255) aaGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCt < 3:198145/139‑1 (MQ=255) cTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCaa < 8:82812/139‑1 (MQ=255) ttATTGCGATATGCGGAAGTAGTGCAGTGGTAGACCACCACCTTGCCACGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCTATTGGa < 4:174347/138‑1 (MQ=21) tGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATagag > 3:47193/1‑139 (MQ=255) tGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATagag > 1:71928/1‑139 (MQ=255) ggAAGTTGTTCAGTGGTATAACACCACATTGCCATGGGGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTa < 2:151192/139‑1 (MQ=25) ggAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTACCGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGTATAGAGCGTTTGACTa < 8:41042/139‑1 (MQ=255) gAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTAc < 1:251526/139‑1 (MQ=255) tagtagtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGATTCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTAc < 7:142187/137‑1 (MQ=255) gtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGAtg > 3:121556/1‑137 (MQ=255) gtTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCGTCCGCTGCATTT‑ATACTAACGGGATATTTTAAATTGAAATGCGCCCGTAGCTCAATTGGATAGAGCGTTTGAGTAAGGATc > 8:163173/1‑138 (MQ=21) gtTCAGTGGTACAACACCACCTTACCAAAGTAGGGGTACAGGCTTCGAATAACCAATTCCCCTCCATTT‑ATACGAACGCGATATTGAAAATTTCAATCCGCCCGTATCTCAATTGGAGTCAGCGTTTGACTACGGATc > 8:74060/1‑138 (MQ=255) ttCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCaa < 3:227724/139‑1 (MQ=255) ttCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCaa < 4:216193/139‑1 (MQ=255) tCAGTGGTAGAACACCACCTGGCCAAGGTGGGTGTCGCGGGTTCGTAGCGCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCaa < 5:97698/138‑1 (MQ=37) gTGGTAGAACACCACCTTGCCAAAGTGGGGGTCGCGGGTTCGAATCCCGTGGTCCGCGCCATTT‑ATCCTAACGGGATATTGTTAATTTAATTTCGCCCGTTGCTCAATTGGTTAGCGCGTTTGGCTACGGTTCAAAAgg > 7:125677/1‑139 (MQ=21) tAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAg < 6:239897/139‑1 (MQ=255) cacgTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTTACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGa < 7:163173/139‑1 (MQ=255) caccTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGa < 4:90256/139‑1 (MQ=255) cAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTGCATTT‑ATACTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGTTAGAGCGTTTGACGACGGATCAAAAGGTTAGGGGTTCGACTCCgctc > 6:303782/1‑139 (MQ=25) ggggTCGCGGGTTCGAAACCCGTCTTACGCTCCATTT‑ATACTTAAGGGATATTTTAAATTTAATTTCTCCCGTAGCTCAATTGGAAAGAGCGTTTGACAACGGATCAAAAAGATTGGGTTTGGTCTGCACACGGgcgcg > 6:212613/1‑139 (MQ=255) gTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTTCT‑CCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCa < 5:98063/139‑1 (MQ=255) gcgGGTTCGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCata < 5:303782/138‑1 (MQ=255) cgGGTTCGAATCCCGTGTTACCCGCCAATT‑ATGCATAAGGGATATTGTAAATTGAATTTCGCCCCTAGCACAATTTGAATGAGCGGTTTTCTACAGATCAAAAGGTTAGGGGGTCGACGCCGCGCGGGCGCGCCatg > 2:170882/1‑136 (MQ=255) cgGGTTCGAATCCCGTCTTACGCTCCATTT‑ATACTAACGGGATATTGTTAATTGAATTGCGCCAGTAGCTCAATTGGATATAGCGTTTGACTAGGGAGAAAAAGGTTAGGGGTTCGACTCCTATCGGGCGGGCCagata > 8:137267/1‑138 (MQ=18) cGAATCCCGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCATTTaa < 5:220590/139‑1 (MQ=255) cGTCTTCCGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTaatagaata > 7:32692/1‑139 (MQ=255) cGTCTTCCGCTCCATTT‑ATCCAAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTaatagaat > 7:37311/1‑138 (MQ=255) ccGCTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGAGTTTGACTACGGATCAATAGGTTAGGGTTTCGACTCCTCTCGGGCGCGCCATATCTTTTTATAGAatagatat > 4:29288/1‑137 (MQ=17) cTCCATTT‑ATCCTAACGGGATATTGTAAATTGATTTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCCTATCTTTTAATTGAATAGATAGGaaa > 2:137119/1‑139 (MQ=17) cTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGtt < 1:97056/85‑1 (MQ=255) cTCCATTT‑ATCCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGtt > 2:97056/1‑85 (MQ=255) | GCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAGTTCAGTGGTAGAACACCACCTTGCCAAGGTGGGGGTCGCGGGTTCGAATCCCGTCTTCCGCTCCATTT‑A‑CCTAACGGGATATTGTAAATTGAATTGCGCCCGTAGCTCAATTGGATAGAGCGTTTGACTACGGATCAAAAGGTTAGGGGTTCGACTCCTCTCGGGCGCGCCATATCTTTTAATAGAATAGATAGGAAA > BS168_reference_genome/165900‑166166 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |