| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,234,608 | G→C | 30.5% | intergenic (‑368/‑247) | rrsD ← / → purH | 16S ribosomal RNA/fused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,234,608 | 0 | G | C | 30.5% | 39.2 / 21.5 | 36 | intergenic (‑368/‑247) | rrsD/purH | 16S ribosomal RNA/fused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (10/15); new base C (3/8); total (13/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.08e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01 | |||||||||||
CTCCGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTGCCCA > minE/2234546‑2234675 | ctcCGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGc > 1:939844/1‑68 (MQ=255)ctcCGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGc > 1:835118/1‑68 (MQ=255)ctcCGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGc > 1:2026180/1‑68 (MQ=255) ccGCAGGCCGCAATAGAACAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCtt < 1:1581612/68‑1 (MQ=255) ccGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCtt < 1:418375/68‑1 (MQ=255) ccGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCtt < 1:144702/68‑1 (MQ=255) ccGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCtt < 1:2697788/68‑1 (MQ=255) ccGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCtt < 1:1765533/68‑1 (MQ=255) ccGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGAtt < 1:326604/68‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:965746/61‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:810866/61‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:732437/61‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:395421/61‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:3252332/61‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:2211342/61‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:1816313/61‑1 (MQ=255) cGCAATAGAAAAATTGCATAAAAATGACTGACTGCTGCATTCCCCAGCAAAACCCCGCttt < 1:1729028/61‑1 (MQ=255) gAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCtt < 1:1139918/53‑1 (MQ=255) ttGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGt < 1:1568160/69‑1 (MQ=255) tGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGtt < 1:2341606/69‑1 (MQ=255) aaaaTGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCcac > 1:110579/1‑69 (MQ=255) aaaaTGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCcac > 1:1598999/1‑69 (MQ=255) aaaTGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCcaca < 1:1773456/69‑1 (MQ=255) aaaTGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCcaca < 1:483525/69‑1 (MQ=255) ctgcATTCCCCAGCAAAACCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTaa > 1:1960500/1‑69 (MQ=255) ctgcATTCCCCAGCAAAACCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTaa > 1:1320738/1‑69 (MQ=255) ctgcATTCCCCAGCAAAACCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTaa > 1:1825267/1‑69 (MQ=255) ttCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTc > 1:1609614/1‑68 (MQ=255) ttCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTc > 1:2761643/1‑68 (MQ=255) ttCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTc > 1:1194484/1‑68 (MQ=255) cAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTg < 1:140644/69‑1 (MQ=255) cAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTg < 1:476623/69‑1 (MQ=255) cAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTg < 1:2056929/69‑1 (MQ=255) aaaGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTGcc < 1:2387593/69‑1 (MQ=255) aGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTGCCCa > 1:1297600/1‑69 (MQ=255) aGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTGCCCa > 1:2021830/1‑69 (MQ=255) | CTCCGCAGGCCGCAATAGAAAAATTGCAGAAAAATGACTGACTGCTGCATTCCCCAGCAAAAGCCCGCTTTATACCTTTTTACGCACAGAGTTATCCACAATCATCAATGTAATTTCTGTATTTTGCCCA > minE/2234546‑2234675 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |