| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,296,060 | 0 | C | T | 33.3% | 61.7 / 52.7 | 54 | intergenic (+266/+376) | gltP/yjcO | glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO |
| Reads supporting (aligned to +/- strand): ref base C (21/15); new base T (6/12); total (27/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.48e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGT > NC_000913/4295983‑4296139 | aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGcc > 1:805065/1‑81 (MQ=33)aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGcc > 1:805066/1‑81 (MQ=33) gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt > 2:805122/1‑80 (MQ=21) ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa < 2:805027/81‑1 (MQ=18) cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt < 1:805028/81‑1 (MQ=18) cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt > 2:805059/1‑81 (MQ=33) gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc < 2:805029/81‑1 (MQ=33) gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc > 2:805086/1‑81 (MQ=33) aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg < 1:805030/81‑1 (MQ=18) aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc < 1:805059/81‑1 (MQ=18) aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc > 1:805084/1‑81 (MQ=33) cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGc < 1:805031/81‑1 (MQ=33) gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt > 2:805119/1‑81 (MQ=18) gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt < 2:805032/81‑1 (MQ=33) tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc < 2:805033/81‑1 (MQ=255) tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc < 2:805035/81‑1 (MQ=38) tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc < 2:805034/81‑1 (MQ=38) tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc > 1:805039/1‑81 (MQ=38) cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt < 2:805036/81‑1 (MQ=255) gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt > 2:805064/1‑81 (MQ=255) gATTAAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt > 2:805120/1‑81 (MQ=255) aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt < 2:805038/81‑1 (MQ=255) tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt > 1:805085/1‑80 (MQ=255) aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt > 1:805133/1‑79 (MQ=255) aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc < 2:805039/80‑1 (MQ=255) gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa < 2:805040/80‑1 (MQ=255) gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt < 2:805076/53‑1 (MQ=255) gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt > 1:805076/1‑53 (MQ=35) aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt < 1:805041/80‑1 (MQ=255) ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct > 2:805123/1‑78 (MQ=255) ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg < 1:805042/80‑1 (MQ=255) tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct > 1:805118/1‑77 (MQ=255) gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt < 2:805084/80‑5 (MQ=255) gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt < 1:805083/80‑5 (MQ=255) gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc > 2:805068/1‑79 (MQ=255) gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt < 2:805078/52‑1 (MQ=255) gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt > 1:805078/1‑52 (MQ=35) ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg > 2:805069/1‑81 (MQ=255) ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc < 1:805043/80‑1 (MQ=255) aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATg < 1:805077/37‑1 (MQ=255) aaGGCGTTTACGCCGCATCCGACATCAACGCCTGATg > 2:805077/1‑37 (MQ=40) ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg < 1:805044/81‑1 (MQ=255) ttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGa > 2:805088/1‑81 (MQ=255) tACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt < 1:805037/47‑1 (MQ=12) tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt > 2:805037/1‑47 (MQ=255) gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa < 2:805046/81‑1 (MQ=38) ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg < 1:805047/81‑1 (MQ=38) gcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg > 1:805082/1‑53 (MQ=40) cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt > 1:805081/1‑39 (MQ=40) aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg < 1:805048/81‑1 (MQ=38) gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAttt > 2:805083/1‑81 (MQ=38) aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg > 1:805090/1‑80 (MQ=39) aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg > 2:805028/1‑81 (MQ=33) aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt > 1:805094/1‑81 (MQ=32) | ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGT > NC_000913/4295983‑4296139 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |