Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,562,562	A→G	11.1%	A505A (GCT→GCC)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,562,562	0	A	G	11.1%	218.9 / 15.7	81	A505A (GCT→GCC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base A (40/32);  new base G (5/4);  total (45/36)
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

TGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGA  >  minE/2562497‑2562631                                                                  |                                                                      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:2653374/71‑1 (MQ=255)          gTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCgg                                                         >  1:2822104/1‑71 (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)               gCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCggg                                                        >  1:206254/1‑67 (MQ=255)                cGGAACAGTTCAGATTCTCGGCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGat                                                   <  1:576580/71‑1 (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)                   aaCAAATCAGACTCTTTGCCCACCAGGCCTACCGCCACACAGGCGATAGCGATAGATTGCgggg                                                       <  1:784215/64‑1 (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)                          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:2459157/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:677748/1‑71 (MQ=255)                                 tttGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGc                                  >  1:2717157/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:2893795/1‑69 (MQ=255)                                                         acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc            >  1:3236348/1‑69 (MQ=255)                                                         acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc            >  1:2817278/1‑69 (MQ=255)                                                         acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc            >  1:2224424/1‑69 (MQ=255)                                                         acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCa          <  1:2734771/71‑1 (MQ=255)                                                         acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCa          <  1:2730025/71‑1 (MQ=255)                                                           aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGAc                                    >  1:787385/1‑43 (MQ=255)                                                           aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGAc                                    >  1:3203072/1‑43 (MQ=255)                                                           aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                       >  1:620107/1‑56 (MQ=255)                                                           aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                       >  1:27069/1‑56 (MQ=255)                                                           aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAg         >  1:3752057/1‑70 (MQ=255)                                                           aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc        <  1:123004/71‑1 (MQ=255)                                                            ggCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAg         >  1:1264295/1‑69 (MQ=255)                                                            ggCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc        >  1:1758728/1‑70 (MQ=255)                                                             gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa      <  1:2127498/71‑1 (MQ=255)                                                              cgatagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa      >  1:3319831/1‑70 (MQ=255)                                                              cgatagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa      >  1:3345587/1‑70 (MQ=255)                                                               gatagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcggc               <  1:2875224/60‑1 (MQ=255)                                                                 tagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGa  <  1:3778737/71‑1 (MQ=255)                                                                  |                                                                      TGTTTGACAGTAAAGCGGAACAAATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGA  >  minE/2562497‑2562631

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: ‑