Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,141,214 | 0 | T | G | 57.9% | 10.9 / 18.4 | 19 | G20G (GGT→GGG) | yegD | molecular chaperone |
Reads supporting (aligned to +/- strand): ref base T (2/6); new base G (11/0); total (13/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATG > NZ_CP009273/2141126‑2141300 | cccTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTa < 1:946071/90‑1 (MQ=255) ttcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGGaaaa > 2:617570/1‑89 (MQ=255) ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCCTTGGTGCCGGGAAAACACCCATGCTaaaaaatga > 1:790159/1‑86 (MQ=255) attTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGaaa < 1:617570/90‑1 (MQ=255) ttATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCa < 2:603086/90‑1 (MQ=255) tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCCTTTTCTAAAAATGGAAAAAGACAGCACGctgctg > 1:63846/1‑90 (MQ=255) cAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAAAAGAAAAAGACAACACGCCGCTGGCTTCAATGCttt > 1:557548/1‑90 (MQ=255) cTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTGCTAAAAATGGGAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGcc > 2:29546/1‑90 (MQ=255) gTTCAGTGGCGGTCATTCGGGGCGGGAAAACGCCTTTTTTAAAAAAAGAAAAAAACAACACGCCGCTGCCTTCCATTCTTTTTGCGCCaa > 1:419497/1‑90 (MQ=255) gTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCATTTGCTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCaa > 2:726172/1‑90 (MQ=255) gTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAAAGGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTGTGCGCCaa > 2:951758/1‑90 (MQ=255) ttCAGTGGCGGTCATGCGTGGCGGGAAAACCCCTTTTTTAAAAAAGGAAAAAAAAAACACGCCGCTGCCCTCCAATCTTTTTGCGCCCAc > 1:925532/1‑90 (MQ=255) ttCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTGCTAAAAAAGGAAAAAGACAGCACCCTGCTGCCTTCAATGCTTTTCGCGCCAAc > 2:93829/1‑90 (MQ=255) tCAGTGGCGGGCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAAAGAAAGCACGCGCCCGCCTTCAATGCTTTTTGCGCCAAcg > 1:484787/1‑90 (MQ=255) gTGGCGGTCATGCGTGGCGGGAAAACGCCTTTTCTAAAAATGGGAAAAAACAGCACGCGGCTGCCCTCCAAGTTTTTTGCGCCAAcgcct > 2:455465/1‑88 (MQ=255) tGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTg < 2:925532/90‑1 (MQ=255) ggCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGa < 2:1050722/90‑1 (MQ=255) ggTCATGCGTGACGGTAAACCGCATTTTCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGTGaaac > 1:769394/1‑88 (MQ=255) tGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAATCGGTAAGCGa < 1:951758/90‑1 (MQ=255) tAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg < 2:673728/87‑1 (MQ=255) tAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg > 1:673728/1‑87 (MQ=255) | CCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATG > NZ_CP009273/2141126‑2141300 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |