| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 950,078 | C→T | 19.0% | Q253Q (CAG→CAA) | dos ← | cAMP phosphodiesterase, heme‑regulated |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 950,078 | 0 | C | T | 19.0% | 120.8 / 16.5 | 78 | Q253Q (CAG→CAA) | dos | cAMP phosphodiesterase, heme‑regulated |
| Reads supporting (aligned to +/- strand): ref base C (35/28); new base T (10/5); total (45/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.64e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01 | |||||||||||
GGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAA > minE/950028‑950144 | ggCGGGCTGCTGCACATGGCGGCGAGAATGTTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata < 1:1254843/69‑1 (MQ=255)ggCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata < 1:1853020/69‑1 (MQ=255)ggCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata < 1:3375929/69‑1 (MQ=255) ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat > 1:1091189/1‑66 (MQ=255) ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat > 1:3434758/1‑66 (MQ=255) ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat > 1:1622521/1‑66 (MQ=255) ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat > 1:3254591/1‑66 (MQ=255) ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat > 1:1594702/1‑66 (MQ=255) ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata > 1:1381104/1‑65 (MQ=255) tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa > 1:2984494/1‑68 (MQ=255) tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa > 1:1363648/1‑68 (MQ=255) tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa > 1:29316/1‑68 (MQ=255) tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa > 1:1469358/1‑68 (MQ=255) gcACATGGCGGCGTGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt < 1:3019002/69‑1 (MQ=255) gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt < 1:401562/69‑1 (MQ=255) gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt < 1:1589840/69‑1 (MQ=255) gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt < 1:3162779/69‑1 (MQ=255) gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt < 1:1381073/69‑1 (MQ=255) gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt < 1:2937796/69‑1 (MQ=255) gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt < 1:295681/69‑1 (MQ=255) ggcggcGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc < 1:3325624/69‑1 (MQ=255) ggcggcGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc < 1:2007924/69‑1 (MQ=255) gcggcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc < 1:2090764/68‑1 (MQ=255) gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg < 1:3231645/69‑1 (MQ=255) gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg < 1:417908/69‑1 (MQ=255) gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg < 1:2223975/69‑1 (MQ=255) cgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg < 1:305443/68‑1 (MQ=255) cgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg < 1:717137/68‑1 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTc > 1:1166038/1‑63 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAg > 1:2785020/1‑68 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAg > 1:2474796/1‑68 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:1883510/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:875451/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:2876978/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:705224/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:2563839/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:2110202/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:3427993/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:2260567/1‑69 (MQ=255) ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa > 1:2081042/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCg > 1:1783806/1‑42 (MQ=37) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGa > 1:2201782/1‑68 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:2325594/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:3042303/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:2122124/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:1679364/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:3263500/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:1283537/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:1407879/1‑69 (MQ=255) tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt > 1:2806684/1‑69 (MQ=255) cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc < 1:2089766/68‑1 (MQ=255) cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc < 1:2905556/68‑1 (MQ=255) cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc < 1:1419189/68‑1 (MQ=255) cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc < 1:2754517/68‑1 (MQ=255) cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc < 1:2328465/68‑1 (MQ=255) tCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaact < 1:1855073/44‑3 (MQ=255) cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa > 1:2574235/1‑60 (MQ=255) cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa > 1:1986410/1‑60 (MQ=255) cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa > 1:1812446/1‑60 (MQ=255) cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg > 1:970133/1‑68 (MQ=255) cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg > 1:1205740/1‑68 (MQ=255) cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg > 1:1230445/1‑68 (MQ=255) cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg > 1:2587904/1‑68 (MQ=255) cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc > 1:2770031/1‑69 (MQ=255) cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc > 1:467994/1‑69 (MQ=255) cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc > 1:1706527/1‑69 (MQ=255) cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc > 1:394926/1‑69 (MQ=255) cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc > 1:1964826/1‑69 (MQ=255) cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc > 1:2483638/1‑69 (MQ=255) attCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:338896/67‑1 (MQ=255) aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:2107135/69‑1 (MQ=255) aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:1861930/69‑1 (MQ=255) aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:1477346/69‑1 (MQ=255) aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:2639793/69‑1 (MQ=255) aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:2838490/69‑1 (MQ=255) aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:1130328/69‑1 (MQ=255) aTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATaa < 1:577141/69‑1 (MQ=255) aTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa < 1:2497709/68‑1 (MQ=255) | GGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAA > minE/950028‑950144 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |