| 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.6% | 40.4 / 48.4 | 31 | T275P (ACC→CCC) | htpX | protease HtpX |
| Reads supporting (aligned to +/- strand): ref base T (2/13); new base G (16/0); total (18/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.09e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTT > NZ_CP009273/1905924‑1906098 | ggACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCATACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 1:813656/90‑1 (MQ=255)ggACGCGCTGTTTAGTATTTACTTCATATTCCTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 2:911477/90‑1 (MQ=255) gcgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGaa < 1:1082578/90‑1 (MQ=255) gcgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGaa < 2:1285602/90‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaac > 1:179625/1‑90 (MQ=255) aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg < 1:833882/90‑1 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGc > 2:568192/1‑90 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:658994/90‑1 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGt > 1:297440/1‑90 (MQ=255) cAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACGACTCACTGAGCGATTTTGACTTACCGtt > 1:155037/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGTTTTGACTTACCGTTa > 1:1319773/1‑90 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAACCCCAAAAAGATTTCGACTTACCGttta > 1:296150/1‑88 (MQ=255) aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGc < 2:497707/90‑1 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCa < 1:1338712/90‑1 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 1:286710/90‑1 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGc > 1:405326/1‑90 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGGTTTTGACTTACCGGTAATGGAGAGAGc > 2:129194/1‑90 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTTCCGTTTATGCAGAGAGc > 2:445111/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc < 2:287878/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc > 1:274158/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGTTTCGACTTACCCTTAATGCAGAGAGcc > 1:582910/1‑90 (MQ=255) cAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACCCACCGAGCGATTTCGACTTTCCGTTTATGCCGAAAGCcaccat > 2:1223998/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCCCTGAGCGGTTTTGACTTACCCTTAATGCAGAGAGGCATCATGc > 2:1265805/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCACTGAGCGGTTTTGACTTACCGTTAATGGAGAGAGGCAACATGc > 1:668989/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACACACCGAGCGGTTTTGACTTACCGTTAATGCAGAGAGGCATCATGc > 2:1356503/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGCGGGGCCCGAACAACTCACTCCGCCATTTTGCCTTACCCTTAATGCAGCGAACCCTCAAGc > 1:949911/1‑90 (MQ=255) tCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCACTTAGCGATTTTGCCTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGgttc > 2:100956/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:582910/90‑1 (MQ=255) ttATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTACCGTTTATGCAGAGAGCCATCATGCTGGTTTcttctttc > 1:860388/1‑88 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 2:297440/90‑1 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 1:881024/90‑1 (MQ=255) ggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGtttt > 2:132246/1‑90 (MQ=255) | GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTT > NZ_CP009273/1905924‑1906098 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |