Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,008,951 T→A 26.6% G95G (GGT→GGA folM → dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,008,9510TA26.6% 134.4 / 66.2 117G95G (GGT→GGAfolMdihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (45/41);  new base A (18/13);  total (63/54)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-04

Only 100 of 117 total aligned reads displayed.

CACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAAC  >  minE/1008891‑1009017
                                                            |                                                                  
cACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACt                                                            >  1:932150/1‑69 (MQ=255)
   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                                                         <  1:3336428/69‑3 (MQ=255)
   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                                                         <  1:1425/69‑3 (MQ=255)
   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                                                         <  1:1829886/69‑3 (MQ=255)
   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                                                         <  1:641244/69‑3 (MQ=255)
   ccATGGTCTGCGTGCTATTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTgac                                                         <  1:938481/69‑3 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGacgacg                                                    >  1:3369426/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGacgacg                                                    >  1:3365545/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:1502683/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:10038/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:1753047/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:1934117/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:219871/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:2233183/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:246253/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    <  1:2646691/69‑1 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:2827655/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:615102/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:511645/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:417861/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:454110/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:1318979/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:1068601/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:1351346/1‑69 (MQ=255)
        gTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    <  1:1161178/69‑1 (MQ=255)
        gTCTGCGTGCTATTTTGCATAAAGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    >  1:31869/1‑69 (MQ=255)
         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAACCCGGGTGCGCCACTGGCCGACGt                                                   >  1:2717526/1‑69 (MQ=255)
         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                                                   >  1:157002/1‑69 (MQ=255)
         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                                                   >  1:504548/1‑69 (MQ=255)
         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                                                   >  1:3406288/1‑69 (MQ=255)
         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                                                   >  1:2099654/1‑69 (MQ=255)
         tCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGt                                                   >  1:2097180/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)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:2976054/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:3154673/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1234763/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:343235/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:2838185/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1243061/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1060405/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1458919/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1301324/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1450130/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1809470/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1922855/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:2269492/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1326434/1‑69 (MQ=255)
                   aTTTTGCATAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGc                                         >  1:1327178/1‑69 (MQ=255)
                      ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat                                       <  1:486359/68‑1 (MQ=255)
                      ttGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCat                                       <  1:2411775/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)
                           tAACGCCAGCGCGTGGATGGCGGAAAAACCGGGAGCGCCACTGACCGACGTACTGGCTTGCATGATGc                                  <  1:1369277/68‑1 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCg                                                       >  1:2018385/1‑46 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                     >  1:2107127/1‑64 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                     >  1:895277/1‑64 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                     >  1:3377129/1‑64 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                     >  1:1371336/1‑64 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                     >  1:507561/1‑64 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                     >  1:2625120/1‑64 (MQ=255)
                            aaCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCatga                                     >  1:922923/1‑64 (MQ=255)
                             aCGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGa                               >  1:1942517/1‑69 (MQ=255)
                               gCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc                             <  1:709539/69‑1 (MQ=255)
                               gCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc                             <  1:463419/69‑1 (MQ=255)
                                     gCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                       <  1:1072164/69‑1 (MQ=255)
                                      cGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                       >  1:2281631/1‑68 (MQ=255)
                                      cGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                       >  1:2453051/1‑68 (MQ=255)
                                      cGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                       >  1:1630892/1‑68 (MQ=255)
                                      cGTGGACGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGtt                       >  1:710634/1‑68 (MQ=255)
                                       gTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACg                                                    <  1:54321/38‑1 (MQ=255)
                                        tGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                     <  1:3408865/68‑1 (MQ=255)
                                        tGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                     <  1:1651166/68‑1 (MQ=255)
                                        tGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                     <  1:2461879/68‑1 (MQ=255)
                                        tGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                     <  1:1253185/68‑1 (MQ=255)
                                        tGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                     <  1:460537/68‑1 (MQ=255)
                                        tGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                     <  1:933087/68‑1 (MQ=255)
                                        tGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTaa                     <  1:3243847/68‑1 (MQ=255)
                                         ggATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATa                   <  1:2952070/69‑1 (MQ=255)
                                          taTGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAc                  <  1:973623/68‑1 (MQ=255)
                                          gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAc                  <  1:419245/69‑1 (MQ=255)
                                          gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAc                  <  1:2950725/69‑1 (MQ=255)
                                          gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAc                  <  1:2799668/69‑1 (MQ=255)
                                          gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAc                  <  1:2795693/69‑1 (MQ=255)
                                          gATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAc                  <  1:2201045/69‑1 (MQ=255)
                                            tGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATAcc                 <  1:2363228/68‑1 (MQ=255)
                                                 gAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATa            <  1:511124/68‑1 (MQ=255)
                                                 gAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATa            <  1:1269850/68‑1 (MQ=255)
                                                  aaaaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATaac          <  1:707690/69‑3 (MQ=255)
                                                  aaaaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAcc          <  1:3217273/69‑1 (MQ=255)
                                                  aaaaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAcc          <  1:1143851/69‑1 (MQ=255)
                                                    aaaCCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATAcc          <  1:3312329/67‑1 (MQ=255)
                                                        cGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCa    <  1:3115554/69‑1 (MQ=255)
                                                          ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAAc  <  1:2467522/69‑1 (MQ=255)
                                                          ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAAc  <  1:2133853/69‑1 (MQ=255)
                                                          ggTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAAc  <  1:1854273/69‑1 (MQ=255)
                                                            |                                                                  
CACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATCCACGTTAATACCCCATACCTGCTCAAC  >  minE/1008891‑1009017

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: