| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 608,412 | A→G | 15.2% | G614G (GGT→GGC) | pflB ← | pyruvate formate lyase I |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 608,412 | 0 | A | G | 15.2% | 186.5 / 22.6 | 79 | G614G (GGT→GGC) | pflB | pyruvate formate lyase I |
| Reads supporting (aligned to +/- strand): ref base A (32/35); new base G (7/5); total (39/40) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.46e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01 | |||||||||||
GCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGT > minE/608342‑608481 | gCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa > 1:4342200/1‑71 (MQ=255) gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa < 1:1332168/71‑1 (MQ=255) aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa > 1:4628740/1‑71 (MQ=255) cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg < 1:3273572/71‑1 (MQ=255) cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg > 1:5220220/1‑71 (MQ=255) aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt < 1:1479695/71‑1 (MQ=255) cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt > 1:2916497/1‑70 (MQ=255) cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg < 1:4386550/71‑1 (MQ=255) cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg < 1:1196942/71‑1 (MQ=255) ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc > 1:321984/1‑53 (MQ=255) cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt > 1:3633108/1‑70 (MQ=255) cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg < 1:899285/71‑1 (MQ=255) aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg < 1:2862864/64‑1 (MQ=255) aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg < 1:5010659/64‑1 (MQ=255) aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt > 1:3351370/1‑71 (MQ=255) cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg > 1:3361145/1‑71 (MQ=255) gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt > 1:1669343/1‑70 (MQ=255) gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt > 1:2058523/1‑70 (MQ=255) gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt > 1:2981590/1‑70 (MQ=255) gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca < 1:3983295/49‑1 (MQ=255) gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca < 1:2340756/49‑1 (MQ=255) gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt < 1:210484/55‑1 (MQ=255) gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt < 1:4115035/55‑1 (MQ=255) cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg > 1:2490241/1‑70 (MQ=255) gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa > 1:781163/1‑71 (MQ=255) gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac < 1:1390336/44‑1 (MQ=255) gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg < 1:392881/57‑1 (MQ=255) gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa > 1:1605881/1‑71 (MQ=255) gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa > 1:2419006/1‑71 (MQ=255) gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa > 1:1599268/1‑71 (MQ=255) cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac < 1:3516012/43‑1 (MQ=255) cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCa > 1:2961006/1‑63 (MQ=255) cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGAc < 1:5056468/71‑1 (MQ=255) cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg < 1:985849/71‑1 (MQ=255) cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg < 1:1131484/71‑1 (MQ=255) gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTATaa > 1:964193/1‑50 (MQ=255) gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt > 1:779657/1‑69 (MQ=255) acgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg < 1:5523749/47‑1 (MQ=255) cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt > 1:1411661/1‑71 (MQ=255) gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTc < 1:2463970/71‑1 (MQ=255) acCGGCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg < 1:998088/71‑1 (MQ=255) ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTa > 1:2208333/1‑42 (MQ=255) ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg < 1:1327464/71‑1 (MQ=255) ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg > 1:3422674/1‑71 (MQ=255) cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg < 1:3416414/71‑1 (MQ=255) cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg < 1:3959433/71‑1 (MQ=255) cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg < 1:858281/71‑1 (MQ=255) gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa > 1:1857874/1‑58 (MQ=255) gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa > 1:5010318/1‑58 (MQ=255) tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa < 1:3858501/53‑1 (MQ=255) tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa < 1:3540171/53‑1 (MQ=255) ggggTGTTACCCGTTTTCTTACCATACACAACGTTAg < 1:1368028/37‑1 (MQ=255) ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGt > 1:5545369/1‑35 (MQ=12) ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa < 1:18036/37‑1 (MQ=12) ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGta > 1:2900815/1‑56 (MQ=255) ttACCCGTTTTCTTACCATACACAACGTTAGAAGt < 1:3801679/35‑1 (MQ=255) ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa < 1:1265460/46‑1 (MQ=255) ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg > 1:4274323/1‑71 (MQ=255) ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGCAt < 1:4533783/62‑1 (MQ=255) tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg < 1:4226166/58‑1 (MQ=255) tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt < 1:3955736/66‑1 (MQ=255) aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGAta > 1:4473590/1‑36 (MQ=12) aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa > 1:1242288/1‑35 (MQ=255) aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa > 1:3691654/1‑35 (MQ=255) aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg > 1:115970/1‑70 (MQ=255) gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg > 1:1062633/1‑66 (MQ=255) gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg > 1:3634530/1‑66 (MQ=255) gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt > 1:4710746/1‑71 (MQ=255) gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt > 1:2879597/1‑71 (MQ=255) gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt > 1:2541054/1‑71 (MQ=255) gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt > 1:2102100/1‑71 (MQ=255) ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa < 1:38625/41‑1 (MQ=255) ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa < 1:2378405/41‑1 (MQ=255) ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa < 1:5667625/41‑1 (MQ=255) ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc < 1:4317281/63‑1 (MQ=255) tttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACt > 1:4245046/1‑44 (MQ=12) ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt < 1:4831545/70‑1 (MQ=255) ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg > 1:1910693/1‑71 (MQ=255) ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg > 1:3934166/1‑55 (MQ=255) tCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa < 1:2128709/40‑1 (MQ=255) aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa < 1:3568046/68‑1 (MQ=255) aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt > 1:1920407/1‑70 (MQ=255) aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt < 1:4449503/70‑1 (MQ=255) | GCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGT > minE/608342‑608481 |
| 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: ‑ |