Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 700,026 | 0 | T | G | 54.5% | 16.1 / 24.3 | 22 | G209G (GGT→GGG) | nagE | PTS N‑acetyl glucosamine transporter subunit IIABC |
Reads supporting (aligned to +/- strand): ref base T (5/5); new base G (12/0); total (17/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.57e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑C‑GGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCC > NZ_CP009273/699937‑700103 | atcaACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGt < 2:536863/90‑1 (MQ=255) gTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑C‑GGgtt > 2:210302/1‑90 (MQ=255) ctgATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑C‑GGTTTTCCAc < 1:494563/90‑1 (MQ=255) gATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGGTTTTCACgg > 2:494563/1‑90 (MQ=255) cTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑C‑gggt > 1:864185/1‑67 (MQ=255) cTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑C‑GGtt < 2:864185/69‑1 (MQ=255) aGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGTTTTTCACGGTGACATTAACCGCTTCTATg > 1:358857/1‑90 (MQ=255) aGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGGTTTCCACGGGGACATTTACCGCTTCTATg > 1:835919/1‑90 (MQ=255) aGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGGTTTCCACGGGGACATTAACCGCTTCTATg > 2:311155/1‑90 (MQ=255) ggTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑c > 1:549618/1‑57 (MQ=255) ggTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑c < 2:549618/57‑1 (MQ=255) cTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCGCCAACGCGGCGGGGACC‑GGTTTTCC‑CGGGGAAATTTACCGCTTCTATTCCggg > 1:876854/1‑89 (MQ=255) cTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑CGGGTTTTCCCCGG‑GACATTAACCCCTTCTATGCCggg > 1:666597/1‑89 (MQ=255) cTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGGTTTCCACCGGGACATTAACCGCTTCTATGCCggg > 2:691078/1‑89 (MQ=255) acacCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGGTTTCCACCGGGACATTAACCGCCTCTATGCCGGGGACg > 1:887726/1‑90 (MQ=255) tCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGGTTTCCACGGTGACATTAACCGCTTCTATGCCGGGGACGGCACcg > 2:854368/1‑90 (MQ=255) gCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑G‑GGTTTTCCCCGGGGACATTAACCGCCTCTTTTTCGGTGGCGGCTCcacg > 2:925546/1‑90 (MQ=21) ggTTCCAGATTGGTGAATTCACCCACGCGGCGGGTA‑C‑GGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGa < 1:250569/90‑1 (MQ=255) gTTCCAGATTGGTGAATTCACCAACGCGGCGGGGA‑C‑GGGTTTCCCCGGGGACATTAACCGCTTCTATGCCGGGGACGGCACCCCgggggt > 2:581216/1‑88 (MQ=255) ggTGAATTCACCAACGCGGCGGGTA‑C‑GGGTTTCCACGGGGACATTAACCGCTTCTATGCCGGGGACGGCACCCCGGGGGTGTTCATGTcc > 2:792101/1‑90 (MQ=255) aTTCACCAACGCGGCGGGTA‑C‑GGTTTTCCACGGTGACATTAACCGCTTCTATTCCGGAGACGGCACCGCGGGGAAGTTCCTGTTCAGctt > 2:818481/1‑90 (MQ=255) tCACCAACGCGGCGGGGA‑C‑GGGTTTCCACCGGGACATTAAACGCTTCTATGCCGGGGGCGGGACCGCGGGGGTGTTTATGTCCGGcttct > 1:411434/1‑90 (MQ=255) cAACGCGGCGGGTA‑C‑GGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTccc < 1:581216/90‑1 (MQ=255) | ATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTA‑C‑GGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCC > NZ_CP009273/699937‑700103 |
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: ‑ |