| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 890,115 | 0 | G | C | 41.7% | 30.6 / 11.1 | 26 | G9G (GGG→GGC) | ybjM | inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base G (11/3); new base C (0/10); total (11/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTT > NC_000913/890055‑890182 | aTGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATct < 1:632530/68‑1 (MQ=255) ttCGATTGACGCAGTTTCAGGAGGATCGAGTTAAACATAAACTTCGTTGGGCGGTTTCAATCTGCTGt < 2:739524/68‑1 (MQ=255) tCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGtt > 2:1167328/1‑68 (MQ=255) cGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATCTGCTGttt < 2:242330/68‑1 (MQ=255) cGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCctc > 1:713738/1‑68 (MQ=255) gCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATCTGCTGTTTTGTCctct < 2:832510/68‑1 (MQ=255) gCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATCTGCTGTTTTGTCctct < 1:998769/68‑1 (MQ=255) agcgtCAGTAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCTAAATCTGCTGTTTTGTCCTCTTc < 2:961045/65‑1 (MQ=255) aGTGTCAGGAGGAACGAGTGAAACATAATCAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTc < 1:455318/68‑1 (MQ=255) aGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATCTGCTGTTTTGTCCTCTTc < 1:253137/68‑1 (MQ=255) aGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATCTGCTGTTTTGTCCTCTTc < 2:346872/68‑1 (MQ=255) aGTGGCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATCTGCTGTTTTGTCCTCTTc < 2:422818/68‑1 (MQ=255) gtgtCAGGAGGAACGAGTGAAACATAAACAACGTTGGGC‑GGGGCAATCTGCTGTTTTGTCCTCTTCAt < 2:701222/68‑1 (MQ=255) tgtCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCAt > 1:599884/1‑68 (MQ=255) gaggaACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTgtggt > 2:991781/1‑68 (MQ=255) gaggaACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTATGTCCTCTACATTgtggt > 2:1138104/1‑68 (MQ=255) gaACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGtgtg < 1:1167330/68‑1 (MQ=255) gaACGAGTGAAACATAAACAACGTTGGGCGGGCGCAATCTGCTGTTTTGTCCTCTTCATTGTGGtgtg < 1:349955/68‑1 (MQ=255) taCGAGTGAAACATAAACAACGTTGGGCGTTTTCAATCTTCTGTTTTGTCCTCTTCATTGTGgtgttc < 1:283424/67‑3 (MQ=255) cGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCt > 1:49311/1‑68 (MQ=255) acaacGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGc > 2:356862/1‑68 (MQ=255) acaacGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGc > 2:894844/1‑68 (MQ=255) acaacGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCATTTTCTGGCGACGc > 2:473416/1‑68 (MQ=255) caacGTTGGGCGGGGGCAAACTGCTGTTTAGTCCTCTTCATTGTGGTGTGCCTTATTCTGGCGACGca > 2:449013/1‑68 (MQ=255) caacGTTGGGCGGGCGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGca < 2:830835/68‑1 (MQ=255) gTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCATTTTTCTGGCGACGCACATg > 1:93186/1‑68 (MQ=255) ggCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCtt > 1:296473/1‑68 (MQ=255) | ATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTT > NC_000913/890055‑890182 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |