| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 868,851 | A→C | 7.1% | K189T (AAA→ACA) | sohB → | predicted inner membrane peptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 868,851 | 0 | A | C | 7.1% | 145.1 / 6.2 | 56 | K189T (AAA→ACA) | sohB | predicted inner membrane peptidase |
| Reads supporting (aligned to +/- strand): ref base A (28/24); new base C (2/2); total (30/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01 | |||||||||||
CCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTT > minE/868784‑868915 | ccTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACaaaa < 1:773491/70‑1 (MQ=255) cTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAAtt > 1:944664/1‑71 (MQ=255) aaCTGTTACGGTAGCCAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGttt > 1:1549437/1‑71 (MQ=255) cTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTcc < 1:1261991/71‑1 (MQ=255) ttACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCg < 1:1173116/69‑1 (MQ=255) aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACaa > 1:3440219/1‑56 (MQ=255) aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAAt > 1:644012/1‑59 (MQ=255) aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAAt > 1:2420967/1‑59 (MQ=255) aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTATGTGGCGGACaa > 1:1075562/1‑56 (MQ=255) cGGTAGACAAAGTCGCTGCCAGCGGCGGTTATATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCg > 1:1018428/1‑71 (MQ=255) cGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCg > 1:149972/1‑71 (MQ=255) gTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGtt < 1:3323325/71‑1 (MQ=255) tAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGttt < 1:2924684/71‑1 (MQ=255) tAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGttt > 1:2281860/1‑71 (MQ=255) tAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGttt > 1:2203802/1‑71 (MQ=255) aGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtg > 1:3525435/1‑71 (MQ=255) ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTc > 1:3316596/1‑70 (MQ=255) ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTc > 1:492027/1‑70 (MQ=255) ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTc > 1:928853/1‑70 (MQ=255) ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGAAACAATTTTGTCCg < 1:1728659/48‑1 (MQ=25) ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGAAACAATTTTGTCCg < 1:3508469/48‑1 (MQ=25) aGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCAt < 1:2645047/71‑1 (MQ=255) aGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCAt < 1:562443/71‑1 (MQ=255) gcggcgGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtg < 1:2329611/60‑1 (MQ=255) gcggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGtt < 1:2093305/63‑1 (MQ=255) cggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCg < 1:2392541/71‑1 (MQ=255) ttACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtgg > 1:2977573/1‑71 (MQ=255) tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCCCCGTTTGCTATtgtg > 1:26578/1‑54 (MQ=255) tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg > 1:2386917/1‑54 (MQ=255) tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg > 1:2110890/1‑54 (MQ=255) tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg > 1:1827127/1‑54 (MQ=255) tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg > 1:1294628/1‑54 (MQ=255) tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTgg > 1:3551744/1‑65 (MQ=255) aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtggtg > 1:2495669/1‑71 (MQ=255) aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtggtg > 1:1736901/1‑71 (MQ=255) aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtggt < 1:2234813/70‑1 (MQ=255) tgGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa < 1:1073524/71‑1 (MQ=255) tgGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa < 1:2927328/71‑1 (MQ=255) tgGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCa > 1:1263836/1‑70 (MQ=255) ggCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaaa > 1:1160323/1‑71 (MQ=255) ggCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa < 1:1078019/70‑1 (MQ=255) tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGc > 1:328421/1‑70 (MQ=255) tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGc > 1:3550952/1‑70 (MQ=255) gtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTggg > 1:1267936/1‑43 (MQ=255) tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGccc > 1:488293/1‑70 (MQ=255) tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGccc < 1:753912/70‑1 (MQ=255) tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCa < 1:255261/71‑1 (MQ=255) gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg < 1:3481199/49‑1 (MQ=255) gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg < 1:3000332/49‑1 (MQ=255) gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg < 1:2409160/49‑1 (MQ=255) gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg < 1:3130057/49‑1 (MQ=255) gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg < 1:3063692/49‑1 (MQ=255) gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg < 1:3059359/49‑1 (MQ=255) tgCGGAAACAATTTTGTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACt > 1:133263/2‑71 (MQ=38) tgCGGAAACAATTTTGTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACt > 1:2686665/2‑71 (MQ=38) gCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACtt < 1:971128/71‑1 (MQ=255) cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGgcgc < 1:955983/57‑1 (MQ=255) cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACttt < 1:2163775/71‑1 (MQ=255) gACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa > 1:1108831/1‑57 (MQ=255) aaaaTTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa < 1:1191850/54‑1 (MQ=255) aaaaTTGTTTCCGCACCGTTTGCTATTGTGGGGTCCATTGGCGTGGTGGCGCaa < 1:397670/54‑1 (MQ=255) aaaTTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa < 1:992382/53‑1 (MQ=255) | CCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTT > minE/868784‑868915 |
| 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: ‑ |