| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,778,686 | 0 | T | G | 58.8% | 7.1 / 14.5 | 17 | intergenic (‑316/+435) | BW25113_RS25705/BW25113_RS13825 | DUF5507 domain‑containing protein/tRNA‑Ile |
| Reads supporting (aligned to +/- strand): ref base T (2/5); new base G (10/0); total (12/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.39e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-02 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778604‑2778774 | gAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTgg < 2:24713/90‑1 (MQ=255)gAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTgg < 2:529664/90‑1 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGGCAGTAGTAGATGGTTaaag > 2:489330/1‑88 (MQ=255) aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTCATc > 1:55736/1‑90 (MQ=255) tGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGGAGTAGATGTTTTAAGCGGGGGGGAGAAAATTTATTCTTActc > 2:689600/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGGCAGGAGTAGGTGTTTTAGGGGTGGGAGAGAAATTTTATCCTTACTCTTCGGc > 1:303413/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTAAAGCGTGGGAGAGACATTTTATCCTTTCTCTACGGc > 2:245875/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAAATGTTTTAAGGGTGGGGGGGACATTTTATCCTTTCTTTTCGGc > 2:307980/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTACTCTTCGGc > 1:348661/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGTCAGTAGTAGATGTTTTAGGGGTGGGAGAGACATTTTTTCCTTAGTCTACGGc > 1:489109/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 1:173548/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 2:489181/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTTGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTTCTCTACGGCAt > 1:442739/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTTAGGGGGGGCAGAGACATTTTATCCTTACTCTACCGcct > 1:699199/1‑88 (MQ=255) cTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTTAGGGGTGGGAGAGACATTTTATCCTTTCTCTTCGGCCTTGtttttc > 2:232126/1‑86 (MQ=255) cATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacat < 1:232126/90‑1 (MQ=255) cggcggGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCTTACTCTACGGCATTGTTTTACATACATTGGTtg > 2:75834/1‑90 (MQ=255) gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc < 2:699199/90‑1 (MQ=255) | GAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778604‑2778774 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |