| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,716,860 | A→C | 25.8% | intergenic (+118/+80) | yphH → / ← glyA | predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,716,860 | 0 | A | C | 25.8% | 98.8 / 46.9 | 66 | intergenic (+118/+80) | yphH/glyA | predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (24/25); new base C (5/12); total (29/37) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.57e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
ACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAA > minE/1716793‑1716927 | aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt > 1:2554609/1‑71 (MQ=255) cGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtt > 1:3195851/1‑71 (MQ=255) cGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtt > 1:1456999/1‑71 (MQ=255) caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg < 1:2660934/70‑1 (MQ=11) aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg > 1:1727867/1‑71 (MQ=255) aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg < 1:2893220/71‑1 (MQ=255) agcttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg < 1:3602840/67‑1 (MQ=11) aaaTTGCCTGATGCGCTGCGCTTATCAGGCCTACAAGTTGAATTGCAATTTATTGAATTTGCACATTTTTg < 1:1597560/71‑1 (MQ=2) aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGt < 1:1290110/71‑1 (MQ=255) gcttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg > 1:3071000/3‑69 (MQ=11) aaTTGCCTGATGCGCTGCGCTTATCAGGCCTACAAGTTGAATTGCAATTTATTGAATTTGCACATTTTTGt < 1:3293701/71‑1 (MQ=0) aaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTa < 1:2306066/71‑1 (MQ=255) aaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTa < 1:2052954/71‑1 (MQ=255) cttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAg < 1:1647948/70‑1 (MQ=11) cttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAg < 1:1199022/70‑1 (MQ=11) cttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAg < 1:198412/70‑1 (MQ=11) ttGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAgg < 1:2345146/71‑1 (MQ=255) ttGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAgg < 1:1777636/71‑1 (MQ=11) gCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc > 1:533538/1‑71 (MQ=255) gCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc > 1:318986/1‑71 (MQ=255) cTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg < 1:2518958/71‑1 (MQ=255) cTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg > 1:1166864/1‑70 (MQ=255) tGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGa > 1:896482/1‑71 (MQ=11) gATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGa > 1:2193709/1‑70 (MQ=11) aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt > 1:1539611/1‑70 (MQ=255) aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt < 1:1708868/70‑1 (MQ=255) aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa < 1:3098876/71‑1 (MQ=255) aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa < 1:3632001/71‑1 (MQ=255) aTGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGAt < 1:2435653/70‑1 (MQ=11) tGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGAt > 1:83255/1‑69 (MQ=11) gcgcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg < 1:1786627/58‑1 (MQ=255) cgcTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAg < 1:920009/70‑1 (MQ=11) gcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc < 1:2621248/71‑1 (MQ=255) gcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc > 1:2081729/1‑71 (MQ=255) gcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc < 1:3736786/71‑1 (MQ=255) tACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg > 1:3572999/1‑67 (MQ=255) gatgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt < 1:3627719/68‑1 (MQ=11) atgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGtt < 1:76576/69‑1 (MQ=11) atgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGtt < 1:3319468/69‑1 (MQ=11) cGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt < 1:650953/70‑1 (MQ=255) cGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTc < 1:3132762/71‑1 (MQ=255) cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATaa > 1:1928570/1‑64 (MQ=2) cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATaa > 1:2901644/1‑64 (MQ=2) cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCgc < 1:152775/69‑2 (MQ=11) cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTc > 1:1089706/1‑71 (MQ=2) ttATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTcgcg > 1:3753857/1‑71 (MQ=2) aTCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCgct < 1:1792260/65‑3 (MQ=2) tCAGGCCTACACGGAG‑ATTGCAATATATTTAATTTGCAAAGTTTTGt > 1:108144/1‑47 (MQ=255) tCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTa > 1:3556636/1‑48 (MQ=255) aGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt > 1:2792880/1‑61 (MQ=11) ggAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt > 1:3395276/1‑62 (MQ=255) ggAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt > 1:2561440/1‑62 (MQ=255) aag‑aTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCg < 1:2694764/39‑1 (MQ=11) ag‑aTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac < 1:1855622/71‑1 (MQ=11) g‑aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcg > 1:2952213/1‑53 (MQ=255) caTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaaca < 1:1112107/70‑1 (MQ=11) ttGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaacaat > 1:345471/1‑70 (MQ=11) atatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca > 1:736985/1‑71 (MQ=255) atatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca > 1:2907564/1‑71 (MQ=255) tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca < 1:2241551/70‑1 (MQ=255) atTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAt < 1:3313617/70‑1 (MQ=255) ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAtt < 1:376076/70‑1 (MQ=255) gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAtt > 1:235304/1‑68 (MQ=255) gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa < 1:3779635/70‑1 (MQ=255) ttGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa > 1:2126778/1‑68 (MQ=255) cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCaaa < 1:2845940/70‑1 (MQ=255) aaaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTg < 1:588361/61‑1 (MQ=255) | ACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAA > minE/1716793‑1716927 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |