Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,269,399 | 0 | T | G | 53.3% | 12.6 / 18.1 | 15 | G331G (GGT→GGG) | ghxP | guanine/hypoxanthine transporter GhxP |
Reads supporting (aligned to +/- strand): ref base T (1/6); new base G (8/0); total (9/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGC > NZ_CP009273/4269273‑4269520 | tGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGG‑GGGG‑GGGAAAAACGGgttt > 1:1751/1‑138 (MQ=255) atcatcGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGC‑GGGG‑GGGGAAAACGGTGTGACGGcga > 1:93426/1‑137 (MQ=255) tGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGGCGGCGGG‑GGGCGGGGAAAAACGG‑TTGGCGGCCATC‑ACC‑GGTGGc > 2:115731/1‑139 (MQ=255) agcagcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGTGACGGCGGC‑GGGG‑GGGAAAACACGGTTGTCGGCGATT‑ACC‑GCTTGGGTGCTTTTGCTCCTGATTCTGTTCctct > 1:357327/1‑139 (MQ=255) gcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGC‑GGGG‑GGGGAAACCCGGTTGTCGGCTATCAACC‑GTTGG‑GGGGTGGTTCTCCTTGTTTTTTTTCTCTCTcc > 2:412049/1‑139 (MQ=255) cGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGGG‑GGGAAAACCGGTTTGA‑GGCGAAC‑ACC‑GTTGGTGGGCTGGTTCTCCCGGTTTTTTTTCTCTCTccc > 1:57669/1‑138 (MQ=255) cGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCg < 2:99164/139‑1 (MQ=255) gTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGCGGGGG‑GGTAAAACGGGTTTGAGGGC‑ATT‑ACCGGTGGGGGGGGTGTTCTCCCTG‑TTTTTTTTCTCTCTCCGc > 2:425442/1‑139 (MQ=255) ggCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGC‑GGGG‑GGGAAAACCGGGTTGACGGCTCTC‑ACC‑CTTGGGGGGCTGTTCCTCCTGATTCTGTTTCTCTCTCCGCTCTCTTAcc > 1:389218/1‑139 (MQ=255) ggCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGCGGC‑GGGC‑GGGGAAAACGGGTTGGAGGCCGTC‑CCC‑CGTGGGGTGCGGTTTCTCCTGGTTCTGTTTCTCTCTCCCCTCTCTTAcc > 2:149149/1‑139 (MQ=255) ggtgCGGCTCCGGCAGCGTTATACATCGAGTCTGCGGCGGGTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCg < 1:15124/138‑1 (MQ=255) gCGGTATACATCGAGTCTGCGGCGGGTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTc < 1:241351/139‑1 (MQ=255) tCGAGTCTGCGGCGGGTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGAt < 1:115731/139‑1 (MQ=255) ggTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCctgc < 2:437965/139‑1 (MQ=255) ggTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCctgc < 2:54533/139‑1 (MQ=255) | TGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGC‑GGGC‑GGTAAAACCGGTTTGACGGCTATC‑ACC‑GTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGC > NZ_CP009273/4269273‑4269520 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |