| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,906,011 | 0 | T | G | 51.7% | 21.6 / 29.7 | 30 | T275P (ACC→CCC) | htpX | protease HtpX |
| Reads supporting (aligned to +/- strand): ref base T (8/6); new base G (15/0); total (24/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.32e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG > NZ_CP009273/1905922‑1906093 | cAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtgggt < 2:787006/90‑1 (MQ=255) aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg < 1:881874/90‑1 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGc > 1:216942/1‑90 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGc > 1:308782/1‑90 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 2:1027348/66‑1 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCa > 1:1027348/1‑66 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGttttg > 1:1746339/1‑88 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATTAACAACTCACTGGGGGATTTTGTCATACCGt > 1:410056/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAAAAACTCACTGAGCGATTTTTACTTACCgg > 2:1603407/1‑89 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACACACTGAGCGATTTTGACTTACCGt > 2:1677240/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGATTTTGACTTTCCgg > 2:1280105/1‑89 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATTGACAAACCACCGAGCGATTTTGTCTTACCCTTa > 2:1285680/1‑90 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 1:568705/54‑1 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa > 1:379970/1‑90 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACACACTGAGCGATTTCGACTTACCGTTaa > 1:1756102/1‑90 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCa > 2:568705/1‑54 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 1:341129/90‑1 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACCCACTGAGCGGTTTTGACTTACCGGTTAAGGAGAGAGc > 2:213966/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc > 1:1800021/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAAAAAATCACTGAGCGATTTTGACTTTCCGTTAATGCagagaggc > 1:1664356/1‑88 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTCGACTTACCGTTAATGCAGAGAGcc > 1:738349/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc > 1:626062/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc > 1:673186/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCCTCATGc > 1:240072/1‑90 (MQ=255) gTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGcttctt > 1:54425/1‑90 (MQ=255) gTTTATCCAGCGGCGGGGGGGAGATGAACAACCCACTGAGCGGTTTCGAATTACCGTTTATGCAAAGAGCCATTCTTCTTGGTGcttctt > 1:1278500/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:1746339/90‑1 (MQ=255) aTCCAGCGGCGGGTGGGTCATGAACAACTCCCTGAGCGATTTTGCCTTACCGTTAATGCAGAGGGCCATCAAGCTGGTTTCTTCTTTCgg > 1:29148/1‑90 (MQ=255) ccAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGtt < 2:54425/90‑1 (MQ=255) gcggcgGGGGGGGCATGAACAACTCACTTAGCGATTTTGACTTTCCCGTAATGCAGAGAGCCACCATGCTGGTTGCTTTTTTCGGGTcct > 2:1660139/1‑88 (MQ=255) gcggcgGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAt > 2:1035207/1‑90 (MQ=255) gggtggggCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTg > 2:746241/1‑90 (MQ=255) | CAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG > NZ_CP009273/1905922‑1906093 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |