Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,906,011 | 0 | T | G | 59.1% | 10.4 / 24.3 | 22 | T275P (ACC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (13/0); total (14/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.81e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAT > NZ_CP009273/1905929‑1906088 | cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGaaa > 1:597674/1‑89 (MQ=255) tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGa > 2:839272/1‑90 (MQ=255) aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg < 1:16087/90‑1 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACCCACTTAGCGa > 1:310193/1‑90 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:501394/90‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGGTTTTTACTTAcc > 1:904995/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGATTTTGACTTACCGt > 2:655138/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGTTTTGACTTACCGTTa > 1:633568/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 2:75154/90‑1 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATAAACAAATCACTTAGCGATTTCTCCTTACCCTTAAAGCag > 2:478078/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTTATGAAAAACCCACAGAGGGATTTTGCCTTTCCGTTAATGCagagaggc > 2:490563/1‑88 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc > 2:100681/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc < 2:633568/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc > 2:859597/1‑90 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACACCCTCAGCGATTTTGACTTACCCTTAATGCAGAGAGGCATCATGCTg > 1:521744/1‑90 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTCAGCGGTTTTGCCTTTCCCTTTATGCAGAGAGCCATCATGCTg > 1:825051/1‑90 (MQ=255) tCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGCGCCCACATGCTGGTTGcttc > 1:372808/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:934275/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:69073/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:324991/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:409581/90‑1 (MQ=255) gcggcgGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTTCCCTTTATGCAGAGAGCCATCATGCTGGGTGCTTCTTTCGGTTcct > 1:336495/1‑88 (MQ=255) | CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAT > NZ_CP009273/1905929‑1906088 |
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: ‑ |