| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,088,380 | 0 | A | C | 56.2% | 11.6 / 15.5 | 16 | Y122S (TAC→TCC) | hisH | imidazole glycerol phosphate synthase subunit HisH |
| Reads supporting (aligned to +/- strand): ref base A (5/2); new base C (0/9); total (5/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.81e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCA > NZ_CP009273/2088299‑2088463 | tGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCgcagg > 2:297316/1‑90 (MQ=255) tcGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGcc < 2:8396/90‑1 (MQ=255) aaGAGGTGCAAAAAATGCCCGACTTTGGTCTGCCCCTCCCCCTTAGGGGCGGGAACCGAGTTTCCCCGAAGGCAGGCAACCGCCTGTTTc < 1:265823/90‑1 (MQ=255) cGTGCCAAAATGGCCCGCCTTTGGTCTGCCACTGCCACATATGGGCGGGACCCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAggg < 1:234708/90‑1 (MQ=255) tGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGa > 1:422726/1‑90 (MQ=255) cttttGTTCTGCCACTCCCACATAGGGGCTGGAACCGCGTTTCCCCGCAGGAAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGcgcgc < 1:495326/88‑2 (MQ=255) aCTTTGGTCTGCCACTGCCCCATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGt < 2:9933/90‑1 (MQ=255) ggTCTGCCCCTGCCACATATGGGGTGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGtacttt < 1:295075/90‑1 (MQ=255) gcccctgccacATATGGGCGGAACCCGCTTTTCCCCGCGGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGtacttttactt < 1:522936/90‑1 (MQ=255) tgccacATATGGGCTGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTc < 2:191150/90‑1 (MQ=255) tgccaaATAGGGGCTGGAACCGCTTTTCCCCGCGGGAAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTc < 1:419828/90‑1 (MQ=255) ggCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCa > 1:479505/1‑90 (MQ=255) gggACCCGCGTTTCCCCGCAGGCGGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATg < 2:420150/89‑1 (MQ=255) gAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGcc > 2:358167/1‑90 (MQ=255) gcgTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCa < 2:373938/90‑1 (MQ=255) gcgTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCa > 1:176264/1‑90 (MQ=255) | TGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCA > NZ_CP009273/2088299‑2088463 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |