Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,975,994 | 0 | C | G | 58.8% | 8.9 / 20.8 | 18 | intergenic (‑126/+41) | otsB/araH | trehalose‑phosphatase/arabinose ABC transporter permease AraH |
Reads supporting (aligned to +/- strand): ref base C (0/7); new base G (10/0); total (10/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGT > NZ_CP009273/1975921‑1976069 | gACATGTAAGTAAAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGTTA‑TGGTGGGAAAAAAc < 2:86073/90‑1 (MQ=255) gACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTA‑TGGGGGGGAAAAAAACTAAAATGTTGCa > 1:206988/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTG‑TGGGGGGGAAAAAAACTAAATTTTTGGTGaaa > 2:154831/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTA‑TGGGGGGGGAAAAAACTAAAATTTTGTAGaaa > 2:30051/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTA‑TGGGGGGGAAAAAAACTAAAATTTTTGAGaaa > 2:488989/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTA‑TGGGGGGGAAAAAAACTAAAATTTTGCAGaaa > 2:60052/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTA‑GGGTGGGGGAAAAAGGGAAAATGGTGGAGaaa > 2:561847/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTT‑TGGGGGGGGAAAAAACTAAAATTTTTCAGaaa > 1:120360/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTA‑TGGTGGGGGAAAAAAGAAAAATGTTGTAGaaa > 1:361657/1‑90 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATTGGGGGGAAAAAAAACAAAATTTTTGAGGaaa > 2:691298/1‑90 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 1:569843/90‑1 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 2:401906/90‑1 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 2:567961/90‑1 (MQ=255) tgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTA‑TGGGGGGGAAAAAAACCAAAATTTTTCAGaaaaaa > 2:569843/1‑90 (MQ=255) agaaCAGCCAACCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCtttt < 1:131428/90‑1 (MQ=255) aaCCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTaa < 1:356152/90‑1 (MQ=255) aCCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAc < 1:497301/90‑1 (MQ=255) gTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGt < 1:101602/90‑1 (MQ=255) | GACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTA‑TGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGT > NZ_CP009273/1975921‑1976069 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |