| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 937,169 | 0 | T | G | 44.8% | 22.8 / 17.4 | 29 | V252G (GTG→GGG) | dmsA | dimethylsulfoxide reductase subunit A |
| Reads supporting (aligned to +/- strand): ref base T (3/13); new base G (13/0); total (16/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.13e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGC > NZ_CP009273/937081‑937233 | ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg < 1:809694/90‑1 (MQ=255) cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa > 1:999386/1‑89 (MQ=255) cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa < 2:999386/89‑1 (MQ=255) aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc < 1:513349/90‑1 (MQ=255) aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGc > 1:95965/1‑90 (MQ=255) aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAGGCACGCCAGaaa > 2:961531/1‑90 (MQ=255) ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTTTCTTTAACAGGCACGCCCGAAAACTaa > 1:465351/1‑90 (MQ=255) gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAGGGACGCCAGAAAACTAAt > 1:805182/1‑90 (MQ=255) aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAGGCACGCCAGAAATCTAATGc > 1:237481/1‑90 (MQ=255) gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg < 2:666099/90‑1 (MQ=255) tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa < 2:643950/90‑1 (MQ=255) tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCATg > 1:415664/1‑90 (MQ=255) ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa < 2:88373/59‑1 (MQ=255) ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGaa > 1:88373/1‑59 (MQ=255) aaCCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGGTTTCTATCTTTCAAAAGCACGCCAGAAAAGTAATGGCCGCACGgtcagcatc > 2:812924/1‑90 (MQ=255) cttgCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc < 2:97172/87‑1 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa < 2:774101/61‑1 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATg < 2:61096/69‑1 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc < 2:465351/90‑1 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAACAAGCACGCCCGAAAACTAATGCCCGCCTGATCATCATCGATc > 2:348756/1‑90 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGaa > 1:774101/1‑61 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATg > 1:61096/1‑69 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCAAAATTTAATGCCCCCATGATCTTCATCGATcc > 1:97172/1‑90 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc > 1:248402/1‑90 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTTAACAGGCACCCCCGAAAATTAATGCCCCCCGGATCATCATCggttc > 1:611146/1‑86 (MQ=255) ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg < 2:611146/90‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc < 2:805182/90‑1 (MQ=255) aCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc < 2:272147/72‑1 (MQ=255) aCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc > 1:272147/1‑72 (MQ=255) | GGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGC > NZ_CP009273/937081‑937233 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |