| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,180,588 | 0 | A | G | 60.0% | 12.7 / 17.8 | 20 | K165R (AAA→AGA) | yegU | ADP‑ribosylglycohydrolase family protein |
| Reads supporting (aligned to +/- strand): ref base A (4/4); new base G (8/4); total (12/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.48e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCA > NC_000913/2180472‑2180719 | caacGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGgc > 1:52024/1‑139 (MQ=255) aacGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGgcg < 1:704406/139‑1 (MQ=255) cGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGcgcg < 1:622967/139‑1 (MQ=255) ggCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCgg < 2:62688/139‑1 (MQ=255) gCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGt < 1:128925/139‑1 (MQ=255) gCGTGACCAACGGCGCAGCGATGCCCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATGCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGt > 2:572079/1‑139 (MQ=255) cGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCa > 1:23600/1‑139 (MQ=255) tttGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACg < 2:150886/131‑1 (MQ=255) tttGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACg > 1:150886/1‑131 (MQ=255) ttgcTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATATGGCGGTTGCAGGCGCGGTAGTCATCGCATGGTCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTg > 2:414205/4‑139 (MQ=255) tgatgaTGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAAtt > 2:176403/1‑139 (MQ=255) atgatgTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTgc < 2:441717/139‑1 (MQ=255) atgTAGCGCTGGCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGcgcg > 2:569256/1‑139 (MQ=255) tAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGAt > 2:159118/1‑109 (MQ=255) tAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGAt < 1:159118/109‑1 (MQ=255) gCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGc > 1:444248/1‑139 (MQ=255) gCGTCCAGCCCGACCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGc > 2:465803/1‑139 (MQ=255) ccGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCa < 1:176403/139‑1 (MQ=255) aCCCATAGATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCa > 2:59296/1‑139 (MQ=255) aCCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCa > 1:482696/1‑139 (MQ=255) | CAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTCCTTTATTGATGATGTAGCGCTGGCGTCCAGCCCGACCCATAAATCCGATCTGGCGGTTGCAGGCGCGGTAGTCATCGCATGGGCGATTTCTCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCA > NC_000913/2180472‑2180719 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |