Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,112 A→G 42.6% S319S (TCT→TCC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,1120AG42.6% 27.7 / 76.3 68S319S (TCT→TCCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (24/15);  new base G (16/13);  total (40/28)
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

GGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCG  >  minE/1015044‑1015181
                                                                    |                                                                     
ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                                                      >  1:5459541/1‑70 (MQ=255)
ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                                       >  1:2699696/1‑69 (MQ=255)
ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                                       >  1:3118813/1‑69 (MQ=255)
ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTAGGTcagcag                                                                      >  1:2488653/1‑69 (MQ=255)
 gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                                    >  1:3866169/1‑71 (MQ=255)
 gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                                    >  1:4853656/1‑71 (MQ=255)
 gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                                    >  1:4189154/1‑71 (MQ=255)
 gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                                    <  1:4837800/71‑1 (MQ=255)
 gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTAGCTTTGATATTACGGTCagcaga                                                                     <  1:489632/70‑1 (MQ=255)
   ttATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagc                                                                   <  1:2001343/70‑1 (MQ=255)
   ttATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                                    >  1:943428/1‑69 (MQ=255)
    tgtGTTCCAGTTTTTCGATGCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagc                                                                   <  1:4692495/67‑1 (MQ=255)
     gtGTTCCAGTTTTTCGATCCATATACCTTCGCGGTTGATTTTCGCTTTAATGTTACCGTCagcggagcag                                                                 <  1:3984026/69‑1 (MQ=255)
     gtGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                                 <  1:2617848/69‑1 (MQ=255)
     gtGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                                 <  1:812238/69‑1 (MQ=255)
      tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                               <  1:4598355/71‑1 (MQ=255)
      tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                               <  1:5354398/71‑1 (MQ=255)
      tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                               <  1:1212523/71‑1 (MQ=255)
      tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                               <  1:3690239/71‑1 (MQ=255)
      tGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGga                                                               <  1:19175/71‑1 (MQ=255)
       gTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                               <  1:2179569/70‑1 (MQ=255)
       gTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGga                                                               <  1:1146097/70‑1 (MQ=255)
        ttCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                                                             >  1:1604107/1‑71 (MQ=255)
        ttCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                                                             >  1:3603061/1‑71 (MQ=255)
        ttCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                                                             <  1:1525592/71‑1 (MQ=255)
        ttCCAGTTTTTCGATCCAGATACCTTAGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                                                             <  1:4348262/71‑1 (MQ=255)
         tCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                                            >  1:3339198/1‑71 (MQ=255)
          ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                                            <  1:5284528/70‑1 (MQ=255)
          ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACg                                                           >  1:55464/1‑71 (MQ=255)
          ccAGTTTTTCGATCCAGATCCCCTGACGGGTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                                            <  1:99413/70‑1 (MQ=255)
          ccAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACg                                                           >  1:1805971/1‑71 (MQ=255)
           cAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                                                                 <  1:425225/64‑1 (MQ=255)
           cAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGGCAGCGGAGCAGGAGACGc                                                          <  1:4591599/71‑1 (MQ=255)
            aGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagca                                                                  >  1:2883009/1‑62 (MQ=255)
             gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                                 <  1:988283/62‑1 (MQ=255)
             gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                                 <  1:3605047/62‑1 (MQ=255)
             gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                                 <  1:1152191/62‑1 (MQ=255)
              tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                                                      <  1:1178992/56‑1 (MQ=255)
              tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                               <  1:3194806/63‑1 (MQ=255)
                tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagca                                                                  >  1:195876/1‑58 (MQ=255)
                tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                                                             >  1:1768927/1‑63 (MQ=255)
                tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCAtg                                                     >  1:5106137/1‑70 (MQ=255)
                   cGATCCAGATACCTTCGGGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCTGGAGACGCCCATgcc                                                   <  1:959164/70‑1 (MQ=255)
                   cGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgcc                                                   >  1:631225/1‑70 (MQ=255)
                    gATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccga                                                 <  1:2196846/71‑1 (MQ=255)
                       ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                                       >  1:1112260/1‑46 (MQ=255)
                       ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                                       >  1:1884346/1‑46 (MQ=255)
                       ccAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccgaccg                                              >  1:2889007/1‑71 (MQ=255)
                          gATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                                       <  1:4953286/43‑1 (MQ=255)
                             ccccTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg                                         >  1:2236504/1‑70 (MQ=255)
                             ccccTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCACATACCGACCGGGCAg                                         >  1:2294138/1‑70 (MQ=255)
                              cccTCACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGa                                       <  1:3037613/71‑1 (MQ=255)
                                cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGCCCGGGCAGGAt                                      >  1:5521129/1‑70 (MQ=255)
                                cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATg                                     >  1:3316979/1‑71 (MQ=255)
                                cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATg                                     >  1:4567470/1‑71 (MQ=255)
                                 tGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                                    >  1:2567914/1‑71 (MQ=255)
                                  gACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgag                                                              >  1:5123485/1‑44 (MQ=255)
                                  cgcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAgg                                                                >  1:2607614/3‑42 (MQ=25)
                                  cgcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAgg                                                                >  1:1368071/3‑42 (MQ=25)
                                   gcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgcgc                                  >  1:678811/2‑71 (MQ=255)
                                   gcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgc                                    <  1:3486527/68‑1 (MQ=255)
                                   aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGcc                                                         <  1:5279903/48‑1 (MQ=255)
                                     ggTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCg                                <  1:2419461/71‑1 (MQ=255)
                                      gTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                               <  1:1919808/71‑1 (MQ=255)
                                       ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCa                                                       >  1:2768811/1‑46 (MQ=255)
                                       ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCa                                                       >  1:492989/1‑46 (MQ=255)
                                       ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCa                                                       >  1:4669302/1‑46 (MQ=255)
                                       ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                               >  1:938789/1‑70 (MQ=255)
                                       ttGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                               >  1:4734037/1‑70 (MQ=255)
                                        tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                             >  1:4244530/1‑71 (MQ=255)
                                        tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                             >  1:1126322/1‑71 (MQ=255)
                                        tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                             >  1:2493085/1‑71 (MQ=255)
                                         gATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGAc                            >  1:5690256/1‑71 (MQ=255)
                                         gATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGAc                            >  1:5197910/1‑71 (MQ=255)
                                          aTCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgaccg                                              >  1:4283364/1‑52 (MQ=255)
                                           ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGa                                       <  1:1891375/56‑1 (MQ=255)
                                           ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGt                          <  1:519775/69‑1 (MQ=255)
                                               cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCa                                                       >  1:226671/1‑38 (MQ=255)
                                               cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGt                          <  1:2708153/67‑1 (MQ=255)
                                                gCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGAc                            >  1:1503710/1‑64 (MQ=255)
                                                gCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCa                      <  1:820874/70‑1 (MQ=255)
                                                     gATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGa                 <  1:5118496/70‑1 (MQ=255)
                                                      atatTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAAt               >  1:5249545/1‑71 (MQ=255)
                                                        atTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCa             >  1:4396330/1‑71 (MQ=255)
                                                         ttACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAc            >  1:4556881/1‑71 (MQ=255)
                                                           aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGAt                                      >  1:2811572/1‑43 (MQ=255)
                                                           aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                             >  1:146559/1‑52 (MQ=255)
                                                           aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                             >  1:3645640/1‑52 (MQ=255)
                                                           aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                             >  1:2761736/1‑52 (MQ=255)
                                                           aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgc          >  1:4825718/1‑71 (MQ=255)
                                                            cGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAcgcg         <  1:4992893/71‑1 (MQ=255)
                                                             ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgc          <  1:124909/69‑1 (MQ=255)
                                                             ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCgg        <  1:5372228/71‑1 (MQ=255)
                                                                 agcggagcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGt    >  1:5579161/1‑71 (MQ=255)
                                                                 agcagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATg     <  1:4931657/70‑1 (MQ=255)
                                                                   cagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgcg         <  1:746093/64‑1 (MQ=255)
                                                                   cagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCg  <  1:3972203/71‑1 (MQ=255)
                                                                    agagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCg  <  1:4041384/70‑1 (MQ=255)
                                                                    agagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCg  <  1:1080131/70‑1 (MQ=255)
                                                                    |                                                                     
GGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCG  >  minE/1015044‑1015181

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: