Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,008,927 T→C 34.7% S87S (AGT→AGC folM → dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,008,9270TC34.7% 63.1 / 53.1 92S87S (AGT→AGCfolMdihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (43/17);  new base C (19/13);  total (62/30)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

GTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATC  >  minE/1008859‑1008990
                                                                    |                                                               
gtgATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGt                                                                 >  1:3227910/1‑69 (MQ=255)
gtgATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGt                                                                 >  1:2048267/1‑69 (MQ=255)
gtgATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGt                                                                 >  1:1237848/1‑69 (MQ=255)
gtgATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGt                                                                 >  1:2417255/1‑69 (MQ=255)
gtgATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGt                                                                 >  1:1292450/1‑69 (MQ=255)
gtgATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGt                                                                 >  1:1319137/1‑69 (MQ=255)
     ggCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTg                                                            <  1:3116149/69‑1 (MQ=255)
     ggCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTg                                                            <  1:2746461/69‑1 (MQ=255)
           tGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATgg                                                       >  1:700896/1‑68 (MQ=255)
                     gTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAcc                                             <  1:1527269/68‑1 (MQ=255)
                     gTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAcc                                             <  1:2523131/68‑1 (MQ=255)
                     gTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAcc                                             <  1:1571248/68‑1 (MQ=255)
                      tACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTg                                                            <  1:3276392/52‑1 (MQ=255)
                      tACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCgg                                           >  1:3333870/1‑69 (MQ=255)
                      tACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCg                                            >  1:2007075/1‑68 (MQ=255)
                                cACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACt                                 >  1:932150/1‑69 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:245597/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:1829886/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:2093638/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:1422577/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:938481/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:641244/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:3336428/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:1425/69‑3 (MQ=255)
                                   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                              <  1:559519/69‑3 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGacgacg                         >  1:3365545/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGacgacg                         >  1:3369426/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:246253/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:615102/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:219871/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:2233183/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:511645/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:454110/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:417861/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         <  1:2646691/69‑1 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:2827655/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:3316496/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:10038/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1318979/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1934117/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1351346/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1307946/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1068601/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1502683/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1753047/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         <  1:1161178/69‑1 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:1610763/1‑69 (MQ=255)
                                        gTCTGCGTGCTATTTTGCATAAAGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                         >  1:31869/1‑69 (MQ=255)
                                         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAACCCGGGTGCGCCACTGGCCGACGt                        >  1:2717526/1‑69 (MQ=255)
                                         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                        >  1:3406288/1‑69 (MQ=255)
                                         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                        >  1:504548/1‑69 (MQ=255)
                                         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                        >  1:157002/1‑69 (MQ=255)
                                         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                        >  1:2097180/1‑69 (MQ=255)
                                         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                        >  1:2099654/1‑69 (MQ=255)
                                          cTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTa                       <  1:1003554/69‑1 (MQ=255)
                                                   aTTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTg               <  1:1548996/68‑1 (MQ=255)
                                                   aTTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTg               <  1:2879330/68‑1 (MQ=255)
                                                   aTTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTg               <  1:3147421/68‑1 (MQ=255)
                                                   aTTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTg               <  1:3256018/68‑1 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAAcc                                             >  1:2495191/1‑38 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTg               >  1:1159466/1‑68 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1234763/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1326434/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1327178/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:343235/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1301324/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1450130/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1060405/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:3158689/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:3154673/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1458919/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:2976054/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1922855/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:2838185/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:2269492/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1243061/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc              >  1:1809470/1‑69 (MQ=255)
                                                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCAACTGACCGACGTACTGGCTTGc              >  1:976639/1‑69 (MQ=255)
                                                      ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat            <  1:2411775/68‑1 (MQ=255)
                                                      ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat            <  1:486359/68‑1 (MQ=255)
                                                           tAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGCATGATGc       <  1:1369277/68‑1 (MQ=255)
                                                           tAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGCATGATGc       <  1:1927501/68‑1 (MQ=255)
                                                           tAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGCATGATGc       <  1:1630897/68‑1 (MQ=255)
                                                           tAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGCATGATGc       <  1:3305723/68‑1 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCg                            >  1:2018385/1‑46 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:2946697/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:507561/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:2691728/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:2625120/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:1371336/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:3377129/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:895277/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:899144/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:922923/1‑64 (MQ=255)
                                                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga          >  1:2107127/1‑64 (MQ=255)
                                                             aCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGa    >  1:1942517/1‑69 (MQ=255)
                                                               gCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc  <  1:463419/69‑1 (MQ=255)
                                                               gCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc  <  1:709539/69‑1 (MQ=255)
                                                                    |                                                               
GTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATC  >  minE/1008859‑1008990

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: