| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,906,007 | 0 | T | G | 53.3% | 13.3 / 18.0 | 15 | H276P (CAC→CCC) | htpX | protease HtpX |
| Reads supporting (aligned to +/- strand): ref base T (1/6); new base G (8/0); total (9/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905943‑1906092 | tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg > 1:339458/1‑66 (MQ=255)tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 2:339458/66‑1 (MQ=255)tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAAAAAATCACCGAGCGa > 2:97942/1‑90 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 1:719302/90‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTAcc > 1:904152/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTa > 1:312245/1‑90 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 1:681876/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGATTTTGACTTTCCGTTAATGCagagaggc > 2:432584/1‑88 (MQ=255) aCGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcat > 1:608910/1‑90 (MQ=255) cAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa < 1:457942/41‑1 (MQ=255) cAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGa > 2:457942/1‑41 (MQ=38) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGc > 1:822966/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGc > 1:250640/1‑90 (MQ=255) tttATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTg > 1:646369/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:881887/90‑1 (MQ=255) tATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCg < 2:667442/90‑1 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 1:97942/90‑1 (MQ=255) | TACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905943‑1906092 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |