Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,612,067 | 0 | A | G | 66.7% | ‑5.4 / 16.7 | 27 | E161G (GAG→GGG) | yjjJ | type II toxin‑antitoxin system HipA family toxin YjjJ |
Reads supporting (aligned to +/- strand): ref base A (3/6); new base G (18/0); total (21/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.84e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGT > NZ_CP009273/4611991‑4612155 | ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGGATTAttaa > 2:647426/1‑87 (MQ=255) cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg < 2:169085/90‑1 (MQ=255) acaggaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 1:707450/86‑1 (MQ=255) tgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCacaac > 1:651447/1‑90 (MQ=255) gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATTAGCGATGGGTTACTGCacaaca > 1:589631/1‑90 (MQ=255) aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATTAGCGATGGATTACTGCACAACACCCTg > 2:616098/1‑90 (MQ=255) aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGGGGGGGGGGGGGATTATAAGCGAGGGAATACTGCACAACAcctgt > 2:707450/1‑87 (MQ=255) cccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGGGGGGGGGGGGATTTTCAGCGATGGATTAGTGGACAACAACCTGCgg > 2:604102/1‑90 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTTGCTGGGCGGGGGGGGGGATTATCAGCGATGGGTTTCTGCACAACACCCCTCGGaa > 1:681370/1‑90 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag < 1:570981/48‑1 (MQ=255) cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATTAGCGATGGGTTTCTGCACAACACCCCGCggga > 2:224612/1‑88 (MQ=255) tGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTTTCAGCGATGGGTTAATGCACAACACCCTGCGGgaa > 2:463804/1‑90 (MQ=255) gTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTCct > 2:589879/1‑90 (MQ=255) aTTTAACGGTGAATATACTGGCGGTTGGTTTGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTTctct > 1:378855/1‑90 (MQ=255) aaCGGTGAATATACTGGCGGTTGGTTTGTCGGGGGGGGGGATTATTAGCGATGGATTACTGCACAAAACCCTGCGGGAATTCCTCTGGAt > 2:725426/1‑90 (MQ=255) aCGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGAAATTTCTCTGGATc > 1:720953/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTCCTCTGGATCaaaaa > 2:651369/1‑90 (MQ=255) gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaaaaa > 1:575143/1‑90 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 1:463804/90‑1 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 1:647426/90‑1 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 1:604102/90‑1 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 2:681370/90‑1 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCCGCGGGAATTCCTCTGGATCAAAAACTCAcccc > 2:599682/1‑89 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCa > 1:118310/1‑90 (MQ=255) tGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTCCTCTGGATCAAAAAATCACCCa > 2:176585/1‑90 (MQ=255) gCGGTTGGTTGGGCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTTCTCTGGATCAAAAAATCACCCAtt > 2:209816/1‑90 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTAc > 2:30894/1‑90 (MQ=255) ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCCGCGGAAATTCCTCTGGATCAAAAACTCACCCATTAc > 1:362599/1‑90 (MQ=255) ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa < 2:362599/90‑1 (MQ=255) gaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGt < 1:30894/90‑1 (MQ=255) | CTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGT > NZ_CP009273/4611991‑4612155 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |