Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,439,558 A→G 60.5% I362I (ATT→ATC 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,5580AG60.5% 42.0 / 57.9 86I362I (ATT→ATCilvCketol‑acid reductoisomerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (7/27);  new base G (35/17);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 2.47e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TCTTCAATGATGCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCG  >  minE/2439493‑2439621
                                                                 |                                                               
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACTCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:1684764/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:1111401/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:3077571/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:103095/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:2519116/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:2402000/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:744707/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:753982/1‑69 (MQ=255)
tctTCGATGATGCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGatc                                                              >  1:1262569/1‑69 (MQ=255)
tctTCAATGATGCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatc                                                              >  1:135680/1‑69 (MQ=255)
tctTCAATGATGCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatc                                                              >  1:3144318/1‑69 (MQ=255)
tctTCAATGATGCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAat                                                               >  1:3181218/1‑68 (MQ=255)
      atgatgCCGGAATCGAGCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:2956986/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:1423978/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:847354/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:1546740/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:1325765/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:3168586/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:1350197/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:752662/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:1773915/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:2767484/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:274342/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:751958/66‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAg                                                         <  1:3185525/68‑1 (MQ=255)
      atgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAg                                                         <  1:747094/68‑1 (MQ=255)
       tgatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTa                                                       <  1:1643091/69‑1 (MQ=255)
        gatgCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTa                                                       <  1:1973441/68‑1 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:1811497/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:1422337/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:822403/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:1557789/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:2987486/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:1321545/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:3256387/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:447391/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:567570/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:1149744/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:595857/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:604949/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:2037425/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:2761832/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:1882685/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:2582695/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:2571512/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:246220/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:950326/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:64500/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:662912/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:722120/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:21591/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:2143389/1‑69 (MQ=255)
        gatgCCGGAATCAACCATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:885345/1‑69 (MQ=255)
        gatgCCGGAATCAACAATGGTTTCAAACGCCAACTCAACGCCTGCTTTTACCATCGCGATCATCAGTAc                                                      >  1:3226208/1‑69 (MQ=255)
              ggAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAatcatc                                                           <  1:2855503/58‑1 (MQ=255)
                      ccATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAg                                         <  1:2931062/68‑1 (MQ=255)
                      ccATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAg                                         <  1:1268706/68‑1 (MQ=255)
                      ccATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAg                                         <  1:451198/68‑1 (MQ=255)
                      ccATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAg                                         <  1:63579/68‑1 (MQ=255)
                      ccATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAg                                         <  1:2829540/68‑1 (MQ=255)
                           gTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTcc                                   >  1:3405036/1‑69 (MQ=255)
                           gTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTcc                                   >  1:354062/1‑69 (MQ=255)
                             ttCGACCGCCAGTTCAACGCCCGCTTTCACCATCGCAATCCTCAGTACGCCTTTATCGAAGTACTCCTg                                 >  1:3101100/1‑69 (MQ=255)
                              tCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATAGTACGCCTTTATCGAAGTACTCCTGCt                               >  1:454727/1‑69 (MQ=255)
                                   cGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTa                                                       <  1:531034/41‑1 (MQ=255)
                                                cccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCa              <  1:2303700/69‑1 (MQ=255)
                                                cccGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCa              <  1:1714603/69‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:1973008/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:972887/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:3372345/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:861423/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:1449830/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:1507216/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:1556060/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:430240/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:751126/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:2132531/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:2908527/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:2349725/67‑1 (MQ=255)
                                                 ctgCTTTTACCATCGCGATCATCAGTACGCCTTTATCGAAGTACTCCTGTTCACCGATTTTGCCTTCAt             <  1:2884817/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)
                                                                 |                                                               
TCTTCAATGATGCCGGAATCGACCATGGTTTCGAACGCCAGTTCAACGCCCGCTTTCACCATCGCAATCATCAGTACGCCTTTATCGAAGTACTCCTGCTCGCCGATTTTGCCTTCATACTGCGGCGCG  >  minE/2439493‑2439621

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: