| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,190,452 | 1 | . | T | 22.2% | 34.1 / 15.3 | 18 | intergenic (‑5/+164) | phoP/purB | response regulator in two‑component regulatory system with PhoQ/adenylosuccinate lyase |
| Reads supporting (aligned to +/- strand): ref base . (9/5); new base T (2/2); total (11/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATATCCGGTATATGTTCATTGAGATAATAATCGGCTTCTTTGGCATCTTCTGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAACCAGCGGGGAGGGGAGGTAAAGTGAAAAAATAAAAAGCGGATAATCTTAATAAG > NC_000913/1190307‑1190576 | aaTATCCGGTATATGTTCATTGAGATAATAATCGGCTTCTTTGGCATCTTCTGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTa > 2:109211/1‑146 (MQ=255)aaTATCCGGTATATGTTCATTGAGATAATAATCGGCTTCTTTGGCATCTTCTGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTa < 1:109211/146‑1 (MQ=255) ggCTTCTTTGGCATCTTCTGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAgt > 3:74314/1‑142 (MQ=255) ggCTTCTTTGGCATCTTCTGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAgt < 4:74314/142‑1 (MQ=255) tcttctGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAAcgcg > 1:26350/1‑137 (MQ=255) tcttctGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAAcgcg < 2:26350/137‑1 (MQ=255) aCCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACAtt > 2:61433/1‑148 (MQ=255) aCCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATta < 1:61433/149‑1 (MQ=255) cAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAAc < 6:96979/146‑1 (MQ=255) cAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAAc > 5:96979/1‑146 (MQ=255) ggtgACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTATTTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATtaaa < 3:95735/115‑1 (MQ=255) ggtgACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTATTTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATtaaa > 4:95735/1‑115 (MQ=255) ggtgACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTATTTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATtaaa > 4:94149/1‑115 (MQ=255) ggtgACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTATTTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATtaaa < 3:94149/115‑1 (MQ=255) aCGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAACCAGCGGGGAGGGGAGGTAAAGTGAAAAAATaaa > 4:87586/1‑148 (MQ=255) tCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAACCAGCGGGGAGGGGAGGTAAAGTGAAAAAATAAAAAGCGGATAATCTtaa > 5:111056/1‑148 (MQ=255) caacaaCCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAACCAGCGGGGAGGGGAGGTAAAGTGAAAAAATAAAAAGCGGATAATCTTAATAAg > 7:49436/1‑148 (MQ=255) aacaacCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAACCAGCGGGGAGGGGAGGTAAAGTGAAAAAATAAAAAGCGGATAATCTTAATAAg > 2:33968/1‑147 (MQ=255) caacCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACa < 4:78520/82‑1 (MQ=255) caacCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACa > 3:78520/1‑82 (MQ=255) | AATATCCGGTATATGTTCATTGAGATAATAATCGGCTTCTTTGGCATCTTCTGCGTCATCGACCTGATGACCAGCATCCTGAATCTGAACTTTAAGGTGGTGACGTAACAACGCATTGTCTTCAACAACCAGTACGCGCATTTTTA‑TTTCTCCCTGTCTTAATTATTAAAATAGTGTAACGCGATTATGTGGTTATGGGGGTAAACATTAAATAAACCAGCGGGGAGGGGAGGTAAAGTGAAAAAATAAAAAGCGGATAATCTTAATAAG > NC_000913/1190307‑1190576 |
| 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: ‑ |