| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,778,686 | 0 | T | G | 61.1% | 7.0 / 17.0 | 18 | 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 (11/0); total (13/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.45e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATT > NZ_CP009273/2778627‑2778763 | ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTTGGTC‑GTAGTAGAT‑GTTTTAggggg > 2:587241/1‑87 (MQ=255) aaTGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTT‑TGGGCTGGAGTAGAT‑GTTTTAGgtgtctct > 2:661755/1‑86 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTT‑TGGTCAGTAGTAGAT‑GTTTTAGGCGTGGCAGAGACa > 1:453766/1‑90 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTT‑TGGGCAGTAGTAGAT‑GTTTTAGGGGGGGCAGGGACATTTTTTCCt > 2:336762/1‑90 (MQ=255) ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGtagtag > 1:191372/1‑55 (MQ=255) ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGtagtag < 2:191372/55‑1 (MQ=255) ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGGGGCAGAGACATTTTATCCtttc > 1:493243/1‑88 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCATTACTCTACGGc > 1:124800/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTT‑TGGTCAGTAGTAGAT‑GTTTAAGGCGGGGGAGAGACAATTTATCCTTTCTCTACggg > 2:652327/1‑89 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTT‑TGGTCAGGAGTAGAT‑GTTTAAGGGGGGGGAGAGACAATTTTTCCTTTCTCTTCGGc > 1:658286/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc > 1:385014/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGGGGCAGAGACATTTTATCCTTACTCTACGGc > 1:550488/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGGTT‑TTGTCAGTAGTAGATGGTTTAAGG‑GTGGCAGAGACATTTCATCCTTTTTCTACGGc > 2:392520/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 1:336762/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 2:658286/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCATTACTCTACGGCAt < 2:124800/90‑1 (MQ=255) aCTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCta < 1:652327/90‑1 (MQ=255) gcctgcGGCGGGGGTTC‑TGGTCAGTAGTAGAT‑GTTTTAGGCGTGGCAGAGACATTTTATCCTTTCTCTACGGGATTGTTTTTCATACAtt > 2:287785/1‑90 (MQ=255) | GGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTC‑TGGTCAGTAGTAGAT‑GTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATT > NZ_CP009273/2778627‑2778763 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |