Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,562,553 A→G 10.6% C508C (TGT→TGC lldP ← L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,562,5530AG10.6% 228.9 / 14.7 85C508C (TGT→TGClldPL‑lactate permease
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base G (5/4);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCA  >  minE/2562483‑2562623
                                                                      |                                                                      
tGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGca                                                                        <  1:2385219/71‑1 (MQ=255)
  aaGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGcaca                                                                      <  1:3539052/71‑1 (MQ=255)
  aaGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGca                                                                        >  1:2012181/1‑69 (MQ=255)
   aGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAg                                                                     >  1:2659662/1‑71 (MQ=255)
       cAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcga                                                                 >  1:2126940/1‑71 (MQ=255)
       cAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcga                                                                 >  1:3612896/1‑71 (MQ=255)
              tgtTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgatagcga                                                           <  1:2106502/70‑1 (MQ=255)
              tgtTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgatagcg                                                            <  1:3718222/69‑1 (MQ=255)
              tgtTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgatagcg                                                            <  1:1385297/69‑1 (MQ=255)
                  tGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGa                                                       <  1:900050/70‑1 (MQ=255)
                     cAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGc                                                   <  1:2649008/71‑1 (MQ=255)
                     cAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGc                                                   >  1:2347841/1‑71 (MQ=255)
                      aGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCg                                                  >  1:563070/1‑71 (MQ=255)
                       gTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgg                                                 >  1:2822104/1‑71 (MQ=255)
                       gTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgg                                                 <  1:2653374/71‑1 (MQ=255)
                       gTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgg                                                 >  1:3663710/1‑71 (MQ=255)
                       gTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCCATAGATTGCgg                                                 >  1:2226946/1‑71 (MQ=255)
                        tAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCggg                                                <  1:701858/71‑1 (MQ=255)
                          aaGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgatagcgata                                                         <  1:2055413/60‑1 (MQ=255)
                          aaGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgggg                                               <  1:3463994/70‑1 (MQ=255)
                          aaGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCggg                                                >  1:2317926/1‑69 (MQ=255)
                          aaGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGga                                              <  1:278935/71‑1 (MQ=255)
                           aGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGca                                                                        <  1:264072/44‑1 (MQ=255)
                           aGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                                <  1:321257/52‑1 (MQ=255)
                           aGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                                <  1:873430/52‑1 (MQ=255)
                            gCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCggg                                                >  1:206254/1‑67 (MQ=255)
                             cGGAACAGTTCAGATTCTCGGCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGat                                           <  1:576580/71‑1 (MQ=255)
                             cGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgc                                                                   >  1:919084/1‑47 (MQ=255)
                              ggAACAGTTCAGATTCTCGGCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGat                                           >  1:138170/1‑70 (MQ=255)
                              ggAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGat                                           >  1:2284432/1‑70 (MQ=255)
                               gAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGatca                                         <  1:1998823/71‑1 (MQ=255)
                               gAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGatca                                         <  1:2599496/71‑1 (MQ=255)
                                aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                                <  1:1911568/47‑1 (MQ=255)
                                aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                                <  1:3593175/47‑1 (MQ=255)
                                aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                                <  1:2883883/47‑1 (MQ=255)
                                aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCACACAGGCGATAGCGATAGATTGCgggg                                               <  1:784215/64‑1 (MQ=255)
                                 aCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                                >  1:1191154/1‑46 (MQ=255)
                                 aCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGatcatc                                       >  1:3279114/1‑71 (MQ=255)
                                 aCAAATCAGACACTTTGCCCACCAGGCCTACCGCCGCACAGgcgatag                                                              >  1:1720285/1‑48 (MQ=255)
                                  cAAATCAGACTCTTTGCACACCAGGACTACCGCCGCACAGgcgatag                                                              >  1:15815/1‑47 (MQ=255)
                                    aaTCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCg                                                  <  1:2597273/57‑1 (MQ=255)
                                    aaTCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTa                                    <  1:389391/71‑1 (MQ=255)
                                     ttCAGATTCTCGGCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCttgc                                   <  1:809536/70‑3 (MQ=255)
                                     aTCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCtt                                     >  1:1068239/1‑69 (MQ=255)
                                      tCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgatagcgat                                                          <  1:245443/47‑1 (MQ=255)
                                      tCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTAcc                                  <  1:3769706/71‑1 (MQ=255)
                                       cAGACTCTTTGCCCACCAGGCCTACCGCCGCACAgg                                                                    >  1:1715191/1‑36 (MQ=255)
                                       cAGACTCTTTGCCCACCAGGCCTACCGCCGCACAgg                                                                    >  1:1395249/1‑36 (MQ=255)
                                       cAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgata                                                               >  1:1267607/1‑41 (MQ=255)
                                        aGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCg                                 <  1:1210235/70‑1 (MQ=255)
                                         gACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGt                               <  1:232027/71‑1 (MQ=255)
                                          aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTa                                    >  1:3142298/1‑65 (MQ=255)
                                          aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTa                                    >  1:2317194/1‑65 (MQ=255)
                                          aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                              <  1:2459157/71‑1 (MQ=255)
                                          aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                              <  1:2444038/71‑1 (MQ=255)
                                          aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                              <  1:2736074/71‑1 (MQ=255)
                                           ctctTTGCCCACCAGGCCTACCGCCGCACAGgcgatagcgata                                                         >  1:672374/1‑43 (MQ=255)
                                            tctTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgg                                                 <  1:676199/50‑1 (MQ=255)
                                              tttGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGatcat                                        >  1:3173815/1‑57 (MQ=255)
                                              tttGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGatcat                                        <  1:3885581/57‑1 (MQ=255)
                                              tttGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGc                          >  1:2717157/1‑71 (MQ=255)
                                              tttGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGc                          >  1:677748/1‑71 (MQ=255)
                                              tcggCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGc                          >  1:421311/4‑71 (MQ=255)
                                                tGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCa                        >  1:1028580/1‑71 (MQ=255)
                                                ggCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGccg                        >  1:3720287/2‑70 (MQ=255)
                                                 gCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGtta                             >  1:1182852/1‑63 (MQ=255)
                                                 gCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGcc                         <  1:679284/69‑1 (MQ=255)
                                                  cccACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCAcc                      <  1:2061269/71‑1 (MQ=255)
                                                   ccaccaTGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGccgccg                     >  1:159661/1‑71 (MQ=255)
                                                   ccaccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCg                     <  1:1039501/71‑1 (MQ=255)
                                                   ccaccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCg                     <  1:2547204/71‑1 (MQ=255)
                                                    caccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCgg                    <  1:732329/71‑1 (MQ=255)
                                                    caccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCgg                    <  1:659663/71‑1 (MQ=255)
                                                      ccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCtt                                     <  1:1725721/52‑1 (MQ=255)
                                                       caGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                 <  1:2312697/71‑1 (MQ=255)
                                                       caGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                 >  1:3293905/1‑71 (MQ=255)
                                                       caGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTCACCGGTGGCGTCACCggtaa                 <  1:2461919/71‑3 (MQ=38)
                                                        aGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCt                                      >  1:3529145/1‑49 (MQ=255)
                                                        aGGCCTACAGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCtt                                     >  1:1639258/1‑50 (MQ=255)
                                                         ggCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGgaga                                            <  1:2008618/42‑1 (MQ=255)
                                                              aCCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggc          >  1:1933971/1‑71 (MQ=255)
                                                               ccgccgCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcg         >  1:2456872/1‑71 (MQ=255)
                                                               ccgccgCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcg         >  1:464145/1‑71 (MQ=255)
                                                                cgcggcGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGccag                                <  1:2784748/47‑3 (MQ=16)
                                                                 gccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa               >  1:1912159/1‑63 (MQ=255)
                                                                 gccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcgg        >  1:3036914/1‑70 (MQ=255)
                                                                  ccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGt                               >  1:1236178/1‑46 (MQ=255)
                                                                  ccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcggc       >  1:1666617/1‑70 (MQ=255)
                                                                    gcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc    >  1:3799266/1‑71 (MQ=255)
                                                                     cacaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc    >  1:2603306/1‑70 (MQ=255)
                                                                     cacaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc    >  1:3749515/1‑70 (MQ=255)
                                                                      acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                 <  1:3764947/56‑1 (MQ=255)
                                                                      acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc    >  1:2817278/1‑69 (MQ=255)
                                                                      acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc    >  1:2893795/1‑69 (MQ=255)
                                                                      acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc    >  1:2224424/1‑69 (MQ=255)
                                                                      acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc    >  1:3236348/1‑69 (MQ=255)
                                                                      acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCa  <  1:2730025/71‑1 (MQ=255)
                                                                      acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCa  <  1:2734771/71‑1 (MQ=255)
                                                                      |                                                                      
TGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCA  >  minE/2562483‑2562623

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: