Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,202,022	G→A	17.5%	C392Y (TGC→TAC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,202,022	0	G	A	17.5%	190.6 / 31.5	80	C392Y (TGC→TAC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base A (7/7);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAC  >  minE/1201956‑1202088                                                                   |                                                                   aaaTTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTg                                                                    >  1:1177551/1‑67 (MQ=255) aaaTTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTg                                                                    >  1:1900520/1‑67 (MQ=255) aaaTTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAtt                                                                >  1:3305665/1‑71 (MQ=255)  aaTTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAtt                                                                >  1:1882014/1‑70 (MQ=255)  aaTTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAtt                                                                >  1:3582813/1‑70 (MQ=255)       gCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCg                                                           <  1:312184/70‑1 (MQ=255)       gCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCg                                                           <  1:2622517/70‑1 (MQ=255)        cATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTACAtt                                                                <  1:1758767/64‑1 (MQ=255)         aTGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCATGTGCCGTGCATTGAGCGTaa                                                        <  1:2094581/71‑1 (MQ=255)           ggAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATg                                                      >  1:1734748/1‑71 (MQ=255)            gAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGc                                                     <  1:1883287/71‑1 (MQ=255)            gAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGc                                                     >  1:2658606/1‑71 (MQ=255)            gAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTACATTGAGCGTAATGc                                                     >  1:2140392/1‑71 (MQ=255)             aaCACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTACATTGAGCGTAATGcc                                                    <  1:308729/71‑1 (MQ=255)              acacAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCa                                                   >  1:3434963/1‑71 (MQ=255)              acacAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCa                                                   >  1:2042769/1‑71 (MQ=255)               cacaACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCAt                                                  >  1:1546095/1‑71 (MQ=255)                  aaCCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGc                                               >  1:986116/1‑71 (MQ=255)                    ccTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCct                                             <  1:1949728/71‑1 (MQ=255)                    ccTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCct                                             <  1:2025707/71‑1 (MQ=255)                     cTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCg                                                           <  1:2907493/56‑1 (MQ=255)                     cTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCg                                                           <  1:1748141/56‑1 (MQ=255)                     cTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTACa                                                                  >  1:1426178/1‑49 (MQ=255)                     cTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTACa                                                                  >  1:849738/1‑49 (MQ=255)                      ttGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCctct                                           >  1:573451/1‑71 (MQ=255)                       tGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGc                                               >  1:932264/1‑66 (MQ=255)                       tGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtg                                          >  1:3076461/1‑71 (MQ=255)                        ggTTTAACCTGCGACCCTGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgt                                         >  1:1129389/1‑71 (MQ=255)                        ggTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgt                                         >  1:2319212/1‑71 (MQ=255)                         gTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgt                                         <  1:527658/70‑1 (MQ=255)                           ttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgtg                                        >  1:3094384/1‑69 (MQ=255)                           ttAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgtg                                        >  1:2910635/1‑69 (MQ=255)                            tAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAt                                                                 >  1:3292133/1‑43 (MQ=255)                             aaCCTGCGCCCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAgg                                    >  1:1485045/1‑71 (MQ=255)                              aCCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGc                                               >  1:354495/1‑59 (MQ=255)                              aCCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGc                                   <  1:3207237/71‑1 (MQ=255)                                cTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGa                                 >  1:2830294/1‑71 (MQ=255)                                cTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGa                                 >  1:3145325/1‑71 (MQ=255)                                cTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGa                                 >  1:2117527/1‑71 (MQ=255)                                cTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGa                                 >  1:1188615/1‑71 (MQ=255)                                   cGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAtt                               <  1:1816275/70‑1 (MQ=255)                                   cGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAtt                               <  1:3542975/70‑1 (MQ=255)                                   cGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAtt                               <  1:395518/70‑1 (MQ=255)                                      cccGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAAcg                           <  1:343900/71‑1 (MQ=255)                                      cccGGTTGCAGGGCAGGTTCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAAcg                           >  1:1128140/1‑71 (MQ=255)                                         ggTTGGAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAAcgccg                        >  1:1569949/1‑71 (MQ=255)                                         ggTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAAcgccg                        >  1:1963347/1‑71 (MQ=255)                                           ttGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGt                                                          <  1:1474125/35‑1 (MQ=255)                                              cagggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCcgcgcg                    <  1:257462/70‑1 (MQ=255)                                              cagggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCgg                   <  1:655707/71‑1 (MQ=255)                                                gggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGAt                 <  1:1941471/71‑1 (MQ=255)                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCAt                                                  >  1:3587888/1‑37 (MQ=255)                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCAt                                                  >  1:3540082/1‑37 (MQ=255)                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCAt                                                  >  1:3007648/1‑37 (MQ=255)                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCAt                                                  >  1:1346510/1‑37 (MQ=255)                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATg                >  1:2596258/1‑71 (MQ=255)                                                  gcaggTTCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATgg               >  1:602073/1‑71 (MQ=255)                                                   caggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGc              >  1:432120/1‑71 (MQ=255)                                                   caggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGc              >  1:662610/1‑71 (MQ=255)                                                   caggTTCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGc              <  1:2514035/71‑1 (MQ=255)                                                   caggTTCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGc              <  1:854643/71‑1 (MQ=255)                                                    aggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGct             <  1:809717/71‑1 (MQ=255)                                                    aggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGct             <  1:1128198/71‑1 (MQ=255)                                                    aggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGct             <  1:1464498/71‑1 (MQ=255)                                                    aggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGct             <  1:2560131/71‑1 (MQ=255)                                                      gTTCAGGTGCCGTGCATTGATCGTAATGCCATTGCCTCTGTGAAGGCGa                                 <  1:3491106/49‑1 (MQ=255)                                                      gTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGc                                   <  1:1300024/47‑1 (MQ=255)                                                      gTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGctct           <  1:2674185/71‑1 (MQ=255)                                                      gTTCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAtt                               <  1:2379810/51‑1 (MQ=255)                                                      gTTCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAtt                               <  1:837606/51‑1 (MQ=255)                                                      gTTCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAtt                               <  1:3574546/51‑1 (MQ=255)                                                        tCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTa                              >  1:3188653/1‑50 (MQ=255)                                                        tCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTaa                             >  1:2104603/1‑51 (MQ=255)                                                        tCAGGTGCCGTACATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTaa                             >  1:1362028/1‑51 (MQ=255)                                                         cAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTgc         <  1:3195985/70‑1 (MQ=255)                                                         cAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTgc         <  1:2313521/70‑1 (MQ=255)                                                         cAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTATGTGAAGGCGATTAACGCCGCGCGGATGGCTCTgcg        >  1:3122982/1‑71 (MQ=255)                                                            gTGCCGTGCATTGAGCGTAATGCCATTGCCTATGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGcgccg     >  1:2934163/1‑71 (MQ=255)                                                              gCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCa   >  1:3524533/1‑71 (MQ=255)                                                               ccGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCa   >  1:2100785/1‑70 (MQ=255)                                                                 gTGCATTGAACGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATgg               <  1:2813023/56‑1 (MQ=255)                                                                 gTGCATAGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc  >  1:361137/1‑69 (MQ=255)                                                                   |                                                                   AAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAC  >  minE/1201956‑1202088

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