Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,794,496 | 0 | T | G | 58.8% | 13.1 / 16.0 | 17 | V58G (GTG→GGG) | nrdI | class Ib ribonucleoside‑diphosphate reductase assembly flavoprotein NrdI |
Reads supporting (aligned to +/- strand): ref base T (1/6); new base G (10/0); total (11/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.66e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTCTGGTAATCGCAACTTTGG > NZ_CP009273/2794413‑2794583 | cAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACg > 1:117232/1‑90 (MQ=255) gATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGGACGGCTGGGGCGGGGGCACCACA‑GGGAAATCGCtt > 2:99039/1‑90 (MQ=255) aGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACA‑GGTAATTCGCTTTTTAAacga < 1:99765/90‑1 (MQ=255) aGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGACA‑GGTAATTCGCTTTTTAAacga > 1:180703/1‑90 (MQ=255) gCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGACA‑GGGAATTCGCCTTTTAAAAGACGAGca > 1:433614/1‑90 (MQ=255) gcTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACA‑GGTAATTCGCTTTTTAAACGACGAGcac < 1:267389/89‑1 (MQ=255) tACATCGTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGACA‑GGGAATTCGCTTTTTAAACGACGAGCACAAc > 1:330266/1‑90 (MQ=255) tACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACA‑GGGAATTCGCTTTTTAAACGACGAGCACAAc > 1:122920/1‑90 (MQ=255) tGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGACA‑GGGAATTTGCTTTTTAAAAGACGAGCACAACCgggggg > 2:61492/1‑87 (MQ=255) tGATCGTGCCCTCTTACGGCGGCGGCGGGACGGCTGGCGCGGGGCCCCCACA‑GGGAATTTGCTTTTTTAACGACGAGCACAAACGGGGGt > 2:362510/1‑90 (MQ=255) gATCGTGCCCTCTTACGGCGGCGGCGGGACGGCTGGGGCGGGGCCCCGGCA‑GGGAATTTTCTTTTTTAACGGCCAGGCCAACCggggggc > 2:267389/1‑86 (MQ=255) gTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCCCGACAGGGAAATTCCCTTTTT‑AACGACGGGCACCAACGGGGGGTGCtt > 2:377310/1‑90 (MQ=255) tGCCCTCTTACGGCGGCGGCGGTACGGCTGGGGCGGGGCCACGGCA‑GGGAATTTGCTTTTTTAACGACGAGCACAACCGGGCGGTGCTTc > 2:278320/1‑90 (MQ=255) cccTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgc < 2:330500/90‑1 (MQ=255) cccTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgc < 1:101616/90‑1 (MQ=255) cggcggTACGGCTGGCGCGGTGCCACGACA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTc < 1:160107/90‑1 (MQ=255) ggTGCCACGACA‑GGTAATTCGCTTTTCAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTCTGGTAATCGCAACTTTgg < 1:99039/90‑1 (MQ=255) | CAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACA‑GGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTCTGGTAATCGCAACTTTGG > NZ_CP009273/2794413‑2794583 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |