Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,088,380 | 0 | A | C | 55.0% | 10.1 / 21.5 | 20 | Y122S (TAC→TCC) | hisH | imidazole glycerol phosphate synthase subunit HisH |
Reads supporting (aligned to +/- strand): ref base A (6/3); new base C (0/11); total (6/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.17e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTCGACTTGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCA > NZ_CP009273/2088291‑2088453 | cGTCGACTTGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTa > 2:691959/1‑90 (MQ=255) cTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCgcaggca < 2:724097/90‑1 (MQ=255) tcGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGcc < 2:155688/90‑1 (MQ=255) cGTGCCGAAAATGCCCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAggg < 2:845547/90‑1 (MQ=255) aaTGACCGCCTTGGGTCTGCCGCTGCCACATATGGGCTGGACCCGCGTTTCCCCGCGGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGa < 2:674197/90‑1 (MQ=255) cgccTTGGGTCTGCCACTGCCACATATGGGCTGGACCCCCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGgcgc < 2:643767/87‑1 (MQ=255) cGACTTTGGCCTGCCCCTCCCACATATGGGCGGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGgcgc < 2:643768/90‑1 (MQ=255) gACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGcgcg > 1:879277/1‑90 (MQ=255) tcTTTGGTCTACCACTCCCACTGTGGGGCAGGATCGCGGTTCCCCCGCGGGCCGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGt < 1:977587/89‑1 (MQ=255) cttttGGCCTGCCACTCCCAAATATGGGATGGAACCGCGTTTCCCCGCGGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGt < 1:557945/88‑1 (MQ=255) gccactgccacATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGtacttttactt > 1:733205/1‑90 (MQ=255) gccactcccacATAGGGGCTGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGtacttttactt < 1:830884/90‑1 (MQ=255) tgccacATATGGGTTGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTc < 1:443941/90‑1 (MQ=255) tgccacATATGGGCTGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTc < 2:314715/90‑1 (MQ=255) tgccacATATGGGCTGGAACCGCGTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTc < 1:461443/90‑1 (MQ=255) tcccacATATGGGCTGAAACCGCTTTTCCCCGCGGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTc < 1:691959/88‑1 (MQ=255) gccacATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTca > 2:67289/1‑90 (MQ=255) ccacATAGGGGCGGGACCCGCTTTTCCCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTcac < 1:491047/90‑1 (MQ=255) caTATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGc > 1:407253/1‑90 (MQ=255) ggCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCa > 1:117124/1‑90 (MQ=255) ggCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCa > 2:35518/1‑90 (MQ=255) | CGTCGACTTGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCA > NZ_CP009273/2088291‑2088453 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |