Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,350 G→A 27.1% intergenic (+118/‑1661) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3500GA27.1% 221.7 / 90.6 140intergenic (+118/‑1661)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/66);  new base A (20/18);  total (56/84)
Fisher's exact test for biased strand distribution p-value = 8.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 144 total aligned reads displayed.

GCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCG  >  minE/2571286‑2571420
                                                                   |                                                                      
 cATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                                  <  1:4267138/70‑1 (MQ=255)
   tatTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctc                                                               >  1:2722097/1‑71 (MQ=255)
   tatTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctc                                                               >  1:2352113/1‑71 (MQ=255)
    atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                                                              >  1:4021442/1‑71 (MQ=255)
    atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                                                              >  1:5669681/1‑71 (MQ=255)
         ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATg                                                                        <  1:5365132/56‑1 (MQ=255)
           cATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                                                       <  1:1926925/55‑1 (MQ=255)
            aTTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAt                                                      <  1:327776/71‑1 (MQ=255)
             ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                                                     >  1:598891/1‑71 (MQ=255)
             ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                                                     >  1:4836689/1‑71 (MQ=255)
              tAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGTAtt                                                     <  1:225473/70‑1 (MQ=255)
               aGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAttt                                                    <  1:2451652/70‑1 (MQ=255)
                     gatgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                             <  1:1022779/71‑1 (MQ=255)
                      atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATc                                              >  1:3425292/1‑71 (MQ=255)
                      atgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCAtt                                                          <  1:1042987/57‑1 (MQ=255)
                        gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                             <  1:3671395/68‑1 (MQ=255)
                        gatTTTACTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                                          >  1:734041/1‑71 (MQ=255)
                         atTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCtctc                                                             <  1:2802838/51‑1 (MQ=255)
                          tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                                  <  1:4008930/45‑1 (MQ=255)
                          tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                                  <  1:720328/45‑1 (MQ=255)
                          tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCtt                                        <  1:119029/71‑1 (MQ=255)
                          ttttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt                                         >  1:1490904/1‑71 (MQ=255)
                           tttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt                                         >  1:1847902/1‑70 (MQ=255)
                            tttCTTTAATATTCATTATGTTAATAAGTTCTCGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGt                                      >  1:5527992/1‑71 (MQ=255)
                             ttccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCtt                                        <  1:1754179/67‑1 (MQ=255)
                             ttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa                                                >  1:2233766/1‑60 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:327799/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:529762/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:2989061/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:2353435/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:5635651/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:2033084/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                                       <  1:1522883/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                                       <  1:4254032/68‑1 (MQ=255)
                              tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                                       <  1:4309757/68‑1 (MQ=255)
                              tCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                                     <  1:5695919/70‑1 (MQ=255)
                                cTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                                     <  1:5031128/69‑1 (MQ=255)
                                cTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                                   >  1:4191709/1‑71 (MQ=255)
                                cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                                    >  1:1524999/1‑71 (MQ=255)
                                cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                                    >  1:4842888/1‑71 (MQ=255)
                                   tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGcc                                >  1:837254/1‑71 (MQ=255)
                                   tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGcc                                >  1:109743/1‑71 (MQ=255)
                                   tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTa                                                >  1:2221760/1‑56 (MQ=255)
                                   tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCt                                   >  1:3759709/1‑69 (MQ=255)
                                   tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                                 >  1:5528958/1‑71 (MQ=255)
                                   tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                                 >  1:4743825/1‑71 (MQ=255)
                                   tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                                 >  1:3672754/1‑71 (MQ=255)
                                   tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                                 >  1:2674080/1‑71 (MQ=255)
                                    aaTATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                                   <  1:2164482/67‑1 (MQ=255)
                                     atatTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCg                               <  1:1877765/70‑1 (MQ=255)
                                        ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATc                                              >  1:4038816/1‑52 (MQ=255)
                                        ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATc                                              >  1:5095556/1‑52 (MQ=255)
                                        ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                                   <  1:1611576/63‑1 (MQ=255)
                                        ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                                   <  1:1716614/63‑1 (MQ=255)
                                        ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAgt                           <  1:5305193/71‑2 (MQ=255)
                                        ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                            <  1:1175257/70‑1 (MQ=255)
                                          cATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                         >  1:2803109/1‑71 (MQ=255)
                                           agcaTGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                                    >  1:5207379/4‑60 (MQ=255)
                                           agcaTGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                                    >  1:1731674/4‑60 (MQ=255)
                                           aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                             >  1:482338/1‑50 (MQ=255)
                                           aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                         <  1:3930766/70‑1 (MQ=255)
                                           aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATg                        <  1:4985725/71‑1 (MQ=255)
                                           aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTATATCCGTCTTGTTCTGCCGCAGGAt                         <  1:5026670/70‑1 (MQ=255)
                                            ttATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa                                                <  1:4250179/46‑1 (MQ=255)
                                            ttATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                                     <  1:4167997/57‑1 (MQ=255)
                                            ttATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                                     <  1:127074/57‑1 (MQ=255)
                                               tGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc                    <  1:4784700/71‑1 (MQ=255)
                                               tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:380857/53‑1 (MQ=255)
                                               tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                       <  1:4295897/53‑1 (MQ=255)
                                                 ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                   <  1:596909/70‑1 (MQ=255)
                                                 ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                   <  1:3411049/70‑1 (MQ=255)
                                                  taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCt                                         >  1:479064/1‑47 (MQ=255)
                                                  taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                                   >  1:3661204/1‑53 (MQ=255)
                                                  taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTg                                  >  1:4344942/1‑54 (MQ=255)
                                                  taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                  >  1:4223154/1‑70 (MQ=255)
                                                    ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCg                                            <  1:4732441/42‑1 (MQ=255)
                                                    ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGt                                      <  1:4440026/48‑1 (MQ=255)
                                                    ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                >  1:2060905/1‑70 (MQ=255)
                                                      aaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                       >  1:246516/1‑61 (MQ=255)
                                                        gTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc                    <  1:153790/62‑1 (MQ=255)
                                                           cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc         <  1:5004618/70‑1 (MQ=255)
                                                           cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc         <  1:5097885/70‑1 (MQ=255)
                                                           cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc         <  1:3511341/70‑1 (MQ=255)
                                                             tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc         <  1:3932380/69‑1 (MQ=255)
                                                             tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc         <  1:774512/69‑1 (MQ=255)
                                                               ggATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                       >  1:1598337/1‑38 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                         >  1:1823475/1‑51 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                         >  1:2566562/1‑51 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc         <  1:5436525/67‑1 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag      <  1:3413843/70‑1 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga       <  1:4358569/69‑1 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga       <  1:2161337/69‑1 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga       <  1:4149641/69‑1 (MQ=255)
                                                                gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga       <  1:4143204/69‑1 (MQ=255)
                                                                  tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                            <  1:3844992/46‑1 (MQ=255)
                                                                  tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa              <  1:4685267/60‑1 (MQ=255)
                                                                  tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa              <  1:1392556/60‑1 (MQ=255)
                                                                  tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa              <  1:3179875/60‑1 (MQ=255)
                                                                   gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg  >  1:463623/1‑71 (MQ=255)
                                                                   gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg  >  1:1537116/1‑71 (MQ=255)
                                                                   |                                                                      
GCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCG  >  minE/2571286‑2571420

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: