| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,778,686 | 0 | T | G | 55.6% | 14.5 / 20.3 | 18 | intergenic (‑316/+435) | BW25113_RS25705/BW25113_RS13825 | DUF5507 domain‑containing protein/tRNA‑Ile |
| Reads supporting (aligned to +/- strand): ref base T (1/7); new base G (10/0); total (11/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTG > NZ_CP009273/2778606‑2778764 | aaTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 2:24866/90‑1 (MQ=255)aaTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 2:285704/90‑1 (MQ=255) cAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 2:59905/79‑1 (MQ=255) cAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTc > 1:59905/1‑79 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGGAGATGTTTAAgg > 1:90925/1‑90 (MQ=255) ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCgg > 1:34286/1‑89 (MQ=255) ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGGCAGTAGGAGATGTTTTAGGCgg > 1:200307/1‑89 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCGGGGCAGAGACa > 2:57198/1‑90 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTAAAGCGGGGCAGAGACa > 2:79846/1‑90 (MQ=255) aTAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGTCAGGAGTAGATGTTTTAAGGGGGGGAGAGAAAt > 2:256623/1‑90 (MQ=255) aaCGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTAc < 2:126072/90‑1 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGTGGTTTGGGCAGGAGTAGAGGTTTAAGGGGGGGCCGAGAAATTTCATTCTTACTCTAccgg > 2:414907/1‑87 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTAAGGCGGGGGAGAGGCATTTTATTCTTTCTCTACggg > 1:321352/1‑89 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc > 2:139813/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc > 1:92627/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 1:542848/90‑1 (MQ=255) tgcctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTacatacat < 2:34286/90‑1 (MQ=255) cctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTg < 2:142671/90‑1 (MQ=255) | AATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTG > NZ_CP009273/2778606‑2778764 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |