Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 106,295 | 0 | T | G | 53.3% | 16.7 / 20.7 | 15 | G510G (GGT→GGG) | secA | preprotein translocase subunit SecA |
Reads supporting (aligned to +/- strand): ref base T (1/6); new base G (8/0); total (9/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAA > NZ_CP009273/106209‑106373 | aaaTTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACa < 1:434294/90‑1 (MQ=255) cGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGGCGGGGGGCAGAGAtttg > 1:859031/1‑88 (MQ=255) cAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCt < 2:108817/90‑1 (MQ=255) ttgttgCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTggc < 1:1701/90‑1 (MQ=255) gTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGGCGGGGGGCAGATATTTTGGTCGGGGGGGGCTGGGAGGCAGAAGATGGcg > 2:174662/1‑90 (MQ=255) tATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCcgcg < 1:628571/90‑1 (MQ=255) ccGGCTGCGGTGACTATCGCGACCAATATGGCGGGGGGGGGGGCAGATATTTTGGTCGGGGGGGGGTGGGGGGGAGAAGGTGgcgcgccg > 2:374177/1‑88 (MQ=255) ccGGCTGCGGTGACTATCGCGACCAATATGGCGGGGCGGGGGGCAGAAATTTTGGTCGGGGGGGGGTGGGAGGGAGAAGGTGGCGCGCTg > 2:429998/1‑90 (MQ=255) cTGCGGTGACTATCGCGACCAATATGGCGGGGCGGGGGGCAGAAATTTTGTTCGGGGGGGGGGGGGAGGGAGAAAGTTCCGCGCTGGgaa > 2:399910/1‑90 (MQ=255) aTCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCa < 2:418427/90‑1 (MQ=255) gcgACCAATATGGCGGGTCGTGGGACAGATATTTTGCTCGGGGGGAGCTGGCAGGGAGAAGTTGCCGCGCCGGGAAAACCGACCGCcgag > 2:277719/1‑90 (MQ=255) gcgACCAATATGGCGGGTCGTGGGACAGATATTGTGCTCGGGGGGAGCTGGCAGGCAGAAGGTGCCGCGCTGGGAAAACCGACCGCagag > 2:51042/1‑90 (MQ=255) gcggggcgggGGGCAGATATTGTGCTCGGTGGGAGCTTGCAGGCAGAAGTTGCCGCGCCGGGAAAACCGACCCCCGAGCAAATTTaaaaa > 2:198533/1‑90 (MQ=255) ggTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCcgcg < 1:377853/48‑1 (MQ=255) ggTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCcgcg > 2:377853/1‑48 (MQ=255) | AAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAA > NZ_CP009273/106209‑106373 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |