| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 157,237 | T→C | 52.9% | P17P (CCA→CCG) | btuF ← | vitamin B12 transporter subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 157,237 | 0 | T | C | 52.9% | 6.5 / 58.3 | 70 | P17P (CCA→CCG) | btuF | vitamin B12 transporter subunit |
| Reads supporting (aligned to +/- strand): ref base T (21/12); new base C (13/24); total (34/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.04e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01 | |||||||||||
GGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTT > minE/157177‑157296 | ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa > 1:2025753/1‑68 (MQ=255)ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa > 1:867947/1‑68 (MQ=255)ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa > 1:1115255/1‑68 (MQ=255)ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa > 1:1681464/1‑68 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:1618794/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:2242587/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:195580/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:280808/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:2895811/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:1701905/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:2902772/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:452317/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:589132/69‑1 (MQ=255) ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:1307685/69‑1 (MQ=255) ttCAGTGTTCGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:1889126/69‑1 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGaaaa > 1:1486485/1‑65 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa > 1:3255873/1‑69 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa > 1:2888518/1‑69 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa > 1:2408796/1‑69 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa > 1:1952064/1‑69 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa > 1:865300/1‑69 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa > 1:1103349/1‑69 (MQ=255) cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa > 1:949111/1‑69 (MQ=255) gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:3056123/69‑1 (MQ=255) gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:2331715/69‑1 (MQ=255) gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:329390/69‑1 (MQ=255) gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:2272093/69‑1 (MQ=255) gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:120550/69‑1 (MQ=255) gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:1349291/69‑1 (MQ=255) gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:1789629/69‑1 (MQ=255) tgatgaCGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc < 1:2322932/68‑1 (MQ=255) gatgaCGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa < 1:811066/45‑1 (MQ=38) cggcgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGc < 1:3079822/69‑1 (MQ=255) ggcgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGc < 1:2878466/68‑1 (MQ=255) gcgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCa < 1:1497878/69‑1 (MQ=255) cgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCa < 1:1081336/68‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:21503/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:1113678/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:455007/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:1425628/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:3233735/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:1742025/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:2883458/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:2042231/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:2190159/63‑1 (MQ=255) cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:2358906/63‑1 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:2476514/63‑1 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:3337974/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:2206324/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:80751/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:755123/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:2041870/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:469923/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:460738/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:1038039/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:3250186/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:3244723/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:2387336/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:1808901/2‑67 (MQ=255) attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca > 1:2443549/2‑67 (MQ=255) tGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg < 1:2852459/62‑1 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:1463372/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:1492078/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:1931471/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:528192/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:166052/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:2211823/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:84383/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:1702071/1‑69 (MQ=255) gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt > 1:2723695/1‑69 (MQ=255) | GGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTT > minE/157177‑157296 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |