Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,245,973 | 0 | T | G | 52.4% | 16.1 / 25.0 | 21 | N141T (AAC→ACC) | plsB | glycerol‑3‑phosphate 1‑O‑acyltransferase PlsB |
Reads supporting (aligned to +/- strand): ref base T (2/8); new base G (11/0); total (13/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGC > NZ_CP009273/4245885‑4246062 | tGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGtt < 1:549911/90‑1 (MQ=255) tGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGtt < 1:546030/90‑1 (MQ=255) agaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACCTctc > 2:458107/1‑88 (MQ=255) gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTACTTTGCCTTTTTTTCGCCCCCgcgcgc > 2:140672/1‑90 (MQ=255) gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTTGCCCTTTTTTCTCCCCCgcgcgc > 2:401845/1‑90 (MQ=255) cGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGAc < 2:182197/90‑1 (MQ=255) aGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCCCCCTTTTTACGCCCCCGCGCGCGACCa > 2:442514/1‑90 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg < 1:293519/90‑1 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg < 1:187134/90‑1 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg < 1:674111/90‑1 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTGACTTCTCCTTTTTTCCCCCCCCGCGCGCGGCCAAACAACAccc > 2:593810/1‑89 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCCCTTCGCCCTTTTTTCGCCCCCGCGCGCGCCCCAACAACAccc > 2:442944/1‑89 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCCTTTTTTCCCCCCCGCGCGCGGCCAAACAACAccc > 2:588779/1‑89 (MQ=255) gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac < 1:401845/90‑1 (MQ=255) gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTGCCTTGCCCTTTTCCCCCCCCCGCCCGCGACCAAACATCACCCacac > 1:377194/1‑90 (MQ=255) gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCCGCCAAACAACACCgcccc > 2:433956/1‑86 (MQ=255) aaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACGCCCCCGCGCGCGGCCAAACATCACCGACACTgg > 1:277090/1‑90 (MQ=255) aTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCa < 1:458107/90‑1 (MQ=255) aCGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGc > 1:39595/1‑90 (MQ=255) ttAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATcc < 1:269772/90‑1 (MQ=255) aaGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCGGCCCAACAACACCCACACCGGCACCCCCTGCACAccccg > 2:600352/1‑88 (MQ=255) tACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCGGCGCGCGACCAAACATCACCg > 1:406075/1‑60 (MQ=255) tACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg < 2:406075/60‑1 (MQ=255) ttCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGc < 2:77538/90‑1 (MQ=255) | TGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGC > NZ_CP009273/4245885‑4246062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |