Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 3,701,855 | C→A | 6.8% | P149T (CCT→ACT) | rbsR → | ribose operon repressor (LacI family, D‑ribose) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 3,701,855 | 0 | C | A | 6.8% | 218.2 / 4.3 | 68 | P149T (CCT→ACT) | rbsR | ribose operon repressor (LacI family, D‑ribose) |
Reads supporting (aligned to +/- strand): ref base C (31/32); new base A (3/2); total (34/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01 |
CCTTCAAACCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGACCGATTTAACGGCGCTT > BS168_reference_genome/3701731‑3701990 | ccTTCAAACCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGt > 7:155435/1‑139 (MQ=255) cTTCAAACCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGt > 2:61818/1‑138 (MQ=255) cAAACCTTTAAGCAAAATCTTGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGa < 1:127722/137‑1 (MQ=255) cAAACCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGAc < 8:250035/138‑1 (MQ=255) cAAACCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGAc < 8:155435/138‑1 (MQ=255) cAAACCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGAc < 3:106572/138‑1 (MQ=255) cAAACCTTTAAGCAAAATCATGTCGCAGGAATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGATTTATCCAGGTGTTTTTCTGGACAGAACGCGTGAAGGGGCTCCTTCTGTGTCCAGTGAc < 4:196405/138‑1 (MQ=255) catccTTTAAGCTAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATTCAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGATTGAAGGGGCTCCTTCTGTGTCCAGTGACgg < 8:163382/136‑1 (MQ=255) aaaCCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACgg < 3:52993/139‑1 (MQ=255) ttAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGAc > 4:346089/1‑131 (MQ=255) ttAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGAc < 3:346089/131‑1 (MQ=255) ttAAGCAAAACCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTTGACAGAACGCATGCAGGGGCTCCTTATGTGTCCAGTGACGGCtata > 6:41106/1‑138 (MQ=255) tAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCtata > 2:345450/1‑137 (MQ=255) tAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCtata < 7:85083/137‑1 (MQ=255) tAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCtata < 7:279527/137‑1 (MQ=255) aaTCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTa > 6:104404/1‑138 (MQ=255) aaTCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTa > 3:184337/1‑138 (MQ=255) aGGTATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc < 1:60949/139‑1 (MQ=255) aGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc < 8:251130/139‑1 (MQ=255) gCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCa < 1:345450/139‑1 (MQ=255) tattGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc < 7:322612/133‑1 (MQ=255) tattGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc < 7:31687/133‑1 (MQ=255) tattGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc > 8:31687/1‑133 (MQ=255) tattGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc > 8:322612/1‑133 (MQ=255) cAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCt > 3:48839/1‑86 (MQ=255) cAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCt < 4:48839/86‑1 (MQ=255) cAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTGTTCTGGACAGAACGCTTGAAGGGGCTCCATCAGGGTCCACTGCCGGCTATACACGAGAAAAGTTAGCCGCCCAGGCTAACATACa > 5:123015/1‑139 (MQ=255) aCGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCAGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCAATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACg > 6:200889/1‑139 (MQ=255) aTCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAgcg > 6:106891/1‑139 (MQ=255) gAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAgcg < 6:189547/134‑1 (MQ=255) ggAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACg < 1:61818/139‑1 (MQ=255) cGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTga > 4:153932/1‑137 (MQ=255) cGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTga > 2:46377/1‑137 (MQ=255) cGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTga > 1:171102/1‑137 (MQ=255) gCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAgg > 7:69087/1‑139 (MQ=255) cATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGa > 7:191972/1‑63 (MQ=255) cATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGa < 8:191972/63‑1 (MQ=255) cATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc > 4:139908/1‑91 (MQ=255) cATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAgccgcc < 3:139908/91‑1 (MQ=255) tGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAgg > 5:321060/1‑136 (MQ=255) tGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAgg < 6:321060/136‑1 (MQ=255) aaTTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCg < 5:106891/139‑1 (MQ=255) aTTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCGGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCg < 7:342130/138‑1 (MQ=255) ttATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCt < 8:69087/139‑1 (MQ=255) tATCGAGTTGTTTTTCTGGACAGACAGCTGTTAGGGGCTCCTTCTGTGTCCTGTGACGGCTATACAGGAGTAATGTTAGCCGCCGAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGCGAGGACCCGCt < 6:123015/138‑1 (MQ=255) tCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTACTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCa > 6:243233/1‑96 (MQ=255) tCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTACTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCa > 6:238767/1‑96 (MQ=255) tCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTACTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCa < 5:243233/96‑1 (MQ=255) tCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTACTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCa < 5:238767/96‑1 (MQ=255) ccAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACg > 3:89514/1‑117 (MQ=255) ccAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACg < 4:89514/117‑1 (MQ=255) ttgttTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAg > 3:214864/1‑72 (MQ=255) ttgttTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAg < 4:214864/72‑1 (MQ=255) ttgttTTTCTGGACAGAACGCTTGAAGGGGCTACTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCTGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACATGcc > 7:87581/1‑139 (MQ=255) ttttgCTGGACAGAACTCTTGACCGGGCTCCTTCGGTGGCCAGTGACGGCTATACAGGCGTAACGTTAGCCGCCCCGGCTATCATTCACGGAAAAAGCCCGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGAc < 1:161711/139‑1 (MQ=255) cTGGACAGAACGCTTGCAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCt < 8:369194/138‑1 (MQ=255) cAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGAc < 1:46377/139‑1 (MQ=255) gAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCGTGAGAGGACCCGCTCACCTGCCGACCGCTCAGGACCg > 3:39508/1‑139 (MQ=255) aaCGCTTGAAGGGGGTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGACCGa > 5:310061/1‑139 (MQ=255) aaCGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACg < 4:14071/97‑1 (MQ=255) aaCGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACg > 3:14071/1‑97 (MQ=255) aaCGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGACCGa > 6:144143/1‑139 (MQ=255) aaCGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGACCGa > 6:277070/1‑139 (MQ=255) ttGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGACCGATTTaa > 8:182713/1‑139 (MQ=255) gCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGACCGATTTAACGGCGCtt > 4:192109/1‑139 (MQ=255) gCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCAGCTCACCTGCCGACCGCTCACGACCGATTTAACGGCGCtt > 6:245542/1‑139 (MQ=255) cTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACg > 7:339961/1‑83 (MQ=255) cTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACg < 8:339961/83‑1 (MQ=255) | CCTTCAAACCTTTAAGCAAAATCATGTCGCAGGCATTATTGCCGCAACGAATTACCCGGATCTCGAGGAATACAGCGGCATGAATTATCCAGTTGTTTTTCTGGACAGAACGCTTGAAGGGGCTCCTTCTGTGTCCAGTGACGGCTATACAGGAGTAAAGTTAGCCGCCCAGGCTATCATTCACGGAAAAAGCCAGCGCATTACGCTCTTGAGAGGACCCGCTCACCTGCCGACCGCTCAAGACCGATTTAACGGCGCTT > BS168_reference_genome/3701731‑3701990 |
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: ‑ |