| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,135,975 | 0 | G | C | 32.1% | 39.3 / 21.9 | 28 | A138P (GCC→CCC) | flgH | flagellar protein of basal‑body outer‑membrane L ring |
| Reads supporting (aligned to +/- strand): ref base G (17/2); new base C (0/9); total (17/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.96e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAAC > NC_000913/1135908‑1136037 | gCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCgggg > 1:435409/1‑68 (MQ=255) tGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTGACGCGTTCAACGGAAAGGGCGGGGCCAATg > 2:880980/1‑68 (MQ=255) aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAg > 2:134131/1‑68 (MQ=255) gCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCaa < 1:617593/68‑1 (MQ=255) cTCGAGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAAt > 1:950691/1‑68 (MQ=255) cGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATAc < 2:921532/68‑1 (MQ=255) gCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCACTGCCAGCAATACCtt > 1:730200/1‑68 (MQ=255) gCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCtt > 1:292571/1‑68 (MQ=255) gCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGtt > 2:58088/1‑68 (MQ=255) ccTCCGGTGGTAGCACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTtg > 2:889770/1‑68 (MQ=255) ccTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTtg > 1:699579/1‑68 (MQ=255) ccTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTtg > 2:890891/1‑68 (MQ=255) cTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCTGCACGTtga < 1:883049/68‑1 (MQ=255) cGGTGGTAACCCGTTCACCGGAACGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGCtgacgg < 2:523403/68‑1 (MQ=255) cGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTtgacgg < 2:470216/68‑1 (MQ=255) gtggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCGGCAATACCTTTAGCGGCACGTtgacggtg > 2:896224/1‑68 (MQ=255) ggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTtgacggtgcc > 2:894558/1‑66 (MQ=255) cgcGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCATTAGCGGCACGTTGACGGTGCCGGTTg > 2:44150/3‑68 (MQ=255) cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGCCGGTTg > 2:582551/1‑68 (MQ=255) cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTg > 1:943591/1‑68 (MQ=255) cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTg > 2:826710/1‑68 (MQ=255) ctctCAACGCAAAGGTCGGGCCCAATCCCAGCAATACCTTTAGCGGCACGTTGACGGCGACGGTCGAc < 2:101131/65‑1 (MQ=255) cGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGCCGGTGACGGTTGAc > 1:188469/1‑68 (MQ=255) cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc < 2:943588/68‑1 (MQ=255) cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc < 1:686914/68‑1 (MQ=255) cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc < 1:361471/68‑1 (MQ=255) aaCGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGt < 2:292571/68‑1 (MQ=255) cGGAAAGGGCGGGGCCAATGCCAGCAATACCATTAGCGGCACGTTGACGGTGACGGTTGACCAGGTaa > 1:23416/1‑67 (MQ=255) gCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAAc < 1:732042/68‑1 (MQ=255) | GCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAAC > NC_000913/1135908‑1136037 |
| 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: ‑ |