Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 70,638 | A→G | 6.9% | H34R (CAT→CGT) | spoIIE → | SpoIIAA‑phosphate serine phosphatase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 70,638 | 0 | A | G | 6.9% | 171.3 / 4.3 | 58 | H34R (CAT→CGT) | spoIIE | SpoIIAA‑phosphate serine phosphatase |
Reads supporting (aligned to +/- strand): ref base A (27/27); new base G (2/2); total (29/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 = 8.31e-01 |
GGCAACGAATATAACAGGTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCAGAAGAGACAA > BS168_reference_genome/70505‑70776 | ggCAACGAATATAACAGGTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTg > 4:312926/1‑138 (MQ=255) cGAATATAACAGGTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCa < 6:197157/139‑1 (MQ=255) gTGGGATATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCta > 1:231424/1‑138 (MQ=255) gTGGGAGATGAGAGGAATTGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCAAGAAGATGGTGACTCATCATATGCATTAACTTTTTATCta > 6:406481/1‑138 (MQ=255) gTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCta > 1:357043/1‑138 (MQ=255) gTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCta > 2:27114/1‑138 (MQ=255) gTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCta > 2:245533/1‑138 (MQ=255) gTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCta > 5:36427/1‑138 (MQ=255) gTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCta > 3:239750/1‑138 (MQ=255) gTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGGTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGGTTTTCta > 5:316503/1‑138 (MQ=255) tGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCtata < 4:148377/139‑1 (MQ=255) gaTGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGa < 1:78350/139‑1 (MQ=255) aggagaGTAATGGAAAAAGCAGAAAGAAGAGGGAAGCGGCCAATGGAGGGACAAGCTTTGGTAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATGATTTGCATTCACTGTTTTTCTATAAAGGa < 6:36427/136‑1 (MQ=255) gagGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCAt < 5:163281/139‑1 (MQ=255) aaTGGAAAAAGCAGAAAGATCAGTGAACGCGCCATAGGCGGGACAAGCTTTGGAAAAACTCCAATCGATTTTTAACAGGGGCACGAAGCTGGTCACTCATCATTTGCATTCACTGTTTTTCTATAAAGCATTCAtat > 5:375987/1‑135 (MQ=255) aaTGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTAt > 6:195001/1‑139 (MQ=255) aGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATc > 7:208577/1‑139 (MQ=255) aaagaaGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGAtt > 5:321436/1‑139 (MQ=255) gAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGa > 8:392034/1‑139 (MQ=255) aaCGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGa < 2:231424/138‑1 (MQ=255) aaCGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGAc < 3:392617/139‑1 (MQ=255) aaCGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGAc < 5:195001/139‑1 (MQ=255) aCGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACg < 7:392034/139‑1 (MQ=255) aCGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACg < 1:27114/139‑1 (MQ=255) gggCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACggg > 5:163415/1‑139 (MQ=255) gCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCt < 6:385805/139‑1 (MQ=255) gCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCt > 2:187356/1‑139 (MQ=255) gCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCt < 1:245533/139‑1 (MQ=255) tGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGTTCACGAAGCTGGTGAATCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCat < 4:239750/139‑1 (MQ=255) tGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCat < 3:312926/139‑1 (MQ=255) tGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCa > 8:225605/1‑138 (MQ=255) tGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCa > 4:265078/1‑138 (MQ=255) gCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCatat < 6:321436/139‑1 (MQ=255) ggACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCat > 5:413773/1‑134 (MQ=255) ggACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCat < 6:413773/134‑1 (MQ=255) ggACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCAtattat > 3:351630/1‑139 (MQ=255) ggACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCAtattat > 1:120532/1‑139 (MQ=255) ggACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCAtatta > 2:217732/1‑138 (MQ=255) aCAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATcc > 1:99321/1‑139 (MQ=255) aGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAg < 1:187356/139‑1 (MQ=255) ggAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCt < 6:83280/138‑1 (MQ=255) ggAAAAACTCCAATCGGTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCt < 2:371622/138‑1 (MQ=255) gAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCg > 1:192682/1‑130 (MQ=255) gAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCg < 2:192682/130‑1 (MQ=255) aCTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCGTTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTc > 7:368810/1‑139 (MQ=255) aCTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCGTTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTc > 7:274843/1‑139 (MQ=255) cggatattAACACGGCCACGAAGCTGGTGACTCCTCATTTGCATTCACTGTTTTGCTATAAAGGAGTCATTTAAGGTGTCATCGCATTTTTGCGTGGACGGGCTGTCATATTATCCGAGGTGCTTCCCTTCGCACTTcc < 4:360742/133‑1 (MQ=255) gtattttAACAGGCGCCCGACGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTAGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTATCGCACTTCct < 6:316503/136‑1 (MQ=255) gTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCct < 6:158895/139‑1 (MQ=255) ttttttAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCctt > 6:130270/1‑139 (MQ=255) ttttttAACAGCGGCACGAAGCTGGTGACTCATCATTTGAATTCACTGATTTTCTATAAAGGATTCATTAATGATGTCATCTGATTTTTTCTTGGACGGGCTTTCATATTATCCGAGTTGATTCCTTTCGACCTTcat > 8:18857/1‑136 (MQ=255) aGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTCTCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAg < 5:406481/139‑1 (MQ=255) ggCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGa < 2:357043/139‑1 (MQ=255) gAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGctt < 8:208577/139‑1 (MQ=255) ggTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTAt > 2:131620/1‑138 (MQ=255) aCTCATCGTTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCagaa < 8:274843/139‑1 (MQ=255) aCTCATCGTTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCagaa < 8:368810/139‑1 (MQ=255) aCTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCagaa < 6:163415/139‑1 (MQ=255) tcatcaTTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCagaa > 8:81757/1‑137 (MQ=255) atcatTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCg > 2:367128/1‑99 (MQ=255) atcatTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCg < 1:367128/99‑1 (MQ=255) atcatTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCAGAAgaga > 3:353689/1‑139 (MQ=255) tcatTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAg < 7:72879/116‑1 (MQ=255) tcatTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAg > 8:72879/1‑116 (MQ=255) atTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCAGAAGAGACaa < 4:353689/139‑1 (MQ=255) | GGCAACGAATATAACAGGTGGGAGATGAGAGGAATGGAAAAAGCAGAAAGAAGAGTGAACGGGCCAATGGCGGGACAAGCTTTGGAAAAACTCCAATCGTTTTTTAACAGGGGCACGAAGCTGGTGACTCATCATTTGCATTCACTGTTTTTCTATAAAGGATTCATTTATGTTGTCATCGGATTTTTGCTTGGACGGGCTTTCATATTATCCGAGGTGCTTCCTTTCGCACTTCCTTTCTTTGGAGCGATGCTTCTTATCAGAAGAGACAA > BS168_reference_genome/70505‑70776 |
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: ‑ |