Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,562,703	C→G	9.5%	G458G (GGG→GGC)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,562,703	0	C	G	9.5%	178.7 / 7.9	63	G458G (GGG→GGC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base G (4/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCA  >  minE/2562634‑2562773                                                                      |                                                                       agagACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGAcc                                                                       >  1:37931/1‑71 (MQ=255) agagACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGAc                                                                        <  1:3658920/70‑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:3162486/65‑1 (MQ=255)               gttgtGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAgg                                                               <  1:1949096/65‑1 (MQ=255)                     ctgcGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAg                                                                >  1:1939526/1‑58 (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)                                    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:866392/1‑47 (MQ=255)                                              cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagcc                                                  >  1:988312/1‑47 (MQ=255)                                              cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                               >  1:3869225/1‑66 (MQ=255)                                              cAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                               >  1:3873898/1‑66 (MQ=255)                                                gggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGaa                                        >  1:3295031/1‑55 (MQ=255)                                                 ggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                       <  1:1327871/71‑1 (MQ=255)                                                 ggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGa                       <  1:124496/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:3308563/51‑1 (MQ=255)                                                           gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCgag                                 <  1:1882487/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)                                                                 tCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaaga                             >  1:1271012/1‑49 (MQ=255)                                                                 tCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATg                <  1:1101573/62‑1 (MQ=255)                                                                  ccGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt                 >  1:3584682/1‑60 (MQ=255)                                                                    gACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc    <  1:1735141/71‑1 (MQ=255)                                                                     aCCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc    >  1:154673/1‑70 (MQ=255)                                                                      cccGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa  <  1:2371450/71‑1 (MQ=255)                                                                      cccGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa  <  1:2318651/71‑1 (MQ=255)                                                                      |                                                                       AGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCA  >  minE/2562634‑2562773

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