| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,256,643 | A→T | 5.3% | L132Q (CTG→CAG) | rpoC ← | RNA polymerase, beta prime subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,256,643 | 0 | A | T | 5.3% | 227.2 / 3.5 | 76 | L132Q (CTG→CAG) | rpoC | RNA polymerase, beta prime subunit |
| Reads supporting (aligned to +/- strand): ref base A (27/45); new base T (2/2); total (29/47) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 | |||||||||||
CCAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTCGGGGA > minE/2256575‑2256713 | ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAg > 1:636292/1‑70 (MQ=255)ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGc > 1:165971/1‑71 (MQ=255)ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGc > 1:750534/1‑71 (MQ=255)ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGc > 1:2018698/1‑71 (MQ=255) gTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCatat < 1:23661/70‑1 (MQ=255) gTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATc < 1:1309723/71‑1 (MQ=255) tCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCg < 1:1901823/71‑1 (MQ=255) cATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCg > 1:3527816/1‑70 (MQ=255) aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCa > 1:1123132/1‑58 (MQ=255) aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGag > 1:889063/1‑71 (MQ=255) aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGag > 1:3196960/1‑71 (MQ=255) aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGag > 1:2528396/1‑71 (MQ=255) tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc < 1:1273968/71‑1 (MQ=255) tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc < 1:3640092/71‑1 (MQ=255) tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc < 1:2346255/71‑1 (MQ=255) tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc < 1:2664361/71‑1 (MQ=255) ccGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagca < 1:837355/70‑1 (MQ=255) ccGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagca < 1:2416458/70‑1 (MQ=255) cGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagcagc < 1:2322040/71‑1 (MQ=255) gCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAgcagca < 1:3274822/71‑1 (MQ=255) aTAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGAt < 1:2156873/71‑1 (MQ=255) tAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATa < 1:687848/71‑1 (MQ=255) aaCCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATAc > 1:159683/1‑71 (MQ=255) aaCCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATAc > 1:3023021/1‑71 (MQ=255) aCCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACg > 1:2767317/1‑71 (MQ=255) aCCACATAGGATTCAAAGTACAGTACGCGCTCGATATGCCGCTGCgc < 1:1204921/47‑2 (MQ=21) aCCACATAGGATTCAAAGTACAGTACGCGCTCGATATGCCGCTGCgc < 1:1617076/47‑2 (MQ=21) ccACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg < 1:2382304/71‑1 (MQ=255) ccACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg < 1:1426142/71‑1 (MQ=255) acaTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa < 1:2567060/71‑1 (MQ=255) caTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGc < 1:607181/44‑1 (MQ=255) caTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa > 1:1939258/1‑70 (MQ=255) aTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATc > 1:3602075/1‑48 (MQ=255) aTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATa < 1:464415/64‑1 (MQ=255) aTTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa > 1:742142/1‑64 (MQ=255) aTTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa > 1:496180/1‑64 (MQ=255) ttCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg < 1:1502088/61‑1 (MQ=255) cAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGa < 1:3724320/70‑1 (MQ=255) cAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt > 1:735935/1‑71 (MQ=255) aaaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt < 1:657277/70‑1 (MQ=255) aaaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt < 1:1880169/70‑1 (MQ=255) aaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGAc < 1:741210/49‑1 (MQ=255) aaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAttt < 1:2599756/71‑1 (MQ=255) agtacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATa < 1:2295406/48‑1 (MQ=255) tacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGAc < 1:3396431/41‑1 (MQ=255) tacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAcc < 1:3590960/71‑1 (MQ=255) tacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAc > 1:83095/1‑70 (MQ=255) cgcgTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGTAAc > 1:1220244/1‑68 (MQ=255) gcgTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGa < 1:1139027/65‑1 (MQ=255) gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCg < 1:543964/38‑1 (MQ=255) gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCg < 1:2409202/38‑1 (MQ=255) gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCg < 1:1513901/38‑1 (MQ=255) gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAttt < 1:138096/58‑1 (MQ=255) ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt < 1:2313438/55‑1 (MQ=255) ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCATGAACCAGAtgt < 1:959721/71‑1 (MQ=255) ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAg < 1:1393936/67‑1 (MQ=255) ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAg < 1:304898/67‑1 (MQ=255) ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAg < 1:3736660/67‑1 (MQ=255) ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgt < 1:2965208/71‑1 (MQ=255) ctCGATATGCCGCTGCGCGATATCGAACAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgt > 1:3287004/2‑71 (MQ=21) ctCGATATGCCGCTGCGCGATATCGAACAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgt > 1:944790/2‑71 (MQ=21) cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgtg > 1:390952/1‑70 (MQ=255) cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgc > 1:2152800/1‑71 (MQ=255) cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgc > 1:896976/1‑71 (MQ=255) cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgc > 1:1595031/1‑71 (MQ=255) gATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg < 1:2996281/40‑1 (MQ=255) gATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgcg > 1:2285704/1‑71 (MQ=255) gATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgcg > 1:1532327/1‑71 (MQ=255) tatCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAt > 1:2201370/1‑64 (MQ=255) atCGCGCAGCGGCATATCGAGCAGCAGAGCGATACGGGACGGCAGCGATTTCAGGAAcc < 1:556636/59‑1 (MQ=255) atCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAcc < 1:2727854/59‑1 (MQ=255) atCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAcc < 1:337916/59‑1 (MQ=255) cgcgCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTc > 1:2990736/1‑71 (MQ=255) cgcgCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTc > 1:2186607/1‑71 (MQ=255) gcgcAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTc < 1:2587135/70‑1 (MQ=255) gcAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTCggg < 1:2367150/71‑1 (MQ=255) gcAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTCggg < 1:2371330/71‑1 (MQ=255) aGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTCGGGGa < 1:1938928/71‑1 (MQ=255) aGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTCGGGGa < 1:137341/71‑1 (MQ=255) | CCAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTCGGGGA > minE/2256575‑2256713 |
| 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: ‑ |