Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,906,011 | 0 | T | G | 57.8% | 52.2 / 78.1 | 64 | T275P (ACC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (7/20); new base G (37/0); total (44/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.53e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-09 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCA > NZ_CP009273/1905922‑1906100 | cAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtgggt < 2:1594242/90‑1 (MQ=255) ggACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCa > 1:883771/1‑90 (MQ=255) aCGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATg < 1:1209325/90‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 2:1893788/85‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATg < 2:727134/87‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCa > 1:1893788/1‑85 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATg > 1:727134/1‑87 (MQ=255) aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg < 2:905031/90‑1 (MQ=255) aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg < 2:2714413/90‑1 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGc > 1:1101771/1‑90 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGa > 2:2634925/1‑90 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAGGAACAAATCACTTAGCGa > 2:2567928/1‑90 (MQ=255) cTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATTAAAAACTCACTGAGCGGtt > 2:1768026/1‑90 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAAAAACTCACTTAGCGATTTCg > 2:1769238/1‑90 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 1:1769238/90‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTTAGCGATTTTGACTTAcc > 2:1915071/1‑90 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTCAGCGATTTTGTCTTAcc > 1:1967356/1‑90 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAAACACTGAGCGATTTTGACTTTcc > 1:1619642/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtgggt < 2:2186254/55‑1 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATg < 1:2417370/59‑1 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGt < 2:485435/90‑1 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGt > 2:576443/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGt > 2:527793/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCGTGAACAAAACCCTTAGCGATTTTTACTTACCGt > 1:2648267/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATg > 2:2417370/1‑59 (MQ=255) tCAGGGATTCACCCGTACGCAGGGCTTCAATTCGTTTATCCAGCGGCgggtgggt > 1:2186254/1‑55 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTTACTTACCGTTa > 1:961910/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACCTAGCGATTTCGTCTTTTCGttt > 2:1364088/1‑89 (MQ=255) gTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAAt < 1:1325520/90‑1 (MQ=255) aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAACCACTGGGCGATTTTGACTTACCGTTAATGc > 1:2988401/1‑90 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCa < 2:1111937/90‑1 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAAATCACTGAGCGCTTTCGACTTACCGTTAATGCa > 2:104156/1‑90 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGGTTTTGACTTACCGTTAATGCa > 2:2598460/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 1:666686/90‑1 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCCCTTAGCGGTTTTTACTTTCCGTTAATGGag > 1:2822705/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACCCACCGAGCGATTTTTACTTACCGGTAATGCag > 1:2889073/1‑90 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 1:1731543/90‑1 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGc > 2:713639/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTTGACTTTCCGTTTATGCAGAGAGcc > 2:1530579/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAAAAACTCACTGAGCGATTTCGACTTACCGTTTATGCagagaggc > 2:1291569/1‑88 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc > 2:2373078/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc > 2:537436/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc > 1:729314/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGATTTTGACTTACCCTTTATGCAGAGAGcc > 1:986304/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACCGAGCGATTTTGACTTAACCTTAATTGagagaggc > 1:611104/1‑88 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAAGCACTGAGCGGTTTTTACTTACCGGTAATGGagaggggc > 1:2393076/1‑86 (MQ=255) aCGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcat < 1:2634925/90‑1 (MQ=255) gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatca < 1:2178018/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTTGACTTa > 1:1026488/1‑65 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:1026488/65‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc < 2:2889073/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCTTTAACAACCCCCCGAGCGAGTTTTTCTTTCCGGTCATGCAGAGAGCCATCATGc > 1:1818765/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAAAAACTCACTTAGCGATTTTTACTTTCCCTTAATGCAGAGAGCCAACATGc > 1:1321334/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTACCGTTTATGCAGAGAGCCCTCATGc > 1:937019/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATTCAGAGAGCcctctgtt > 2:2297279/1‑86 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCcaacattc > 1:2821406/1‑88 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGc > 2:2960016/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGc > 2:2797190/1‑90 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGGGGGTGATGGACAACTCACTGAGCGAATTGGGCTTATCGTTAATGAAGAGAGGCATCATGCTg > 2:356810/1‑90 (MQ=255) ttCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGTTTTGACTTACCGTTAATGGAGAGAGCCATCATGCTGGt > 2:1275585/1‑90 (MQ=255) cAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTg < 2:2648267/90‑1 (MQ=255) aTTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTCGACTTACCGTTAATGCAGAGAGGCATCATGCTGGTTGct > 1:886990/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:761236/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:2373078/90‑1 (MQ=255) gcggcgGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCCTTAATGCAGAGAGCCATCATGCTGGTTGCTTTTTTCGGGTCAt > 1:168607/1‑90 (MQ=255) gcggcgGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTTCCCTTAATGCAGAGAGCCCTCATGCTGGGTGCTTCTTGCGGGTCAt > 1:1478405/1‑90 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 2:883771/90‑1 (MQ=255) tCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa < 1:1225797/90‑1 (MQ=255) tCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa < 2:168607/90‑1 (MQ=255) | CAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCA > NZ_CP009273/1905922‑1906100 |
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: ‑ |