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.2% | 29.5 / 55.2 | 59 | T275P (ACC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (6/21); new base G (32/0); total (38/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.70e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-09 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCA > NZ_CP009273/1905927‑1906100 | cgcgCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGa > 1:1712091/1‑90 (MQ=255) gcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGaaaa > 2:1294642/1‑88 (MQ=255) tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGa > 1:897162/1‑90 (MQ=255) tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGGAAAAAACACTGa > 2:351339/1‑90 (MQ=255) aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg < 2:1712091/90‑1 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAACAACCCACCGAGc > 2:1552969/1‑90 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGc > 1:548860/1‑90 (MQ=255) ttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCg < 1:929554/90‑1 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCgg > 2:782159/1‑89 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAAAAAATCACTTAGCGAttttg > 1:1684878/1‑88 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCTTGAACAACTCACTGAGCGGttttg > 1:1775373/1‑88 (MQ=255) tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAc < 2:7799/90‑1 (MQ=255) ttACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACtt < 1:782159/90‑1 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:897162/90‑1 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:749236/90‑1 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 1:26472/90‑1 (MQ=255) aCTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAACCACTGAGCGGTTTTGGCTTAc > 1:1219600/1‑90 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCACTGAGCGATTTTTACTTAcc > 2:1638043/1‑90 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAACCACTGAGCGATTTTGACTTAcc > 2:1219941/1‑90 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAAACACTGAGCGGTTTTTACTTAcc > 2:512360/1‑90 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACACACTGAGCGATTTTGACTTAcc > 2:1631481/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGATCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCCCTGAGCGATTTTGACTTAccct > 1:1241497/1‑88 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGt < 1:1049821/90‑1 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCGTGAAAAACCCCCTGAGCGATTTCTACTTAcccc > 2:567892/1‑88 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGATTTTTACTTTccct > 1:673941/1‑88 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAAACCCCGAGCGATTTTTACTTTCCGt > 2:598770/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTa < 1:1638043/90‑1 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGGGGGTGGGTCATGAACAACTCACTGAGGGATTTCGACTTa < 1:1067786/83‑1 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGGGGGTGGGTCATGAACAACTCACTGAGGGATTTCGACTTa > 2:1067786/1‑83 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCCTTaa > 1:512806/1‑90 (MQ=255) gTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAAt < 2:316533/90‑1 (MQ=255) tATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTTAGCGATTTTGACTTTCCGTTAATg > 1:976917/1‑90 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATgga > 1:1325832/1‑88 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACACCCTGAGCGATTTTGACTTTCCGTTAAttga > 2:12214/1‑87 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 1:598770/90‑1 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCAAGAAAAAATCACTGAGCGCTTTTGACTTACCGTTAATGCag > 1:316281/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACCGAGCGGTTTTGACTTACCGGTTATGGag > 2:1049821/1‑90 (MQ=255) aCCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagag < 1:1552969/90‑1 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 1:685704/90‑1 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAACCACCGAGCGATTTTGACTTTCCGTTTATGGagagg > 2:216761/1‑89 (MQ=255) cccGTACGAAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGGACAACCCCCCGAGCGGTTTTTACTTTCCGTTAATGCagaga > 2:84633/1‑90 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGCTGAACAACACACTTAGCGGTTTTGACTTACCCGTAATGCAGAGGGc > 2:1497998/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc < 1:1631481/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc < 1:512360/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAAAAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc > 1:1457851/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCCCTGAGCGATTTTTACTTTCCGTTTATGCAGAGAGcc > 2:1324106/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTTCCCTTAATGCAGAGAGGCATCATGc > 2:685704/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc > 1:1092313/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAACCACCTAGCGGTTTTGACTTTCCGTTAATGCAGAGGGCcatcattc > 2:381132/1‑88 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTTCCGTTAATGCAGAGAGCCATCATGCTg > 1:1369148/1‑90 (MQ=255) cTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTAAGCGGTTTCGACTTACCGTTAATGGAGAGAGCCATCATGCTgg > 1:461209/1‑90 (MQ=255) cAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTg < 1:426333/90‑1 (MQ=255) cAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTg < 2:384961/90‑1 (MQ=255) gTTTATCCAGCGGCGGGGGGGGCATGAACAACTCCCTGAGCGCTTTTGACTTACCCGTAATGCAGAGAGCCATCCTTCTGGTTGcttctt > 1:746172/1‑90 (MQ=255) gTTTATCCAGCGGCGGGGGGGGCACGAAAAACTCACTCAGCGGGTTTGACTTTCCCTTAATGCAGAGGGGCCTCATGCTGGGTGgttctt > 2:18999/1‑90 (MQ=255) tttATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTTGGCTTACCGTTAATGCAGAGAGCCATCATGCTGGGTGCTTTTTg > 1:1551797/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTga < 2:316281/90‑2 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:1462588/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:1324106/90‑1 (MQ=255) tCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa > 1:1366170/1‑90 (MQ=255) | CGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCA > NZ_CP009273/1905927‑1906100 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |