Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,562,718	C→G	21.3%	G453G (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,718	0	C	G	21.3%	152.5 / 35.3	75	G453G (GGG→GGC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base C (22/37);  new base G (11/5);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCC  >  minE/2562652‑2562784                                                                   |                                                                   tgctgcGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCcca                                                                <  1:2720917/71‑1 (MQ=255)     gcGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCccagc                                                              >  1:1995086/1‑69 (MQ=255)        gTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCagccagcc                                                         >  1:3022863/1‑71 (MQ=255)         tGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCagccagcc                                                         <  1:1658321/70‑1 (MQ=255)          ggCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCagccagcc                                                         >  1:3255476/1‑69 (MQ=255)               gCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGa                                                    <  1:915637/69‑1 (MQ=255)                cAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAAc                                                  <  1:4768465/70‑1 (MQ=255)                cAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCAGTCAGGAATACCCCCAGCCAGCCGAGGAAc                                                  <  1:2596432/70‑1 (MQ=255)                 aGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACgg                                                >  1:383742/1‑71 (MQ=255)                     cggcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCagccagcc                                                         <  1:324851/58‑1 (MQ=255)                     cggcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCg                                              <  1:5575591/69‑1 (MQ=255)                       gcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCcca                                                                >  1:3175738/1‑49 (MQ=255)                       gcgAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaa                                          <  1:2752940/71‑1 (MQ=255)                         gAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaaga                                        <  1:699560/71‑1 (MQ=255)                          aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAgaagaa                                       <  1:1830777/71‑1 (MQ=255)                          aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAgaagaa                                       <  1:1829411/71‑1 (MQ=255)                           aCAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGaagaag                                      >  1:3731184/1‑71 (MQ=255)                           aCAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGaagaag                                      >  1:170310/1‑71 (MQ=255)                             aGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGt                                    <  1:1533219/71‑1 (MQ=255)                             aGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGt                                    <  1:4194809/71‑1 (MQ=255)                             aGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACgg                                                >  1:4234613/1‑59 (MQ=255)                             aGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAgaagaa                                       >  1:1320443/1‑68 (MQ=255)                              gggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTg                                   <  1:3965777/71‑1 (MQ=255)                              gggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTg                                   <  1:260290/71‑1 (MQ=255)                              gggCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCACCAGCCAGCCGAGGAACGGCGagaaga                                        <  1:4815674/66‑1 (MQ=255)                                gCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaa                                 >  1:4260842/1‑71 (MQ=255)                                gCGTTAGAGGAGGTGTCCGAGCCGCTAAGCAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaa                                 >  1:3230097/1‑71 (MQ=255)                                   ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaag                                         <  1:5207988/60‑1 (MQ=255)                                   ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAAt                                <  1:5457942/69‑1 (MQ=255)                                   ttAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAAt                                <  1:3331589/69‑1 (MQ=255)                                   ttAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaacat                              <  1:4839846/71‑4 (MQ=255)                                     agagGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACggt                                               >  1:4210752/1‑51 (MQ=19)                                     aGATGAGGTATCCGACCCGGTCAGGAATACCCCcagccagc                                                          >  1:4126596/1‑41 (MQ=255)                                     aGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCAt                            >  1:3162993/1‑71 (MQ=255)                                        tGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa                          <  1:1856228/70‑1 (MQ=255)                                         gAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGAcc                        >  1:2169980/1‑71 (MQ=255)                                          aGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGt                                    <  1:3951219/58‑1 (MQ=255)                                           ggTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGagaa                                          >  1:5360147/1‑51 (MQ=19)                                               tCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGTGagaag                                         >  1:3513565/1‑48 (MQ=255)                                               tCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGagaag                                         >  1:140899/1‑48 (MQ=255)                                                 cGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaac                                >  1:1739630/1‑54 (MQ=255)                                                 cGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcg                >  1:4110773/1‑71 (MQ=255)                                                  gACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc               <  1:2298358/71‑1 (MQ=255)                                                  gACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc               <  1:2511757/71‑1 (MQ=255)                                                  gACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc               <  1:4141274/71‑1 (MQ=255)                                                   aCCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGcc              <  1:3007749/71‑1 (MQ=255)                                                    cccGGTCAGGAATACCCCCAGCCAGACGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa             <  1:4808242/71‑1 (MQ=255)                                                     ccGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAg            >  1:977483/1‑71 (MQ=255)                                                        gTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATg                               <  1:274018/49‑1 (MQ=255)                                                        gTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa             <  1:2022168/67‑1 (MQ=255)                                                        gTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa             <  1:2223517/67‑1 (MQ=255)                                                        gTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGc         <  1:813294/71‑1 (MQ=255)                                                          cAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGt           >  1:5265785/1‑67 (MQ=255)                                                          cAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGt           >  1:2502727/1‑67 (MQ=255)                                                          cAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGt           >  1:4962244/1‑67 (MQ=255)                                                          aaGGAATACGCCCAGCCAGCCGAGAAACGGTGagaaga                                        >  1:5273705/2‑38 (MQ=21)                                                           aGGAATACGCCCAGCCAGCCGAGAAACGGTGagaaga                                        >  1:3669812/1‑37 (MQ=25)                                                           aGGAATACGCCCAGCCAGCCGAGAAACGGTGagaaga                                        >  1:3549661/1‑37 (MQ=25)                                                           aGGAATACGCCCAGCCAGCCGAGAAACGGTGagaaga                                        >  1:4663760/1‑37 (MQ=25)                                                            ggAATACGCCCAGCCAGCCGAGAAACGGTGagaag                                         <  1:3911580/35‑1 (MQ=21)                                                            ggAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCAAGTGCCAgc     >  1:820881/1‑71 (MQ=255)                                                             gAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGcc              <  1:912513/61‑1 (MQ=255)                                                              aaTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa             <  1:1851296/61‑1 (MQ=255)                                                              aaTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGc         <  1:4579725/65‑1 (MQ=255)                                                              aaTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGc         <  1:4916532/65‑1 (MQ=255)                                                              aaTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGc         <  1:813052/65‑1 (MQ=255)                                                               aTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGAc                         >  1:5112987/1‑48 (MQ=255)                                                               aTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGAc                         >  1:1599010/1‑48 (MQ=255)                                                               aTACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGcc  <  1:306401/71‑1 (MQ=255)                                                                tACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGaa                                 <  1:802741/39‑1 (MQ=16)                                                                tACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTg                                   <  1:1308663/37‑1 (MQ=255)                                                                tACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc               <  1:578560/57‑1 (MQ=255)                                                                tACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTgcgc               <  1:923332/57‑1 (MQ=255)                                                                tACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcgc   >  1:5560770/1‑69 (MQ=255)                                                                tACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcgc   >  1:1001235/1‑69 (MQ=255)                                                                  cccccAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCa             >  1:5599825/1‑57 (MQ=255)                                                                   ccccAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAgcgc   <  1:1302933/66‑1 (MQ=255)                                                                   |                                                                   TGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCC  >  minE/2562652‑2562784

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