| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,975,994 | 0 | C | G | 53.8% | 10.0 / 15.5 | 14 | intergenic (‑126/+41) | otsB/araH | trehalose‑phosphatase/arabinose ABC transporter permease AraH |
| Reads supporting (aligned to +/- strand): ref base C (1/5); new base G (7/0); total (9/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAC > NZ_CP009273/1975924‑1976066 | aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTTGGGTT‑TGGGGGGaaaaaagct > 1:492642/1‑87 (MQ=255) taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTT‑TGGGGGGGAAAAAAAGTAAAatgttt > 1:736424/1‑89 (MQ=255) taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGGTGGGTT‑TGGGGGGGGAAAAAACTAAAatgttt > 1:352771/1‑89 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGGTGGTGA‑TGGGGGGGGAAAAAGACAAATTTGTGCAGcaa > 1:676222/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGGTGGGTATTGGGGGGAAAAAAACCAAAATT‑TTGTAGaaa > 1:284871/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGGTGGGTA‑TTGGGGGGAAAAAAACTAAAATTTTGTAGaaa > 1:215928/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGGTGGGGA‑TGGGGGGGAAAAAAGCTAAATTGTTTTAGaaa > 1:339471/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGGGGGGGT‑TGGGGGGGAAAAAAACAAAAATTTTGTAGaaa > 1:775367/1‑90 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGTA‑TGGGGGGGAAAAAAGCCAAAATTTTGGAGaaaa > 1:133629/1‑90 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 1:368074/90‑1 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 2:47602/90‑1 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 2:628680/90‑1 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 2:696000/90‑1 (MQ=255) aCCGTTATGGC‑GGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAc < 2:489096/90‑1 (MQ=255) | ATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAC > NZ_CP009273/1975924‑1976066 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |