Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 937,169 | 0 | T | G | 55.3% | 10.1 / 36.9 | 38 | V252G (GTG→GGG) | dmsA | dimethylsulfoxide reductase subunit A |
Reads supporting (aligned to +/- strand): ref base T (1/16); new base G (21/0); total (22/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.64e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCG > NZ_CP009273/937087‑937258 | aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTTTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc < 2:861296/90‑1 (MQ=255) aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc < 2:981393/90‑1 (MQ=255) gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt > 1:1616869/1‑90 (MQ=255) tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa < 2:1415362/90‑1 (MQ=255) tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTAGGGTGACTTACTa < 1:743334/81‑1 (MQ=255) tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTAGGGGGACTTACTa > 2:743334/1‑81 (MQ=255) cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc < 1:150346/90‑1 (MQ=255) gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg < 2:163239/69‑1 (MQ=255) aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc < 1:1022826/90‑1 (MQ=255) cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTTAAAAGGGACAccc > 2:150346/1‑89 (MQ=255) aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAg < 2:1831032/90‑1 (MQ=255) tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTATTATCTTTAACAAGGACGCCAGaa > 1:1002779/1‑90 (MQ=255) tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTGAAAAGGCACCCCaaaa > 2:915803/1‑90 (MQ=255) tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTGAACAGGCACCCCCGaa > 2:860786/1‑90 (MQ=255) ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAACTaa > 2:1862627/1‑90 (MQ=255) gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCGTTCTATCTTTGACAGGGACGCCCAAAAACTAAt > 2:181239/1‑90 (MQ=255) gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg < 2:1002780/90‑1 (MQ=255) cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTGAAAAAGCACGCCAGAAAACTAATGCCCGc > 2:502781/1‑90 (MQ=255) gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt < 1:607649/90‑1 (MQ=255) tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATTTTTAACAGGCACGCCAGAAAACTAATTCCCCCATg > 2:1615822/1‑90 (MQ=255) ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATCTTTAAAAGGGACGCCCGAAAACTAATGCCCGCATGa > 1:1184793/1‑90 (MQ=255) agcccTGTCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc < 2:1218031/88‑1 (MQ=255) aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCCTGatcatcatc > 2:1131430/1‑90 (MQ=255) aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACCCCAGAAATCTAATGCCCGCATGatcatcatc > 2:205786/1‑90 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc < 2:168539/90‑1 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc < 2:1616870/90‑1 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa < 2:563094/39‑1 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAATTAATTCCCCCATGATCCTCatccatcc > 2:276611/1‑90 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCCGAAATCTAATGCCCCAAAGATCATCATCGATcc > 2:516921/1‑90 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc > 2:1722555/1‑90 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTa > 1:563094/1‑39 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATTTTTGAAAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCTATcc > 2:808269/1‑90 (MQ=255) tGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTGAACAGGCACGCCCAAAAACTTATGGCCCCACGGTTATCATCGAttc > 1:331322/1‑88 (MQ=255) cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg < 1:1332623/90‑1 (MQ=255) aCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCCGAAAACTAATGCCCGCCTGGTCATCATCGATCCgcgccct > 1:1218030/1‑87 (MQ=255) cgAATGAGTGGCGGTGGGGGGACCTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCCTGATCCTCATCGATCCGCGCTATAcc > 2:777423/1‑90 (MQ=255) tGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCCGAAAACTAAAGCCCCCCTGATCATCATCGATCCCCGCCATACCCACACCgg > 1:714793/1‑90 (MQ=255) gggggggggACTTTCTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCCCGATCATCATCGATCCGCGCTATACCCACACCCGGGcc > 1:1595644/4‑90 (MQ=255) ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCggg < 2:403140/90‑1 (MQ=255) tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg < 2:823029/90‑1 (MQ=255) | AGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCG > NZ_CP009273/937087‑937258 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |