| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,883,667 | T→C | 35.8% | Q272Q (CAA→CAG) | rumA ← | 23S rRNA (uracil‑5)‑methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,883,667 | 0 | T | C | 35.8% | 48.6 / 46.1 | 78 | Q272Q (CAA→CAG) | rumA | 23S rRNA (uracil‑5)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base T (24/26); new base C (16/12); total (40/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTAAAAGTTAAGCGCAACCCGTTTGAGT > minE/1883604‑1883734 | cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:2226747/69‑1 (MQ=255)cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:268530/69‑1 (MQ=255)cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:1827694/69‑1 (MQ=255)cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:3208038/69‑1 (MQ=255)cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:2602291/69‑1 (MQ=255)cAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:328782/69‑1 (MQ=255) aGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:2826666/68‑1 (MQ=255) aGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:1804816/68‑1 (MQ=255) aGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:2695644/68‑1 (MQ=255) tCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTcaca > 1:2602469/1‑69 (MQ=255) tCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTcaca > 1:1624337/1‑69 (MQ=255) gTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCg < 1:2121798/69‑1 (MQ=255) gTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCg < 1:2042004/69‑1 (MQ=255) cgcgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:300710/58‑1 (MQ=255) cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCTCTACCATTTTTTGGTTCACACCCGCATTg < 1:2519852/69‑1 (MQ=255) cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCTCTACCATTTTTTGGTTCACACCCGCATTg < 1:2515780/69‑1 (MQ=255) cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg < 1:2263245/69‑1 (MQ=255) cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg < 1:3124318/69‑1 (MQ=255) cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg < 1:2072153/69‑1 (MQ=255) cgATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTg < 1:1506049/69‑1 (MQ=255) tCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:1574612/53‑1 (MQ=255) cTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGAcc > 1:2771824/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg > 1:1031326/1‑67 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg > 1:2880100/1‑67 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg > 1:2758002/1‑67 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg > 1:3152075/1‑67 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTg > 1:662308/1‑67 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:3079580/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:3008571/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:1720517/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:1640269/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:3212404/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:3220960/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:1424200/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:3434285/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:697053/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:908470/1‑68 (MQ=255) cAGGTTCCACATCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGa > 1:2742140/1‑68 (MQ=255) gCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGtt < 1:3378823/43‑1 (MQ=255) aTCCAGCCATTCCAACGCGCGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:1105614/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:687342/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:3154999/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:459159/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:2977331/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:2490703/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:1576171/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:2495664/69‑1 (MQ=255) aTCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:2781495/69‑1 (MQ=255) tCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATcgcgc > 1:2910489/1‑69 (MQ=255) tCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATcgcgc > 1:779165/1‑69 (MQ=255) tCCAGCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:1479676/68‑1 (MQ=255) cAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCgg > 1:277281/1‑69 (MQ=255) gCCATTCCAACGCACGCGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:284817/64‑1 (MQ=255) aTTCCAACGCACTCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaa < 1:1348338/69‑1 (MQ=255) aTTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaa < 1:735606/69‑1 (MQ=255) aTTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaa < 1:3377137/69‑1 (MQ=255) aTTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaa < 1:982148/69‑1 (MQ=255) aTTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaa < 1:3153783/69‑1 (MQ=255) aTTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaa < 1:2430343/69‑1 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGCCCTGAATAAAATCGCGCGGGCTaaaa > 1:1160502/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:503197/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:1334923/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:1813628/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:1250144/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:3160889/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:3090776/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:2346183/1‑68 (MQ=255) ccAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTaaaa > 1:964563/1‑68 (MQ=255) cgcacgcGCTACCATTTTCTGGTTCACACCCGCATTGACCTGAATAAAATcgcg < 1:574066/54‑1 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATaaa > 1:747163/1‑44 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATaaa > 1:2605598/1‑44 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATaaa > 1:554436/1‑44 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATaaa > 1:3333558/1‑44 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATaaa > 1:1501810/1‑44 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATaaa > 1:1903958/1‑44 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATaaa > 1:2624252/1‑44 (MQ=255) cgcgCTACCATTTTTTGGTTCACACCCGCATTGACATGAATaaa > 1:1626627/1‑44 (MQ=255) aCCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATc > 1:2286758/1‑41 (MQ=255) ttGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTAAAAGTTAAGCGCAACCCGTTTGAGt < 1:159504/68‑1 (MQ=255) ttGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTAAAAGTTAAGCGCAACCCGTTTGAGt < 1:2874448/68‑1 (MQ=255) | CAGATCCAGTACGCGATCTTCAGGTTGCACATCCAGCCATTCCAACGCACGCGCTACCATTTTTTGGTTCACACCCGCATTGACCTGAATAAAATCGCGCGGGCTAAAAGTTAAGCGCAACCCGTTTGAGT > minE/1883604‑1883734 |
| 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: ‑ |