Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,953,261 | 0 | T | G | 38.1% | 36.8 / 10.2 | 42 | intergenic (‑88/+111) | ppiC/BW25113_RS25790 | peptidylprolyl isomerase PpiC/type II toxin‑antitoxin system PemK/MazF family toxin |
Reads supporting (aligned to +/- strand): ref base T (6/20); new base G (16/0); total (22/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.77e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATT > NZ_CP009273/3953173‑3953349 | tGATCGTGCCTTGAGTGAATAATTGTCAGGGGGCGATTTTAGCATGATCCACAGATGCTAACGTAATTGTCATTCAGGGAGTGGCTGGTa < 2:934002/90‑1 (MQ=255) tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGTCTGGTa > 2:88836/1‑90 (MQ=255) tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGTCTGGTa < 1:88836/90‑1 (MQ=255) tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTa < 2:820915/90‑1 (MQ=255) tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTa < 2:147085/90‑1 (MQ=255) gTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACatat < 1:1015333/90‑1 (MQ=255) gCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACatatat < 1:123381/90‑1 (MQ=255) tGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCa > 1:895366/1‑90 (MQ=255) gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaa < 2:1036464/90‑1 (MQ=255) gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATGGCaa > 2:220849/1‑90 (MQ=255) gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaa > 1:554721/1‑90 (MQ=255) tGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaaaaa > 2:64370/1‑89 (MQ=255) gAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGc > 1:553098/1‑90 (MQ=255) ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGAATATATGGCAAAAAACCggga > 2:313302/1‑88 (MQ=255) ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCggga > 1:166226/1‑88 (MQ=255) gTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaaa > 1:1053606/1‑90 (MQ=255) cAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATc < 2:323499/90‑1 (MQ=255) cAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATTGCCAAAACCCGGaaaaac > 1:482309/1‑88 (MQ=255) gcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCt < 2:699909/90‑1 (MQ=255) tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAgg < 1:313302/90‑1 (MQ=255) gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTc > 1:360255/1‑90 (MQ=255) gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAAACCGGGCTTCGCAGGTCCCTCAAgggggg > 1:660477/1‑88 (MQ=255) aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga < 2:287116/90‑1 (MQ=255) tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATGGCAAAAACCCGGAAAAAACGGGGTTCGCAGGGCCCTCCAGGGGGGaa > 2:455079/1‑90 (MQ=255) tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGAAATATATGGCAAAAAACCGGGAAAACCGGGGTTCGCCGGGCCCTCCAGGGGAGaa > 2:627630/1‑90 (MQ=255) gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGaaa > 1:405130/1‑90 (MQ=255) gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATGGCAAAAGCCCGGAAAAATCGGGGTTCGCAGGGGCCTCAAGGGGGGaaa > 2:938903/1‑90 (MQ=255) gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGAGAAATTTTGGAAAAGACCGGAAAATCCCGGCTTCGGATGGGCCTCAAGGGGGGaaa > 2:380916/1‑90 (MQ=255) gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCCGAAAATCCGGGCTTCGCAGGGGCCTCAAAGGGAGaaa > 1:78453/1‑90 (MQ=255) tAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCCGGAAAACCGGGGTTCGCAGGGCCCCCCAGGGGGGAAAAc > 2:900547/1‑90 (MQ=255) aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt < 1:777987/90‑1 (MQ=255) gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATTCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg < 1:64370/90‑1 (MQ=255) tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg < 1:152469/90‑1 (MQ=255) tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAATTAgg > 1:410652/1‑90 (MQ=255) aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc < 1:455079/90‑1 (MQ=255) gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc < 1:508511/90‑1 (MQ=255) gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc < 1:900547/90‑1 (MQ=255) gtcagGGAGTGGCTGGGACATATATTGCAAAAACCCGGGAAAACCGGGCTTCGCAGGGCCCTCAAGGGGGGAAAACTTAGGGGCTCTATg > 2:508511/1‑90 (MQ=255) tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg < 2:410652/90‑1 (MQ=255) cagGGAGTGGCTGGTACATATATGGCAAAAGCACGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg < 1:367555/90‑1 (MQ=255) ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt > 1:321725/1‑90 (MQ=255) gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt < 1:917529/90‑1 (MQ=255) | TGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATT > NZ_CP009273/3953173‑3953349 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |