Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,439,594 G→A 73.4% G350G (GGC→GGT ilvC ← ketol‑acid reductoisomerase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,439,5940GA73.4% 97.1 / 35.6 90G350G (GGC→GGTilvCketol‑acid reductoisomerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (25/41);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCAGGTC  >  minE/2439541‑2439660
                                                     |                                                                  
cccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCa                                                     <  1:2303700/69‑1 (MQ=255)
cccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCa                                                     <  1:1714603/69‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:2132531/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:2349725/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:430240/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:972887/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:1973008/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:2908527/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:2884817/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:861423/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:1449830/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:751126/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:1507216/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:1556060/67‑1 (MQ=255)
 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:3372345/67‑1 (MQ=255)
 ccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCAt                                                    <  1:2205888/69‑1 (MQ=255)
 ccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCAt                                                    <  1:857752/69‑1 (MQ=255)
  tgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                    <  1:2055025/67‑1 (MQ=255)
            aTCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGcgcg                                         <  1:2823917/69‑1 (MQ=255)
            aTCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGcgcg                                         <  1:729987/69‑1 (MQ=255)
            aTCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGcgcg                                         <  1:1271245/69‑1 (MQ=255)
                    catcaGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCaaa                                 >  1:2022198/1‑69 (MQ=255)
                    catcaGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCaaa                                 >  1:3108301/1‑69 (MQ=255)
                    catcaGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCaaa                                 >  1:2738909/1‑69 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTcgcgct                                       >  1:1215614/1‑57 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:3363486/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:407626/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:2424372/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:2420361/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:2655557/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:2941171/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:2723211/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:2338262/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:2831171/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:1103788/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:1058972/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:1345210/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:1501191/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                      >  1:705299/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:446514/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:3101044/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:1332611/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:917150/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:2741129/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:837838/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:184368/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:2550862/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:2407302/1‑60 (MQ=255)
                        aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:448255/1‑60 (MQ=255)
                            cGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGttt                         >  1:2193150/1‑69 (MQ=255)
                              ccTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGtttt                        <  1:480427/68‑1 (MQ=255)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:609458/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:414092/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:843674/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:553616/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:847686/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:542605/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:528570/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:984077/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:3346002/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:3272026/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:2037832/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:1594587/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:159864/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:1740737/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:1576420/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:2028245/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:2099218/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:2118455/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:1321396/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:2625171/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:2732209/36‑1 (MQ=37)
                                 ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:3150607/36‑1 (MQ=37)
                                  tATCTAAGTACTCCTGTTCACGGATTTTGCCTTCa                                                     <  1:756715/35‑1 (MQ=16)
                                  tATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                     <  1:959802/35‑1 (MQ=37)
                                        aGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCtt             <  1:2040964/69‑1 (MQ=255)
                                        aGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCtt             <  1:673286/69‑1 (MQ=255)
                                        aGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCtt             <  1:3150380/69‑1 (MQ=255)
                                        aGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCAGTTTTGCCGGTCTCtt             <  1:3366375/69‑1 (MQ=255)
                                        aGTACTCCTGCTCGCAGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCtt             <  1:1381634/69‑1 (MQ=255)
                                         gTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTc            >  1:14385/1‑69 (MQ=255)
                                         gTACTCCTGCTCGCCGATTTTGCCTTCATACGGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTc            >  1:18252/1‑69 (MQ=255)
                                               cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGtt                                      >  1:2032181/1‑37 (MQ=255)
                                               cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCa      >  1:1251794/1‑69 (MQ=255)
                                               cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCa      >  1:1867630/1‑69 (MQ=255)
                                               cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCa      >  1:2361823/1‑69 (MQ=255)
                                               cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCAAGCCa      >  1:1761087/1‑69 (MQ=255)
                                                    cGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCAGGTc  <  1:1764161/68‑1 (MQ=255)
                                                    cGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCAGGTc  <  1:243076/68‑1 (MQ=255)
                                                    cGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCAGGTc  <  1:2928032/68‑1 (MQ=255)
                                                     |                                                                  
CCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCAGGTC  >  minE/2439541‑2439660

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: