Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,439,591 C→T 75.1% E351E (GAG→GAA 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,5910CT75.1% 102.1 / 31.0 88E351E (GAG→GAAilvCketol‑acid reductoisomerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base T (25/41);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCA  >  minE/2439523‑2439656
                                                                    |                                                                 
tCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATAGTACGCCTTTATCGAAGTACTCCTGCt                                                                  >  1:454727/1‑69 (MQ=255)
                  cccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCa                                                 <  1:2303700/69‑1 (MQ=255)
                  cccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCa                                                 <  1:1714603/69‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:972887/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:2349725/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:2132531/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:1973008/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:2884817/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:2908527/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:3372345/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:1556060/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:1507216/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:1449830/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:861423/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:751126/67‑1 (MQ=255)
                   ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt                                                <  1:430240/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:729987/69‑1 (MQ=255)
                              aTCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGcgcg                                     <  1:2823917/69‑1 (MQ=255)
                              aTCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGcgcg                                     <  1:1271245/69‑1 (MQ=255)
                                      catcaGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCaaa                             >  1:2022198/1‑69 (MQ=255)
                                      catcaGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCaaa                             >  1:2738909/1‑69 (MQ=255)
                                      catcaGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCaaa                             >  1:3108301/1‑69 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTcgcgct                                   >  1:1215614/1‑57 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:2831171/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:2941171/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:2723211/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:2655557/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:3363486/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:2424372/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:2420361/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:705299/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:407626/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:1103788/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:1501191/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:2338262/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:1058972/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGTCGCGGtt                                  >  1:1345210/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:2407302/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:917150/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:2550862/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:446514/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:448255/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:3101044/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:837838/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:1332611/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:2741129/1‑60 (MQ=255)
                                          aGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCATACTGCGGCGCGGtt                                  >  1:184368/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:553616/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:542605/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:843674/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:528570/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:847686/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:1321396/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:984077/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:414092/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:1576420/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:2037832/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:159864/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:3346002/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:3272026/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:3150607/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:2118455/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:2732209/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:1740737/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:2099218/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:2625171/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:2028245/36‑1 (MQ=37)
                                                   ttATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:1594587/36‑1 (MQ=37)
                                                    tATCTAAGTACTCCTGTTCACGGATTTTGCCTTCa                                                 <  1:756715/35‑1 (MQ=16)
                                                    tATCGAAGTACTCCTGTTCACCGATTTTGCCTTCa                                                 <  1:959802/35‑1 (MQ=37)
                                                          aGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCtt         <  1:673286/69‑1 (MQ=255)
                                                          aGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCtt         <  1:2040964/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:1867630/1‑69 (MQ=255)
                                                                 cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCa  >  1:1251794/1‑69 (MQ=255)
                                                                 cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCa  >  1:2361823/1‑69 (MQ=255)
                                                                 cTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCAAGCCa  >  1:1761087/1‑69 (MQ=255)
                                                                    |                                                                 
TCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCGGTTTCAAACGCGGTTTTGCCGGTCTCTTCACGCCA  >  minE/2439523‑2439656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: