Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,953,261 | 0 | T | G | 47.8% | 15.7 / 11.0 | 23 | 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 (2/10); new base G (11/0); total (13/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.73e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAG > NZ_CP009273/3953184‑3953344 | tGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATggga > 1:164465/1‑88 (MQ=255) tGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCa > 2:346941/1‑90 (MQ=255) tGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCa > 1:216840/1‑90 (MQ=255) gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATTGCaa > 1:342938/1‑90 (MQ=255) gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATTGCaa > 1:446874/1‑90 (MQ=255) gAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATTGGaaaaac > 1:156163/1‑88 (MQ=255) ggcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGGAAAAACCCCGGAAAACCgggg > 1:234230/1‑89 (MQ=255) gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTc < 2:517418/90‑1 (MQ=255) aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCa < 1:126461/90‑1 (MQ=255) ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg < 1:126745/90‑1 (MQ=255) tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa < 2:71810/90‑1 (MQ=255) gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa > 1:529143/1‑90 (MQ=255) aaTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg < 2:216840/90‑1 (MQ=255) aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc < 2:234230/90‑1 (MQ=255) ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc < 2:144484/90‑1 (MQ=255) cagtcagGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAAACCCGGGTTCGCAGGGCCCTCAAGGGGGGGAAAAATAGGGGCTCTa > 1:240317/1‑90 (MQ=255) agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga > 1:176346/1‑73 (MQ=255) agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga < 2:176346/73‑1 (MQ=255) gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg < 1:346941/90‑1 (MQ=255) ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAGATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa < 2:18859/90‑1 (MQ=255) gTGGCTGGGGCATATATTGCAAAAACCCGCGAAAACCGGGCCTCGCAGGTCCCTCAAGGGGGGAAAAATTATGGCCCCTATCGGCTACaa > 1:539114/1‑90 (MQ=255) gTGGCTGGGACAAATATAGCAAAAACCCGGAAAATCCCGGCTTCGCAGGTCCCTCCAGGGGGGAAAAATTAGGGGCTCTATGGGCTACaa > 1:103663/1‑90 (MQ=255) ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg > 2:524066/1‑90 (MQ=255) | TGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAG > NZ_CP009273/3953184‑3953344 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |