| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,858,845 | G→A | 14.1% | S373S (TCG→TCA) | amiB → | N‑acetylmuramoyl‑l‑alanine amidase II |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,858,845 | 0 | G | A | 14.1% | 210.8 / 22.0 | 78 | S373S (TCG→TCA) | amiB | N‑acetylmuramoyl‑l‑alanine amidase II |
| Reads supporting (aligned to +/- strand): ref base G (31/36); new base A (9/2); total (40/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.81e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGTTGCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGTTTGCTGGC > minE/2858778‑2858913 | agtTGCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCg > 1:1155864/1‑71 (MQ=255)agtTGCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCg > 1:2682520/1‑71 (MQ=255)agtTGCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCg > 1:800713/1‑71 (MQ=255) gCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGata > 1:2598254/1‑71 (MQ=255) gCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGata > 1:2468821/1‑71 (MQ=255) gCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGata > 1:733120/1‑71 (MQ=255) gCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGata > 1:1497047/1‑71 (MQ=255) aCGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATcc > 1:2678530/1‑71 (MQ=255) cGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATcc > 1:1182850/1‑70 (MQ=255) cGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATcc > 1:2074898/1‑70 (MQ=255) cGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATcc > 1:1778450/1‑70 (MQ=255) cATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCa > 1:2021490/1‑70 (MQ=255) cATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCAt > 1:1213662/1‑71 (MQ=255) cATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCAt > 1:3647488/1‑71 (MQ=255) cATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCAAt > 1:3146528/1‑71 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:836616/70‑1 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:1320913/70‑1 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:3216830/70‑1 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:3240367/70‑1 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:3393629/70‑1 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:1916607/70‑1 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:3757760/70‑1 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg > 1:1359501/1‑70 (MQ=255) ggCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAg > 1:3746240/1‑70 (MQ=255) aaTACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:330998/65‑1 (MQ=255) aaTACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:792682/65‑1 (MQ=255) cATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGa < 1:1481491/71‑1 (MQ=255) cATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGa < 1:1321943/71‑1 (MQ=255) cATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGa < 1:1464628/71‑1 (MQ=255) aTAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGaa < 1:3671818/71‑1 (MQ=255) aTAAACGTCGACCAGAACACGCCACCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCTaa < 1:3300810/71‑1 (MQ=255) tAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTAc < 1:681690/62‑1 (MQ=255) tAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGACGGATATCCCATCAGTAc < 1:253035/62‑1 (MQ=255) aaCGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAAcc > 1:998486/1‑71 (MQ=255) aaCGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAAcc > 1:3473046/1‑71 (MQ=255) aaCGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAAcc > 1:188600/1‑71 (MQ=255) aaCGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAAcc > 1:2172149/1‑71 (MQ=255) aCCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGtttt < 1:801082/70‑1 (MQ=255) aCCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGtttt < 1:2796433/70‑1 (MQ=255) aCCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGtttt < 1:2918926/70‑1 (MQ=255) aCCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGtttt < 1:3518120/70‑1 (MQ=255) aCCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGtttt < 1:1990598/70‑1 (MQ=255) aGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTAt > 1:2336354/1‑69 (MQ=255) aaCACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:2705706/45‑1 (MQ=255) aaCACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAg > 1:3080718/1‑70 (MQ=255) aaCACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAg > 1:3787676/1‑70 (MQ=255) acacGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCa > 1:2475468/1‑43 (MQ=255) acacGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCa > 1:2412311/1‑43 (MQ=255) acacGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAg < 1:3241814/44‑1 (MQ=255) cacGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaa > 1:1945470/1‑71 (MQ=255) cacGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaa > 1:1896044/1‑71 (MQ=255) cacGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaa > 1:1286454/1‑71 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTAc < 1:1619132/45‑1 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:2720463/71‑1 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:3215482/71‑1 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:3571778/71‑1 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:1404342/71‑1 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:3591391/71‑1 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:912678/71‑1 (MQ=255) acGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaa > 1:1905619/1‑70 (MQ=255) cGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATc > 1:375721/1‑64 (MQ=255) cGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATc > 1:372637/1‑64 (MQ=255) cGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:4068856/70‑1 (MQ=255) cGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:2510784/70‑1 (MQ=255) cAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTAc < 1:1334424/41‑1 (MQ=255) ccTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGa < 1:1526178/45‑1 (MQ=255) ttGGCGTTCTGCGCTCGCCGGATATCCCATCAGTAc < 1:2821786/36‑1 (MQ=255) ttGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCg > 1:3795023/1‑69 (MQ=255) gCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCgaaga < 1:3753053/70‑1 (MQ=255) gCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCgaaga < 1:968794/70‑1 (MQ=255) cTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGttt > 1:2988173/1‑71 (MQ=255) cTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGttt > 1:2992894/1‑71 (MQ=255) tGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGcaac < 1:2530472/53‑1 (MQ=255) tGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGTTTg > 1:214385/1‑71 (MQ=255) cgcTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAAc > 1:2401836/1‑64 (MQ=255) cTCACCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGTTTGCTgg > 1:2986571/1‑71 (MQ=255) cGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGTTTGCTGgc > 1:580506/1‑70 (MQ=255) cGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGTTTGCTGgc > 1:3578619/1‑70 (MQ=255) | AGTTGCAACGCATTGGCGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTTATCAGCAACAACAGCGAAGAACGTTTGCTGGC > minE/2858778‑2858913 |
| 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: ‑ |