| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,794,496 | 0 | T | G | 53.8% | 7.3 / 12.8 | 13 | V58G (GTG→GGG) | nrdI | class Ib ribonucleoside‑diphosphate reductase assembly flavoprotein NrdI |
| Reads supporting (aligned to +/- strand): ref base T (1/5); new base G (7/0); total (8/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGA‑CA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATT > NZ_CP009273/2794411‑2794560 | cTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCa < 1:301459/90‑1 (MQ=255)cTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGccc > 1:193730/1‑89 (MQ=255) tCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCAc > 1:81348/1‑90 (MQ=255) tGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGA‑Ca < 2:607564/90‑1 (MQ=255) gAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGGACGGCTGGGGGGGGGCCCCGC‑CA‑GGGaaat > 2:42438/1‑88 (MQ=255) ggATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGA‑CA‑GGGAATTCGCt > 1:254850/1‑90 (MQ=255) aGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGA‑CA‑GGTAATTCGCTTTTTAAacga < 1:276753/90‑1 (MQ=255) aGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGGACGGCTGGGGCGGGGGCCCGACCG‑GGGAATTGCCTTTTTAAA‑GACGAGc > 2:276753/1‑90 (MQ=255) ccTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGA‑CA‑GGGAATTCGCTTTTTAAACGACGAGcac > 1:222381/1‑90 (MQ=255) cTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGGCCCGA‑CAGGGAAATTCCCTTTTT‑AACGACGAGGACAACCCGGCg > 2:452450/1‑90 (MQ=255) cTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGA‑CA‑GGGAATTCGCTTTTTAAACGACGAGCACAAACGGGCg > 1:25465/1‑90 (MQ=255) ggcggcggcggTACGGCTGGCGCGGTGCCACGA‑CA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTAtt < 2:170195/90‑1 (MQ=255) ggcggcggcggTACGGCTGGCGCGGTGCCACGA‑CA‑GGTAATTCACTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTAtt < 2:170194/90‑1 (MQ=255) | CTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGA‑CA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATT > NZ_CP009273/2794411‑2794560 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |