Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,562,550 C→G 9.8% A509A (GCG→GCC lldP ← L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,562,5500CG 9.8% 227.8 / 10.6 81A509A (GCG→GCClldPL‑lactate permease
Reads supporting (aligned to +/- strand):  ref base C (32/41);  new base G (5/3);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGC  >  minE/2562481‑2562618
                                                                     |                                                                    
ggTGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTAccgccg                                                                     <  1:3027271/71‑1 (MQ=255)
ggTGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTAccgccg                                                                     >  1:1728095/1‑71 (MQ=255)
 gTGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTAccgccg                                                                     >  1:2053432/1‑70 (MQ=255)
 gTGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTAccgccg                                                                     >  1:2058136/1‑70 (MQ=255)
 gTGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACcgccgc                                                                    <  1:3802873/71‑1 (MQ=255)
  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:1385297/69‑1 (MQ=255)
                tgtTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgatagcg                                                       <  1:3718222/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:3663710/1‑71 (MQ=255)
                         gTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgg                                            <  1:2653374/71‑1 (MQ=255)
                         gTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgg                                            >  1:2822104/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:2599496/71‑1 (MQ=255)
                                 gAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGatca                                    <  1:1998823/71‑1 (MQ=255)
                                  aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                           <  1:2883883/47‑1 (MQ=255)
                                  aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                           <  1:1911568/47‑1 (MQ=255)
                                  aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGgcgat                                                           <  1:3593175/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:1395249/1‑36 (MQ=255)
                                         cAGACTCTTTGCCCACCAGGCCTACCGCCGCACAgg                                                               >  1:1715191/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:2317194/1‑65 (MQ=255)
                                            aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTa                               >  1:3142298/1‑65 (MQ=255)
                                            aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                         <  1:2444038/71‑1 (MQ=255)
                                            aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                         <  1:2736074/71‑1 (MQ=255)
                                            aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                         <  1:2459157/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:659663/71‑1 (MQ=255)
                                                      caccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCgg               <  1:732329/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)
                                                                     |                                                                    
GGTGAAGATCAGGCTGTGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGC  >  minE/2562481‑2562618

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: