Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,290,243 | 0 | T | G | 54.2% | 20.7 / 32.5 | 24 | L87F (TTA→TTC) | ccmB | heme exporter protein CcmB |
Reads supporting (aligned to +/- strand): ref base T (4/7); new base G (13/0); total (17/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.53e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGC > NZ_CP009273/2290155‑2290330 | ttGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTa < 2:566560/90‑1 (MQ=255) aaGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGCAACAACATCaa > 2:455318/1‑90 (MQ=255) ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAacaaca < 1:133499/73‑1 (MQ=255) ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAacaaca > 2:133499/1‑73 (MQ=255) aGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGCAACAACAACAATTTTTCAAGACTGccc > 2:782539/1‑89 (MQ=255) gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt < 2:320545/74‑1 (MQ=255) gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTc < 2:762160/79‑1 (MQ=255) gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAAt > 1:320545/1‑74 (MQ=255) gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGttt > 1:762160/1‑78 (MQ=255) cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCa > 1:101607/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc > 2:256940/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc > 2:253375/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc > 2:171078/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAAACTGCCGTCCCGCCAATCGTCAc > 2:657463/1‑90 (MQ=255) ccATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAATTGTTCAAGACTGCCGTCCCGCAAAACGTCACGGAACa > 1:300482/1‑90 (MQ=255) aCCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGt > 2:338694/1‑90 (MQ=255) gCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCa > 2:566417/1‑90 (MQ=255) cacaACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTGCCAGCGCCAg > 2:786770/1‑90 (MQ=255) aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa < 2:101607/90‑1 (MQ=255) aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa < 2:32659/90‑1 (MQ=255) aCGGCGGGCCAGGGGAACGGCAACAACATCAATTGTTCAAAACTGCCGTCCTGCAAAACGGCCCGGAACAAGCGTTCCCGCGCCAGCAAg > 1:688264/1‑90 (MQ=255) gggCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGa > 1:243913/1‑90 (MQ=255) ggTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAagcagc > 2:415877/1‑90 (MQ=255) ggTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAagcagc < 1:415877/90‑1 (MQ=255) | TTGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGC > NZ_CP009273/2290155‑2290330 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |