Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,562,691	T→G	12.5%	S462S (TCA→TCC)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,562,691	0	T	G	12.5%	134.0 / 18.0	64	S462S (TCA→TCC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base G (5/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

CCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGA  >  minE/2562621‑2562760                                                                       |                                                                      ccAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGAt                                                                       <  1:530358/71‑1 (MQ=255) ccAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGAt                                                                       >  1:2596171/1‑71 (MQ=255) ccAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGAt                                                                       >  1:290210/1‑71 (MQ=255)  cAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATg                                                                      <  1:3571440/71‑1 (MQ=255)  cAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATg                                                                      <  1:1153818/71‑1 (MQ=255)       aCAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTa                                                                 <  1:2141366/71‑1 (MQ=255)       aCAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTa                                                                 >  1:2191229/1‑71 (MQ=255)        cAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTAt                                                                >  1:2017512/1‑71 (MQ=255)        cAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTAt                                                                >  1:3479693/1‑71 (MQ=255)        cAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTAt                                                                <  1:3602733/71‑1 (MQ=255)          tgtCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTc                                                               >  1:3068789/3‑70 (MQ=255)          gATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATc                                                               >  1:3837183/1‑70 (MQ=255)              agagACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGagc                                                          <  1:1432784/71‑3 (MQ=255)              agagACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGagc                                                          <  1:3620180/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)                  aCGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGAt                                                                       <  1:2254696/54‑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)                        aTTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAg                                                                    <  1:3579992/51‑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: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: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)                                                                      aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa  >  1:1929560/1‑71 (MQ=255)                                                                      aTGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGa  >  1:3147750/1‑71 (MQ=255)                                                                       tGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGa                            >  1:2189709/1‑44 (MQ=255)                                                                       |                                                                      CCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCGGTGGCTTGCAGCGCGGCGAACAGGGCGTTAGATGAGGTATCCGACCCGGTCAGGAATACCCCCAGCCAGCCGAGGAACGGCGAGAAGAAGGTGAATGCATGA  >  minE/2562621‑2562760

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