Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,065,625 | C→T | 11.7% | intergenic (+24/+55) | rpoD → / ← ygjF | RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,065,625 | 0 | C | T | 11.7% | 87.4 / 3.2 | 34 | intergenic (+24/+55) | rpoD/ygjF | RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase |
Reads supporting (aligned to +/- strand): ref base C (10/20); new base T (2/2); total (12/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.02e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCA > minE/2065555‑2065693 | gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAc < 1:1513807/71‑1 (MQ=255) gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAc < 1:1708332/71‑1 (MQ=255) gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGACCAGGCGTTAc < 1:3806964/71‑1 (MQ=255) tCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcg < 1:814064/71‑1 (MQ=255) aCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgcc > 1:1999746/1‑71 (MQ=255) aCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc < 1:734280/70‑1 (MQ=255) ccGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgccg > 1:2733620/1‑70 (MQ=255) cGAGCCTTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa < 1:2153378/71‑1 (MQ=255) gAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCa > 1:2945663/1‑70 (MQ=255) ccGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCAGGATCAGGCGTTACGCCGCAcc < 1:520334/69‑1 (MQ=255) gTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAc < 1:156026/59‑1 (MQ=255) cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCAc < 1:3844102/70‑1 (MQ=255) cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCAc < 1:3643103/70‑1 (MQ=255) aaGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACt < 1:945741/68‑1 (MQ=255) aaGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACt < 1:1825606/68‑1 (MQ=255) ctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCct < 1:2598069/70‑1 (MQ=255) ctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc > 1:1226603/1‑71 (MQ=255) gcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctct > 1:3203436/1‑70 (MQ=255) gCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaa < 1:3166237/71‑1 (MQ=255) cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcg < 1:1345343/41‑1 (MQ=255) cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaa < 1:1229329/70‑1 (MQ=255) cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa > 1:3582902/1‑42 (MQ=255) cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa > 1:483004/1‑42 (MQ=255) cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGca > 1:1682844/1‑63 (MQ=255) tGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc < 1:1477831/37‑1 (MQ=255) gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCAcc < 1:1116410/71‑1 (MQ=255) gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACCCCAcc < 1:1591364/71‑1 (MQ=255) cGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtt < 1:851146/71‑1 (MQ=255) gATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt < 1:623467/71‑1 (MQ=255) ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtttt < 1:3051385/70‑1 (MQ=255) ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtttt < 1:264951/70‑1 (MQ=255) aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCAc > 1:1776422/1‑63 (MQ=255) cGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc < 1:3110541/46‑1 (MQ=255) tAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTg < 1:504840/45‑1 (MQ=255) aGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTg > 1:614746/1‑66 (MQ=255) ggCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt < 1:1505723/59‑1 (MQ=255) ggCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtt < 1:107689/70‑1 (MQ=255) gATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGttttt < 1:2218672/68‑1 (MQ=255) aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGcc > 1:3742418/1‑71 (MQ=255) ggCGTTACGCCGCACCCGGCACTAGGCGCTCTGCACAAACGCCACCTTTTCgg < 1:1837046/53‑1 (MQ=255) cGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt < 1:2461272/47‑1 (MQ=255) ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACg > 1:473691/1‑70 (MQ=255) ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACg > 1:2059866/1‑70 (MQ=255) tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCCTTTCGGTGGCGTTTTTTATCGCCCACGCa > 1:3598320/1‑71 (MQ=255) | GTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCA > minE/2065555‑2065693 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |