Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,245,989 | 0 | C | T | 44.7% | 36.4 / 22.5 | 47 | E136K (GAA→AAA) | plsB | glycerol‑3‑phosphate 1‑O‑acyltransferase PlsB |
Reads supporting (aligned to +/- strand): ref base C (3/23); new base T (21/0); total (24/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.81e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATA > NZ_CP009273/4245909‑4246073 | gTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTGACCTCGCCTTTTTTTCGcccccg > 2:608933/1‑88 (MQ=255) gcgACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGgcgc < 1:995226/90‑1 (MQ=255) gcgACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGgcgc < 2:129426/90‑1 (MQ=255) gcccGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGgcgcgc < 2:1443/88‑1 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg < 2:447135/90‑1 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg < 1:16498/90‑1 (MQ=255) cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCGACCCAACAACAccc > 2:5658/1‑89 (MQ=255) gCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTCCGCCCCCGCGCGCGACCAAAAAACACCCac > 1:716823/1‑90 (MQ=255) gCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTCCGCCCCGGCGGGCGACCCAAAATCACCGac > 2:719/1‑90 (MQ=255) gAAAATTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac < 1:769216/90‑1 (MQ=255) gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac < 1:125899/90‑1 (MQ=255) gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac < 1:414590/90‑1 (MQ=255) gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac < 1:719/90‑1 (MQ=255) aaaaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCGGCGCGCGACCAAACATCACCCACACt > 1:256199/1‑90 (MQ=255) aaaaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCCCTTCGCCTTTTTCACGCCCCCGCGCGCGACCAAAAATCACCGACACt > 1:81944/1‑90 (MQ=255) aaaaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCCTTTTTACCCCCCCGCGCGCGACCCAACAACACCCacacc > 1:818912/1‑89 (MQ=255) aaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACAACACCCACACCgg > 2:385426/1‑90 (MQ=255) aCGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCCCCCTTTTTCCACCCCCGCGCCCGGCCAAACATCACCCACACTGGCCCCCTCTGc > 2:986714/1‑90 (MQ=255) ccGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCa > 1:341345/1‑79 (MQ=255) ccGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCa < 2:341345/79‑1 (MQ=255) aTACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAAAATCACCGACACTGGCACCCTCTGCACATCCAGAttt > 1:272253/1‑90 (MQ=255) tACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCGGCGCGCGACCAAACATCACCg > 1:452202/1‑60 (MQ=255) tACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACATCACCGACACTGGCACCCTCTGCACATCCCGATTTg > 1:36822/1‑90 (MQ=255) tACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAAAAACACCCACACTGGCACCATCTGCACATCCAGAtttt > 1:463338/1‑89 (MQ=255) tACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg < 2:452202/60‑1 (MQ=255) tACGCAGCGGCGGGGTTACTTCGCCTTTTTTACGCCCCCGCGCGCGGCCAAACATCACCCCCACCGGCACCCTCTTCACATCCCGATTTg > 1:947166/1‑90 (MQ=255) tACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCGGCCCAAAATCACCGACACCGGCCCCATCTGCACATCCCGATTTg > 2:266700/1‑90 (MQ=255) tACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCGACCAAAAATCACCGACACCGGCACCCTCTTCCCATCCCGAtttt > 2:17280/1‑89 (MQ=255) cGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTtgg > 1:290381/1‑90 (MQ=255) cGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAAAATCACCGACACTGGCACCCTCTGCACATCCAGATTTggg > 1:879059/1‑90 (MQ=255) cGCAGCGGCGGGGTCACTTCGCCTTTTTCACGCCCCCGCGCGCGACCAAACATCACCCACACCGGCACCCTCTGCACATCCAGATTTtgg > 1:853520/1‑90 (MQ=255) gcggcgGGGTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAAAATCACCCACACTGGCCCCCTCTTCCCATCCCGATTTGGGGTgg > 2:345347/1‑89 (MQ=255) cggcggGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCGGCCAAAAAACACCCACACTGGCACCCTCTGCACATCCAGATTTTGGGTGCt > 2:198363/1‑90 (MQ=255) cggGTTCACTTCGCCTTTTTCCCGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACg < 2:879059/90‑1 (MQ=255) cggGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACg < 1:122088/90‑1 (MQ=255) cggGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACg < 1:384197/90‑1 (MQ=255) ggggTCACTTCGCCTTTTTTACGCCCCCGCGCGCGGCCCAAAAACACCGACACCGGCACCATCTGCACATCCAGATTTTGGGTGTTACgg > 2:990781/1‑90 (MQ=255) ggTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGt < 1:986714/90‑1 (MQ=255) cGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac < 1:582240/43‑1 (MQ=255) cGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac > 2:582240/1‑43 (MQ=255) ccTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCgaga < 1:17280/90‑1 (MQ=255) tttttCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATa < 2:36822/90‑1 (MQ=255) tttttCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATa < 1:675160/90‑1 (MQ=255) tttttCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATa < 1:522632/90‑1 (MQ=255) tttttCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATa < 1:362764/90‑1 (MQ=255) tttttCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATa < 1:182057/90‑1 (MQ=255) tttttCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATa < 1:119816/90‑1 (MQ=255) | GTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATA > NZ_CP009273/4245909‑4246073 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |