Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,202,149	C→T	6.5%	D434D (GAC→GAT)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,202,149	0	C	T	6.5%	310.6 / 7.5	92	D434D (GAC→GAT)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base T (3/3);  total (46/46)
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.86e-01

CAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCT  >  minE/1202089‑1202219                                                             |                                                                       cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                              <  1:2452681/71‑1 (MQ=255)  aGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                              <  1:1034642/70‑1 (MQ=255)   gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                                      >  1:2380527/1‑61 (MQ=255)   gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                                      >  1:2898195/1‑61 (MQ=255)   gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                            <  1:2126291/71‑1 (MQ=255)   gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                                             >  1:489676/1‑70 (MQ=255)      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                         <  1:2214042/71‑1 (MQ=255)      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                         >  1:1360214/1‑71 (MQ=255)      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                         >  1:2019384/1‑71 (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)                ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACa                                                                       >  1:2464782/1‑47 (MQ=255)                  cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCa                                             >  1:3246103/1‑71 (MQ=255)                   gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAc                                            <  1:2278689/71‑1 (MQ=255)                    cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:1667805/71‑1 (MQ=255)                    cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:901792/71‑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:310864/40‑1 (MQ=17)                              ttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                                              <  1:2694751/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:3602236/1‑71 (MQ=255)                                   gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                            >  1:4121969/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:3562572/1‑70 (MQ=255)                                      aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                          >  1:1993506/1‑70 (MQ=255)                                       cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                          >  1:4007118/1‑69 (MQ=255)                                       cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                          >  1:2792915/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:318213/71‑1 (MQ=255)                                                aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               <  1:1728372/71‑1 (MQ=255)                                                aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               <  1:2996298/71‑1 (MQ=255)                                                aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               <  1:2920582/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)                                                     ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                               <  1:38333/50‑1 (MQ=255)                                                      gTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:269596/37‑1 (MQ=255)                                                      gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                          <  1:941229/54‑1 (MQ=255)                                                      gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa         <  1:3817928/71‑1 (MQ=255)                                                       tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAAt        <  1:2972439/71‑1 (MQ=255)                                                       tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAAt        <  1:3268370/71‑1 (MQ=255)                                                       tAAGGACATGAACGACAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACt           <  1:1528519/68‑1 (MQ=255)                                                        aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                   >  1:1467427/1‑59 (MQ=255)                                                        aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                   >  1:1630602/1‑59 (MQ=255)                                                        aaGGACATGAAAGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                   >  1:2892022/1‑59 (MQ=255)                                                         aGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGAc            <  1:177403/65‑1 (MQ=255)                                                           gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt    >  1:3551318/1‑71 (MQ=255)                                                           gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt    >  1:1671807/1‑71 (MQ=255)                                                             cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCTATCAAAGTCCAGTGTGACTAATActt    <  1:1510744/69‑1 (MQ=255)                                                             cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGt                      >  1:3754114/1‑51 (MQ=255)                                                             cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct  >  1:45749/1‑71 (MQ=255)                                                             cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct  <  1:89454/71‑1 (MQ=255)                                                             |                                                                       CAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCT  >  minE/1202089‑1202219

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