| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,906,011 | 0 | T | G | 54.5% | 16.1 / 28.0 | 22 | T275P (ACC→CCC) | htpX | protease HtpX |
| Reads supporting (aligned to +/- strand): ref base T (1/9); new base G (12/0); total (13/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.01e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905929‑1906092 | cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGaaa > 2:92315/1‑89 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCGTGAAAAAAACACTTAGCgc > 1:263214/1‑89 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 1:177484/90‑1 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 1:91626/90‑1 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTa > 2:67295/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAACAAACCACTGAGCGATTTCGACTTACCGTTa > 2:224388/1‑90 (MQ=255) aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 2:75370/51‑1 (MQ=255) aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCa > 1:75370/1‑51 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGc < 2:257171/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAACAACTCACCGAGCGATTTCGACTTTCCGGTTATGAAGAGAGcc > 1:349070/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAACAACCCACTGAGCGGTTTTGACTTACCGTTTATGCAGAGAGcc > 1:10556/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGGTTTTGACTTTCCGTTTATGCAGAGAGcc > 1:58442/1‑90 (MQ=255) cAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAAACACTGAGCGGTTTTGACTTACCGTTTATGCAGAGGGCcatcat > 1:375812/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa < 1:263441/39‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc < 2:263214/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc < 2:375812/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGa > 2:263441/1‑39 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc > 2:142699/1‑90 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAAACACTCAGCGATTTTTACTTACCCTTAATGTAGAGAGCCATCATGCTg > 1:358477/1‑90 (MQ=255) cAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTg < 2:58442/90‑1 (MQ=255) gcggcgGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAt > 1:72035/1‑90 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 1:92315/90‑1 (MQ=255) | CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905929‑1906092 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |