Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,290,243 | 0 | T | G | 52.6% | 18.4 / 25.2 | 19 | L87F (TTA→TTC) | ccmB | heme exporter protein CcmB |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (10/0); total (11/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGC > NZ_CP009273/2290157‑2290327 | gCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAc > 1:495046/1‑90 (MQ=255) cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAACGGcaa > 2:256128/1‑90 (MQ=255) cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGcaa > 1:481141/1‑90 (MQ=255) gTGGCGACAGGATGAGCAACGGCAGACCGGTTACCCTCCAGTGCGCCACCCCCTGCGCCAGCACATCGGCGGGCACGGGTAACGGcaacc < 2:158873/90‑2 (MQ=255) aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt < 1:88841/90‑1 (MQ=255) tAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAg < 2:313941/90‑1 (MQ=255) ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaaca < 2:429092/70‑1 (MQ=255) ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGcaaca > 1:429092/1‑70 (MQ=255) ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTg > 1:78505/1‑90 (MQ=255) cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCa > 1:23091/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGGCCTGCAAATCGTCAc > 2:107421/1‑90 (MQ=255) gTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAAACTGCCGTCCTGCAAATCGTCACg > 2:614065/1‑90 (MQ=255) cTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAAACTGCCGTCCTGCAAAACGTCACCGAACAGTCGTTc > 1:464108/1‑90 (MQ=255) cGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAg > 1:435623/1‑90 (MQ=255) gCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCa > 1:64176/1‑90 (MQ=255) aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa < 1:439740/90‑1 (MQ=255) aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa < 2:464108/90‑1 (MQ=255) ggCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGaa < 2:435623/90‑1 (MQ=255) aaGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAagc < 1:324677/90‑1 (MQ=255) | GCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGC > NZ_CP009273/2290157‑2290327 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |