Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,976,004 | 0 | A | G | 57.1% | 2.2 / 14.0 | 15 | intergenic (‑136/+31) | otsB/araH | trehalose‑phosphatase/arabinose ABC transporter permease AraH |
Reads supporting (aligned to +/- strand): ref base A (2/4); new base G (8/0); total (10/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.50e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGG > NZ_CP009273/1975869‑1976140 | ccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaaaa < 1:24475/139‑1 (MQ=255) cGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATGGTGGGgaaa > 2:487585/1‑138 (MQ=255) tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTGTGGGGGGGAAAAAAGCCCAattatt > 2:147279/1‑139 (MQ=255) ttGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCGGTGGGGGGGAAAAAAGACAAATTGTTGTAGAAAAAAacaac‑a > 2:94643/1‑138 (MQ=255) ttGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTTTTGGGGGGAAAAAAACACAAATGTTGGAGcaaaaaaaaa > 2:254660/1‑132 (MQ=255) tACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGG‑AAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCttt < 1:173036/138‑1 (MQ=255) tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGTCTGTGGTTGGGAAAAAAGCCAAAATTTTGTAGAAAAAAAGAAC‑AGACAGGGCGGTTCGCttttttttttt > 2:193851/1‑133 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGCGATGTTGGGGAAAAAAGCGAAAATGTTGTAGAAAAAAACAACAAAACAAGGGGTTT‑GCTTTTTTCTTGTTAAGGTCGCAGAACa > 1:189094/1‑139 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGGGGGGAAAAAAGCTAAAATGTTGCGGAAAAAAACATA‑AGACAGTGGGGTTAGCTTTTTTGTTTTAAAGGTCGGAGAta > 1:192091/1‑137 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATGGGGGGGAAAAAAGCAAAATTGTTGCAGAAAAAAAAAAC‑AGACAGTGCGGTTCGCGTTTTTTTTGTAAAGGGCGAAGATCa > 1:20651/1‑139 (MQ=255) tgACGGGGGTCGCCATTTGCAACGAAGACCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAcc < 2:127015/139‑1 (MQ=255) gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTTGGTGGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACAAC‑AGACAGGGCGTTTCGCTTTTTTTTTGTTAAGGGCGAAGAACACCGCGGcccg > 1:232402/1‑137 (MQ=255) cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc < 1:254660/139‑1 (MQ=255) tGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa < 1:487585/138‑1 (MQ=255) gtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGc < 2:191174/139‑1 (MQ=255) ggAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGACATATTAAGCAgg > 2:297920/1‑139 (MQ=255) ggAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAgg > 2:528832/1‑139 (MQ=255) | CCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATC‑AGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGG > NZ_CP009273/1975869‑1976140 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |