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.3% | 34.1 / 44.1 | 35 | T275P (ACC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (4/12); new base G (19/0); total (23/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.18e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCA > NZ_CP009273/1905924‑1906100 | ggACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 1:876253/90‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCTGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaaa > 2:421423/1‑89 (MQ=255) gcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGaaca < 1:911128/90‑1 (MQ=255) aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGTGGGTGGGTCATGAACAACTCACTGAg < 1:1057093/90‑1 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCCTGAACAACACACTCAGc > 2:778345/1‑90 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGc > 1:418368/1‑90 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACACACTGAGCGa > 2:483709/1‑90 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACCCAGCGAttttt > 2:1262647/1‑88 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAAACACTGAGCGAttttg > 2:157184/1‑88 (MQ=255) tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTGATGAACAACCCACCGAGGGGTTTTGAc > 1:566351/1‑90 (MQ=255) atTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACt < 1:258591/90‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTTGACTTAcc > 1:449614/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGGATTAACAAATCACCGAGCGATTTTTACTTTCCGt > 2:923982/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAAAAAACCACTGAGCGGTTTTTACTTTccct > 2:1060265/1‑88 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTTAGCGATTTTGACTTACCGt > 2:1082/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTTCCGGTa > 2:970082/1‑90 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAAATCACTTAGCGCTTTTGACTTTCCGTTaa > 2:1265440/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 2:1042096/90‑1 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCACTGAGCGGTTTTGACTTTCCGGTTATGCag > 2:416810/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCAATGGGCGAATTGGGCTTACCGGTAATACag > 1:451067/1‑90 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 1:1277806/90‑1 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTTGTGAACAACTCACCGAGCGGTTTCGGCTTAACGTTAATTgagaga > 1:340295/1‑90 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGc < 2:828345/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAAAAAACCCCTGTGCGATTTTCTCTTACCCTTAATTCAGGGAGcc > 2:1094346/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACCGAGCGGTTTTGACTTACCGTTAATGCAGAGAGcc > 2:717603/1‑90 (MQ=255) cGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatc < 1:845277/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGATCAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc < 2:438921/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc > 1:372084/1‑90 (MQ=255) gTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTTAGCGATTTCGACTTACCGTTAATTCAGAGAGCCATCATGCTGGTTGcttctt > 1:1288941/1‑90 (MQ=255) gTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTCGACTTTCCGTTAATGCAGAGAGCCATCATGCTGGGTGcttctt > 2:206248/1‑90 (MQ=255) aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 2:15894/61‑1 (MQ=255) aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag > 1:15894/1‑61 (MQ=255) cAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTc < 2:205948/90‑1 (MQ=255) gcggcgGGGGGGGCATGAACAACACACTGAGCGATTTTGGCTTTACGTTAATGCAGAGGGGCACCATGCTGGTTGGTTCTTGCGGTTCAt > 2:823702/1‑90 (MQ=255) cggcggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATa < 1:878058/90‑1 (MQ=255) ggggggggCCTGAACAACTCACTGAGCGATTTCGACTTACCGGTAATGCAGAGAGCCATCATGCCGGGTGCTTTTTTCGGGTCATATCTg > 2:1232350/1‑90 (MQ=255) tCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa < 1:610524/90‑1 (MQ=255) | GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCA > NZ_CP009273/1905924‑1906100 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |