Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,781,385 | 0 | T | G | 59.1% | 6.4 / 18.0 | 22 | G396G (GGT→GGG) | envC | murein hydrolase activator EnvC |
Reads supporting (aligned to +/- strand): ref base T (5/4); new base G (13/0); total (18/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.72e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTT > NZ_CP009273/3781307‑3781466 | gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTGAGGGGCGGGCTTcccc > 1:338261/1‑87 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGGCTTCAcc > 1:121141/1‑89 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCAct > 2:91144/1‑90 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCCGGGGCGGCCTTCAct > 1:286536/1‑90 (MQ=255) cgcACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTa > 1:126364/1‑90 (MQ=255) gcACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGCCTTCActcttt > 1:413438/1‑88 (MQ=255) cACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTAtt > 2:155656/1‑90 (MQ=255) gTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGaa > 2:130303/1‑90 (MQ=255) gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTTGAAATTc > 2:58041/1‑90 (MQ=255) ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGGCCTCCCTCTATTTTGAAATTcgccgc > 1:456700/1‑90 (MQ=255) gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc < 1:306054/90‑1 (MQ=255) gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc < 2:338261/90‑1 (MQ=255) ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTcg < 1:131026/68‑1 (MQ=255) ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTcg > 2:131026/1‑68 (MQ=255) ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAAATTGCCGCCCGGGTCAGGCGGTCAAt > 1:43219/1‑90 (MQ=255) cAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCca < 2:283706/90‑1 (MQ=255) ggtggGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCAGGGGCAGGGGGGCAATCCACAGCCGTGGTTggg > 1:116572/1‑90 (MQ=255) gggCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGTGGGTGGGGAg > 2:341488/1‑90 (MQ=255) gcagTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAg > 2:126439/1‑90 (MQ=255) agTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCGATCCCCAGCCGTGGTTGGGGAGATAAGtt > 2:396660/1‑90 (MQ=255) agTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCACAGCCGGGGTTGGGGAGATAAGtt > 2:371527/1‑90 (MQ=255) ggTCAGGGGCGGCCTTCACTCTATTTTGAAATTTGCCGCCAGGGGCAGGGGGGCAATCCACAGCCGGGGGTGGGGAGATAAGGTtttttt > 1:410480/1‑90 (MQ=255) | GAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTT > NZ_CP009273/3781307‑3781466 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |