| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,254,785 | A→G | 15.9% | D751D (GAT→GAC) | rpoC ← | RNA polymerase, beta prime subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,254,785 | 0 | A | G | 15.9% | 150.0 / 19.9 | 63 | D751D (GAT→GAC) | rpoC | RNA polymerase, beta prime subunit |
| Reads supporting (aligned to +/- strand): ref base A (23/30); new base G (6/4); total (29/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.92e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 | |||||||||||
AGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAC > minE/2254715‑2254852 | agaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCa < 1:3638102/71‑1 (MQ=255) gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt < 1:3370327/71‑1 (MQ=255) gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt < 1:1619030/71‑1 (MQ=255) aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg > 1:3747987/1‑71 (MQ=255) aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc < 1:3297121/71‑1 (MQ=255) aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc < 1:2367583/71‑1 (MQ=255) ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCg > 1:3662170/1‑69 (MQ=255) ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc > 1:2802122/1‑70 (MQ=255) ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc > 1:4010258/1‑70 (MQ=255) ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc > 1:3499390/1‑70 (MQ=255) tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg < 1:4690708/60‑1 (MQ=255) tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg < 1:3306006/60‑1 (MQ=255) tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg < 1:3641668/60‑1 (MQ=255) aCGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg < 1:1160143/59‑1 (MQ=255) cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt > 1:1581500/1‑60 (MQ=255) cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt > 1:1737997/1‑60 (MQ=255) cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt > 1:4047195/1‑60 (MQ=255) cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt > 1:294702/1‑60 (MQ=255) cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt > 1:2760828/1‑60 (MQ=255) agTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt > 1:461485/2‑60 (MQ=255) ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAcc < 1:1159017/71‑1 (MQ=255) ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc > 1:1283729/1‑70 (MQ=255) ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc > 1:4213329/1‑70 (MQ=255) cAGAGCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc > 1:1040516/1‑71 (MQ=255) cAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc > 1:3658936/1‑71 (MQ=255) aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc < 1:4016088/70‑1 (MQ=255) aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc < 1:3303827/70‑1 (MQ=255) aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACg < 1:3205649/71‑1 (MQ=255) aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATACGGCTTCGCCATCAGACCAc < 1:3090693/70‑1 (MQ=255) cTTCACGGAAGTTCGCGGTGATTGGCGTTTCGAAGATGGAGCCATCCgg < 1:2967508/49‑1 (MQ=255) ttCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg < 1:4973791/48‑1 (MQ=255) tCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa < 1:886499/67‑1 (MQ=255) cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa < 1:4818177/66‑1 (MQ=255) cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc > 1:3823471/1‑69 (MQ=255) cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc > 1:1726566/1‑69 (MQ=255) cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc < 1:4048458/71‑1 (MQ=255) cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc < 1:1202514/71‑1 (MQ=255) cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc < 1:7006/71‑1 (MQ=255) cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc < 1:4027435/71‑1 (MQ=255) cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc < 1:2769271/71‑1 (MQ=255) ggAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc > 1:1557053/1‑70 (MQ=255) aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg < 1:2906056/40‑1 (MQ=255) aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGc > 1:3509438/1‑71 (MQ=255) cgGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGa < 1:1402842/71‑1 (MQ=255) gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt < 1:5072192/71‑1 (MQ=255) gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt < 1:4459727/71‑1 (MQ=255) gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt > 1:4995036/1‑71 (MQ=255) gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt < 1:4762503/71‑1 (MQ=255) gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt > 1:3248259/1‑71 (MQ=255) tGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATc > 1:1870126/1‑71 (MQ=255) gCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGc < 1:5027828/70‑1 (MQ=255) tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt < 1:1519824/54‑1 (MQ=255) tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt < 1:4115351/54‑1 (MQ=255) tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt < 1:2461947/54‑1 (MQ=255) tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt > 1:2925638/1‑61 (MQ=255) tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt > 1:2466236/1‑61 (MQ=255) tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt > 1:5184857/1‑61 (MQ=255) tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt > 1:1859121/1‑61 (MQ=255) tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt > 1:4603414/1‑61 (MQ=255) tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt > 1:433424/1‑61 (MQ=255) tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt > 1:1575597/1‑61 (MQ=255) gCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAc < 1:621537/71‑1 (MQ=255) gCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAc < 1:3800115/71‑1 (MQ=255) ccATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAc < 1:1569903/70‑1 (MQ=255) | AGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAC > minE/2254715‑2254852 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |