Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,141,214 | 0 | T | G | 51.7% | 31.7 / 39.0 | 29 | G20G (GGT→GGG) | yegD | molecular chaperone |
Reads supporting (aligned to +/- strand): ref base T (3/11); new base G (15/0); total (18/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATG > NZ_CP009273/2141129‑2141300 | ttcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGTATGCGTGGCGGGaaaa > 1:153354/1‑89 (MQ=255) tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCa > 2:468758/1‑90 (MQ=255) cACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCAt < 1:852020/90‑1 (MQ=255) cGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGCGGGAAAACGCATTTg > 2:444803/1‑90 (MQ=255) ggggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCATGCGTGGCGGGAAAACGCATTTTTTaaaa > 1:280154/1‑90 (MQ=255) ggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAAt < 1:642039/90‑1 (MQ=255) tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCATGCGTGACGGGAAACCGCCTTTTCTTAAAAAGGaa > 1:150382/1‑90 (MQ=255) cGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGGCGGTAAAACGCCTTTTCTAAAAAATGAAAAAGa > 1:412825/1‑90 (MQ=255) cGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTGCTAAAAATGGAAAAAGa > 2:444653/1‑90 (MQ=255) gtgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCTGTGGTGTGGGGGGGGAAAACCCATTTTTTAAAAATGAaaaaacag > 1:1444854/1‑86 (MQ=255) gtgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATTGGTGGCGGTAAAACGCATTTTCTAAAAAAGGGAAAAGAc > 1:102542/1‑90 (MQ=255) ttGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgc < 1:1506430/90‑1 (MQ=255) ttGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgc < 1:191112/90‑1 (MQ=255) ttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAACGACAGCACGCtgctgc > 2:666667/1‑90 (MQ=255) aaaCTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTgcgc < 1:544118/90‑1 (MQ=255) aaCTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTtcgcg > 1:121398/1‑90 (MQ=255) aaCTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTgcgcg > 1:530445/1‑90 (MQ=255) cTGTTCAGTGGCGGTGAAGCGGGTCTGTAAACCGCAATTGTTAAAAATGGGAAAAGACAGCACGCGGCTGCCTTCAATGCTTTTCGCGcc > 2:824359/1‑90 (MQ=255) cTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTTAAAAAGGAAAAAGAAAGCACGCCGCTGCCTTCAATGCTTTTCGCGcc > 1:44203/1‑90 (MQ=255) tGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCa < 2:153354/90‑1 (MQ=255) gTTCAGTGGCGGTCATGCGTGACGGGAAAACGCCTTTTCTAAAAAAGGAAAAAGACAGCCCGCCGCTGCCCTCCATGCTTTGCGCGccca > 2:494366/1‑88 (MQ=255) gTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTGCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGGGCGCCaa > 2:1242321/1‑90 (MQ=255) gTGGCGGTCATGCGTGACGGGAAAACCCATTTTCTAAAAAAAGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTTCGCGCCCAcgcgg > 1:841701/1‑89 (MQ=255) gCGGTCATGCGGGACGGGAAACCGCATTTTCTAAAAATGGAAAAAAACAGCACGCTGCTGCCTTCCATGCTTTTGGCGCCAACGCGGGaa > 2:1393109/1‑90 (MQ=255) ggTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGc < 2:1508245/90‑1 (MQ=255) tCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCgg < 2:1171123/90‑1 (MQ=255) aTGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTa < 2:866186/90‑1 (MQ=255) aTGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTa < 1:723797/90‑1 (MQ=255) cGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg < 1:681982/90‑1 (MQ=255) | TTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATG > NZ_CP009273/2141129‑2141300 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |