Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,202,140	C→A	7.9%	T431T (ACC→ACA)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,202,140	0	C	A	7.9%	290.1 / 6.7	88	T431T (ACC→ACA)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (41/40);  new base A (3/4);  total (44/44)
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.32e-01

GGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACT  >  minE/1202074‑1202210                                                                   |                                                                       ggCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGt                                                                     >  1:2007158/1‑70 (MQ=255)   ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                                  <  1:3473038/71‑1 (MQ=255)   ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                                  <  1:1463917/71‑1 (MQ=255)   ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                                  >  1:3828835/1‑71 (MQ=255)     ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGa                                                                <  1:712785/71‑1 (MQ=12)     ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                >  1:2208185/1‑71 (MQ=255)                cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                     <  1:2452681/71‑1 (MQ=255)                 aGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                     <  1:1034642/70‑1 (MQ=255)                  gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                             >  1:2898195/1‑61 (MQ=255)                  gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                             >  1:2380527/1‑61 (MQ=255)                  gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                   <  1:2126291/71‑1 (MQ=255)                  gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                                    >  1:489676/1‑70 (MQ=255)                   tGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCgg                                                                      <  1:3799768/51‑1 (MQ=255)                     cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                >  1:1360214/1‑71 (MQ=255)                     cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                >  1:2019384/1‑71 (MQ=255)                     cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                <  1:2214042/71‑1 (MQ=255)                     cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAAAAc                                                >  1:2187312/1‑71 (MQ=255)                       ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                 >  1:2645408/1‑68 (MQ=255)                       ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                               <  1:3132007/70‑1 (MQ=255)                       ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcg                                              >  1:1763813/1‑71 (MQ=255)                         gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                                    >  1:375368/1‑63 (MQ=255)                          cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                                      >  1:326770/1‑60 (MQ=255)                          cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                           <  1:2960531/71‑1 (MQ=255)                               ctcGCTGGATAAGGTCATCGAATCGATGTACGAAACCGGTAAGGa                                                                <  1:2101549/45‑1 (MQ=255)                               ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                <  1:3929155/45‑1 (MQ=255)                               ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACa                                                              >  1:2464782/1‑47 (MQ=255)                                 cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCa                                    >  1:3246103/1‑71 (MQ=255)                                  gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAc                                   <  1:2278689/71‑1 (MQ=255)                                   cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                  <  1:901792/71‑1 (MQ=255)                                   cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                  <  1:1667805/71‑1 (MQ=255)                                     ggATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                <  1:4065938/39‑1 (MQ=255)                                      gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa                                                          >  1:2266177/1‑44 (MQ=255)                                      gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                                <  1:190977/70‑1 (MQ=255)                                      gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                  <  1:3972827/68‑1 (MQ=255)                                        tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg                             <  1:382465/71‑3 (MQ=9)                                        tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg                             >  1:461654/1‑69 (MQ=9)                                         aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                            <  1:355144/39‑1 (MQ=255)                                         aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                                        >  1:1684957/1‑59 (MQ=255)                                         aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg                            >  1:2299572/1‑71 (MQ=255)                                          aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaa                                          <  1:1213629/56‑1 (MQ=255)                                          aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg                            <  1:969779/70‑1 (MQ=255)                                            gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                     <  1:4137919/43‑1 (MQ=255)                                             ttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                                     <  1:2694751/40‑1 (MQ=17)                                             ttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                                     <  1:310864/40‑1 (MQ=17)                                             tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTg                        <  1:2297868/71‑1 (MQ=255)                                               aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                         >  1:3159503/1‑68 (MQ=255)                                                tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                     <  1:3073080/39‑1 (MQ=255)                                                tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCt                                      >  1:3004455/1‑54 (MQ=255)                                                tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                 <  1:142465/59‑1 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                   >  1:298657/1‑39 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                                         >  1:2227986/1‑49 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                                         >  1:2223297/1‑49 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                                        >  1:2673881/1‑50 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                   >  1:3157823/1‑71 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                   >  1:4121969/1‑71 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                   >  1:3602236/1‑71 (MQ=255)                                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACAGCGAAAc                                        >  1:3177867/1‑50 (MQ=255)                                                  gAAACCATGTACGAAACAGGTAAAGATATGAACGcc                                                      >  1:1676812/1‑36 (MQ=12)                                                  gAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt                   >  1:743552/1‑71 (MQ=17)                                                    aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                 <  1:2023539/71‑1 (MQ=255)                                                     aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                                        >  1:2836474/1‑47 (MQ=255)                                                     aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                      >  1:562780/1‑65 (MQ=255)                                                     aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaa                <  1:1724198/71‑1 (MQ=255)                                                     aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                 >  1:1993506/1‑70 (MQ=255)                                                     aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                 >  1:3562572/1‑70 (MQ=255)                                                      cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                 >  1:2792915/1‑69 (MQ=255)                                                      cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                 >  1:4007118/1‑69 (MQ=255)                                                       gATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCt                                      >  1:3024897/1‑47 (MQ=255)                                                        aTGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTg                        <  1:294420/60‑1 (MQ=255)                                                         tGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc            <  1:1014421/71‑1 (MQ=255)                                                          gTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                               <  1:2176499/35‑1 (MQ=255)                                                          gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                    <  1:3760122/62‑1 (MQ=255)                                                           tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTcc           <  1:3625702/70‑1 (MQ=255)                                                            aCGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg         >  1:1288104/1‑71 (MQ=255)                                                             cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt        >  1:3431877/1‑71 (MQ=255)                                                             cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt        >  1:3444359/1‑71 (MQ=255)                                                             cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg         >  1:738033/1‑70 (MQ=255)                                                             cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGATCTGGCAATCa                 <  1:2211954/62‑1 (MQ=255)                                                             agAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt        <  1:2425418/70‑1 (MQ=255)                                                               aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTTGTCTGGCAATCAAAGTCCAgtgt      <  1:432932/71‑1 (MQ=255)                                                               aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt      <  1:2920582/71‑1 (MQ=255)                                                               aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt      <  1:318213/71‑1 (MQ=255)                                                               aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt      <  1:2996298/71‑1 (MQ=255)                                                               aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt      <  1:1728372/71‑1 (MQ=255)                                                               aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt      <  1:2243705/71‑1 (MQ=255)                                                                aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                       >  1:2512439/1‑53 (MQ=255)                                                                aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                      >  1:3434326/1‑54 (MQ=255)                                                                aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc            <  1:937274/64‑1 (MQ=255)                                                                   cGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACt  <  1:1164676/71‑1 (MQ=255)                                                                   |                                                                       GGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACT  >  minE/1202074‑1202210

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