Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,562,697	A→G	13.6%	D460D (GAT→GAC)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,562,697	0	A	G	13.6%	164.3 / 16.8	66	D460D (GAT→GAC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (5/4);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGT  >  minE/2562627‑2562767                                                                       |                                                                       aCAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTa                                                                        >  1:2191229/1‑71 (MQ=255) aCAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTa                                                                        <  1:2141366/71‑1 (MQ=255)  cAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTAt                                                                       <  1:3602733/71‑1 (MQ=255)  cAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTAt                                                                       >  1:3479693/1‑71 (MQ=255)  cAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTAt                                                                       >  1:2017512/1‑71 (MQ=255)    tgtCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTc                                                                      >  1:3068789/3‑70 (MQ=255)    gATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATc                                                                      >  1:3837183/1‑70 (MQ=255)        agagACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGagc                                                                 <  1:3620180/71‑3 (MQ=255)        agagACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGagc                                                                 <  1:1432784/71‑3 (MQ=255)        agagACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGAcc                                                                 >  1:37931/1‑71 (MQ=255)        agagACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGAc                                                                  <  1:3658920/70‑1 (MQ=255)          agaCGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTa                                                                        <  1:54565/62‑1 (MQ=255)          agaCGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCg                                                               <  1:755822/71‑1 (MQ=255)            aCGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCgg                                                              <  1:2340970/70‑1 (MQ=255)            aCGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGt                                                             >  1:1128347/1‑71 (MQ=255)              gCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCa                                                           <  1:3304159/71‑1 (MQ=255)                cAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAg                                                          >  1:1955246/1‑70 (MQ=255)                  aTTTGTTGTGCTGCGGTGTCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                       <  1:3506186/71‑1 (MQ=255)                  aTTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                       >  1:961878/1‑71 (MQ=255)                   tttGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                       >  1:2686338/1‑70 (MQ=255)                   tttGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGaa                                                       >  1:684665/1‑70 (MQ=255)                    ttgttgTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAAt                                                      <  1:514334/70‑1 (MQ=255)                      gttgtGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAgg                                                         <  1:1949096/65‑1 (MQ=255)                      gttgtGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAgg                                                         <  1:3162486/65‑1 (MQ=255)                            ctgcGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAg                                                          >  1:1939526/1‑58 (MQ=255)                             tgcGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTa                                                                        >  1:2398235/1‑43 (MQ=255)                               cGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagcca                                           <  1:2200930/70‑1 (MQ=255)                                     cTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCAAGCCAGCCGAgg                                    <  1:1715006/71‑1 (MQ=255)                                      ttGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccag                                              <  1:1569462/60‑1 (MQ=255)                                      ttGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAgg                                    >  1:1927720/1‑70 (MQ=255)                                        gCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATAc                                                    >  1:1903400/1‑52 (MQ=255)                                        gCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccag                                              >  1:619563/1‑58 (MQ=255)                                        gCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccag                                              >  1:1131862/1‑58 (MQ=255)                                        gCAGCGCGGCGAACAGGGCGTTAGAAGAGGTATCCGa                                                                   >  1:784628/1‑37 (MQ=255)                                           gcgcGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGc                              <  1:1087162/71‑1 (MQ=255)                                            cgcgGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagccag                                          >  1:3417283/1‑58 (MQ=255)                                               ggcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATAccc                                                  >  1:3138321/1‑47 (MQ=255)                                                gcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                         <  1:1758913/71‑1 (MQ=255)                                                   aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACggt                              >  1:1096630/1‑62 (MQ=255)                                                   aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACgg                               >  1:2521122/1‑62 (MQ=255)                                                   aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGagaag                        <  1:83398/69‑1 (MQ=255)                                                     cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagcc                                            >  1:988312/1‑47 (MQ=255)                                                     cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagcc                                            >  1:866392/1‑47 (MQ=255)                                                     cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                         >  1:3873898/1‑66 (MQ=255)                                                     cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                         >  1:3869225/1‑66 (MQ=255)                                                       gggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGaa                                  >  1:3295031/1‑55 (MQ=255)                                                        ggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                 <  1:124496/71‑1 (MQ=255)                                                        ggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                 <  1:1327871/71‑1 (MQ=255)                                                          cGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAATGTGAAt               <  1:505031/71‑1 (MQ=255)                                                          cGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaa                >  1:3290278/1‑70 (MQ=255)                                                            ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                 <  1:1370431/67‑1 (MQ=255)                                                            ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAAt               <  1:2929700/69‑1 (MQ=255)                                                            ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGc             <  1:2770260/71‑1 (MQ=255)                                                             tAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaacat             >  1:1457708/1‑67 (MQ=255)                                                              aGATGAGGTATCCGACCCGGTCAGGAATACCCccagcca                                           >  1:606549/1‑39 (MQ=255)                                                              aGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCa            <  1:850830/70‑1 (MQ=255)                                                              aGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt           >  1:3828736/1‑71 (MQ=255)                                                               gATGAGGTATCCGACCCGGTCAGGAATACCCccagccag                                          <  1:3538920/39‑1 (MQ=255)                                                               gATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt           <  1:995191/70‑1 (MQ=255)                                                                aggAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACg                                <  1:1772796/46‑1 (MQ=19)                                                                aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa         >  1:1929560/1‑71 (MQ=255)                                                                aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa         >  1:3147750/1‑71 (MQ=255)                                                                 tGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGa                                   >  1:2189709/1‑44 (MQ=255)                                                                  gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCgag                           <  1:1882487/51‑1 (MQ=255)                                                                  gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCgag                           <  1:3308563/51‑1 (MQ=255)                                                                  gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGAcc       >  1:618523/1‑71 (MQ=255)                                                                   aGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCg      <  1:2615893/71‑1 (MQ=255)                                                                    ggTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAgg                    >  1:1576350/1‑56 (MQ=255)                                                                    ggTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCgg     <  1:1222223/71‑1 (MQ=255)                                                                      tATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGgtg   <  1:721815/71‑1 (MQ=255)                                                                       aTCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATg          <  1:2799648/63‑1 (MQ=255)                                                                       aTCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGgtgt  >  1:3640938/1‑71 (MQ=255)                                                                       |                                                                       ACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGT  >  minE/2562627‑2562767

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑