Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,562,601 A→G 5.7% G492G (GGT→GGC lldP ← L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,562,6010AG 5.7% 301.4 / 6.1 104G492G (GGT→GGClldPL‑lactate permease
Reads supporting (aligned to +/- strand):  ref base A (40/58);  new base G (4/2);  total (44/60)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 104 total aligned reads displayed.

GCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTG  >  minE/2562532‑2562666
                                                                     |                                                                 
gCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCAc                                                                  >  1:1941283/1‑71 (MQ=255)
gCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCAc                                                                  >  1:5111582/1‑71 (MQ=255)
 cccACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCAc                                                                  <  1:2776999/70‑1 (MQ=255)
  ccaccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCg                                                                >  1:3093636/1‑71 (MQ=255)
      caGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                                                            <  1:1457428/71‑1 (MQ=255)
       aGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                                                            >  1:3101389/1‑70 (MQ=255)
        ggCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                                                          >  1:4464874/1‑71 (MQ=255)
           cTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTg                                                       >  1:4143537/1‑71 (MQ=255)
             gtcgcggcGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGCCGCCGCTGGTGTTtgc                                                     >  1:2997310/3‑71 (MQ=255)
             gtcgcggcGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGCCGCCGCTGGTGTTtgc                                                     >  1:4290080/3‑71 (MQ=255)
              tcgcggcGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGCCGCCGCTGGTGTTtgct                                                    >  1:2850094/2‑70 (MQ=255)
              ccgccgCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggc                                                     <  1:1358432/70‑1 (MQ=255)
               cgccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTgg                                                      <  1:4372768/68‑1 (MQ=255)
                gccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcgg                                                   >  1:4962567/1‑70 (MQ=255)
                   gcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                               <  1:38867/71‑1 (MQ=255)
                   gcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                               <  1:905535/71‑1 (MQ=255)
                     acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                               >  1:4993910/1‑69 (MQ=255)
                      caGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                                                          >  1:3867121/1‑57 (MQ=255)
                      caGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                                                          >  1:1403883/1‑57 (MQ=255)
                       aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc                                           <  1:5458937/71‑1 (MQ=255)
                       aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc                                           <  1:4477575/71‑1 (MQ=255)
                        ggCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggt                                                              >  1:1532170/1‑51 (MQ=255)
                        ggCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc                                           >  1:1784924/1‑70 (MQ=255)
                         gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc                                           <  1:3603062/69‑1 (MQ=255)
                         gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         <  1:3190675/71‑1 (MQ=255)
                         gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         >  1:3198310/1‑71 (MQ=255)
                         gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         >  1:3692022/1‑71 (MQ=255)
                         gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         >  1:3698900/1‑71 (MQ=255)
                         gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         <  1:4956448/71‑1 (MQ=255)
                           gatagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         >  1:3006952/1‑69 (MQ=255)
                             tagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         >  1:760940/1‑67 (MQ=255)
                             tagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGa                                     <  1:4985427/71‑1 (MQ=255)
                             tagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGa                                     <  1:933193/71‑1 (MQ=255)
                             cggcgatagATTGCGGGGAGATCATCTTGCCAGTTACGCCGCCGCTGGTGTTTGCTGCCAc                                               >  1:3575025/3‑61 (MQ=255)
                              agcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtg                                                             >  1:4864617/1‑46 (MQ=255)
                              agcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtg                                                             >  1:884592/1‑46 (MQ=255)
                              agcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         >  1:1040689/1‑66 (MQ=255)
                              agcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                         >  1:2991359/1‑66 (MQ=255)
                               gcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGAt                                    <  1:3799621/70‑1 (MQ=255)
                               gcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGAt                                    <  1:3589432/70‑1 (MQ=255)
                                cgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATc                                   <  1:4468703/70‑1 (MQ=255)
                                cgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:3975860/71‑1 (MQ=255)
                                cgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:2206542/71‑1 (MQ=255)
                                cgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:2402859/71‑1 (MQ=255)
                                cgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:732248/71‑1 (MQ=255)
                                  atagATTGCGGGGAGATCATCTTGCCAGTTACGCCGCCgctggtgt                                                         <  1:5347731/46‑3 (MQ=255)
                                  atagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCaga                                <  1:4903605/71‑1 (MQ=255)
                                  atagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:2099888/69‑1 (MQ=255)
                                  atagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:4125777/69‑1 (MQ=255)
                                   tagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  >  1:3856092/1‑68 (MQ=255)
                                    agATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCagag                               >  1:5168142/1‑70 (MQ=255)
                                    agATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAgaga                              >  1:4826657/1‑71 (MQ=255)
                                    agATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAgaga                              >  1:3199942/1‑71 (MQ=255)
                                     gATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAgaga                              >  1:1835337/1‑70 (MQ=255)
                                     gATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAgaga                              >  1:4528684/1‑70 (MQ=255)
                                      aTTGCGGGGAGATCATCTTGCCAGTTACGCCGCCGCTGGTGTTTGCTGCCACCAGCa                                          <  1:4719050/57‑1 (MQ=255)
                                      aTTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:521037/65‑1 (MQ=255)
                                      aTTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:4274733/65‑1 (MQ=255)
                                      aTTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                                  <  1:466454/65‑1 (MQ=255)
                                      aTTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACg                            >  1:2290931/1‑71 (MQ=255)
                                         gCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGcc                          <  1:4432399/70‑1 (MQ=255)
                                             ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAAttt                     <  1:5347637/71‑1 (MQ=255)
                                             ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAAttt                     <  1:1706535/71‑1 (MQ=255)
                                             ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAAttt                     <  1:2742185/71‑1 (MQ=255)
                                             ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAAttt                     >  1:3305907/1‑71 (MQ=255)
                                             ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAAttt                     <  1:2518998/71‑1 (MQ=255)
                                             ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACAAGCAACAGATCAGAGACGCCAAttt                     <  1:4530105/71‑1 (MQ=255)
                                              gagaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttg                    <  1:725154/71‑1 (MQ=255)
                                              gagaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttg                    <  1:5653578/71‑1 (MQ=255)
                                              gagaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttg                    <  1:5386954/71‑1 (MQ=255)
                                               agaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgt                   <  1:1250462/71‑1 (MQ=255)
                                               agaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgt                   <  1:3876918/71‑1 (MQ=255)
                                                 atcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGcc                          >  1:1369123/1‑62 (MQ=255)
                                                 atcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGc                           <  1:5144430/61‑1 (MQ=255)
                                                  tcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACg                            <  1:2444955/59‑1 (MQ=255)
                                                  tcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACg                            <  1:1897474/59‑1 (MQ=255)
                                                  tcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACg                            <  1:1601814/59‑1 (MQ=255)
                                                   catcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg                 <  1:2369473/69‑1 (MQ=255)
                                                   catcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTtgttgt                >  1:5348102/1‑70 (MQ=255)
                                                   catcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTtgtg               >  1:1459403/1‑71 (MQ=255)
                                                     tcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAAttt                     <  1:5438080/63‑1 (MQ=255)
                                                      cTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg                 >  1:4879143/1‑66 (MQ=255)
                                                      cTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGtgc              >  1:2737973/1‑69 (MQ=255)
                                                       ttACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGtgctgc           >  1:411249/1‑71 (MQ=255)
                                                        tACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttg                    >  1:3090777/1‑61 (MQ=255)
                                                         aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGtgctgc           <  1:514004/69‑1 (MQ=255)
                                                         aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGtgctgc           <  1:2011742/69‑1 (MQ=255)
                                                         aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGtgct             <  1:2085545/67‑1 (MQ=255)
                                                         aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGtgct             <  1:2268682/67‑1 (MQ=255)
                                                         aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCgg         >  1:1921218/1‑71 (MQ=255)
                                                         aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCg          <  1:2193849/70‑1 (MQ=255)
                                                           cGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg                 <  1:2623471/61‑1 (MQ=255)
                                                           cGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg                 <  1:4991421/61‑1 (MQ=255)
                                                           cGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg                 <  1:4908609/61‑1 (MQ=255)
                                                           cGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg                 <  1:1003998/61‑1 (MQ=255)
                                                           cGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTg       <  1:4458571/71‑1 (MQ=255)
                                                             gTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCag                                 <  1:1497158/43‑1 (MQ=255)
                                                             gTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGc     <  1:1204056/71‑1 (MQ=255)
                                                                aCGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCtt   <  1:4462267/70‑1 (MQ=255)
                                                                  gCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTg  >  1:4052917/1‑69 (MQ=255)
                                                                     |                                                                 
GCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTG  >  minE/2562532‑2562666

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: