| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,290,243 | 0 | T | G | 61.1% | 25.5 / 41.1 | 37 | L87F (TTA→TTC) | ccmB | heme exporter protein CcmB |
| Reads supporting (aligned to +/- strand): ref base T (1/13); new base G (22/0); total (24/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.06e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGC > NZ_CP009273/2290156‑2290330 | tGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTaa < 2:270938/90‑1 (MQ=255) gCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAAc > 1:225496/1‑90 (MQ=255) gTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAAAGGccaaa > 1:1087572/1‑86 (MQ=255) aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt < 2:294798/90‑1 (MQ=255) aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt < 1:991212/90‑1 (MQ=255) gATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAAAAACAACAAttt > 1:95901/1‑89 (MQ=255) tGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGtt > 2:602412/1‑90 (MQ=255) ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAAAGGCAACAACAACAAATTTTTAAAACTg > 1:297763/1‑90 (MQ=255) ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAAATGTTCAAAACTg > 2:1060273/1‑90 (MQ=255) aGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGCAAAAAAACCAATTGTTCAAAACTGCCg > 1:1106155/1‑90 (MQ=255) gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGt > 1:196769/1‑90 (MQ=255) ccGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTcc < 1:663994/90‑1 (MQ=255) gTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTTAAAACTGCCCGCCTGc > 1:391512/1‑90 (MQ=255) tACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAAAATCAATTGTTCAAGACTGCCGTCCTGCaa > 2:125823/1‑90 (MQ=255) cATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATc < 2:196769/90‑1 (MQ=255) tCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAATTTTTCAAAACCGCCGTCCCGCAAAACGt > 1:981730/1‑90 (MQ=255) cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaa < 2:996631/50‑1 (MQ=255) cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGcaa > 1:996631/1‑50 (MQ=255) cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCCACAACATCAATTGTTCAAAACTGCCGTCCTGCAAAACGACa > 1:308824/1‑90 (MQ=255) cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAAACTGCCGTCCTGCAAATCGTCa > 1:126154/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc > 2:421105/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc > 1:520154/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAATTGTTCAAAACTGCCGTCCTGCAAATCGTCAc > 1:1085991/1‑90 (MQ=255) gTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAGACCGCCGGCCTGCAAAACGTCACg > 1:837328/1‑90 (MQ=255) gTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAAACTGCCGTCCTGCAAATCGTCACg > 1:92898/1‑90 (MQ=255) cgcCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGaa < 2:622658/90‑1 (MQ=255) cgcCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGaa < 2:753382/90‑1 (MQ=255) ccATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGGCCTGCAAAACGTCACGGAACa > 2:819146/1‑90 (MQ=255) aTCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGt < 1:119224/90‑1 (MQ=255) aTCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGt < 2:464302/90‑1 (MQ=255) cGCCAGCACAACGGCGGGCAAAGGAAACCGCAACAACATCAATTTTTCCAAACTGCCCTCCCGCCAAAGAGCACGGGAAAGACGTTCCAg > 1:595978/1‑90 (MQ=255) gCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCa > 2:196710/1‑90 (MQ=255) gCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCa > 1:234221/1‑90 (MQ=255) gCACAACGGCGGGCAAGGGGAACGGCAAAAACAACAATTGTTCAAGACTGCCGTCCTGCAAAACGGCACGGGACAGTCGTTCCAGCGCCa > 1:625894/1‑90 (MQ=255) caACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCa < 2:595978/90‑1 (MQ=255) aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa < 1:233860/90‑1 (MQ=255) ggTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAagcagc < 2:95901/90‑1 (MQ=255) | TGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGC > NZ_CP009273/2290156‑2290330 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |