Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,202,161	A→G	17.5%	K438K (AAA→AAG)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,202,161	0	A	G	17.5%	198.7 / 30.0	80	K438K (AAA→AAG)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base G (11/3);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAT  >  minE/1202091‑1202230                                                                       |                                                                      gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                                       >  1:2548114/1‑71 (MQ=255)   gagccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTa                                                                     >  1:1773298/4‑71 (MQ=9)    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                                    >  1:2251408/1‑71 (MQ=255)    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                                    >  1:2919357/1‑71 (MQ=255)    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                                    >  1:1981145/1‑71 (MQ=255)     aCCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                                                   <  1:2127135/71‑1 (MQ=255)      ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcg                                                                  >  1:1472404/1‑71 (MQ=255)         cgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGa                                                               <  1:3355599/71‑1 (MQ=9)         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                                               <  1:3422741/71‑1 (MQ=255)             tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAcc                                                           >  1:3313604/1‑71 (MQ=255)                cGCTGGATAAGGTCATCGAAACGATCTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                                                             >  1:499541/1‑66 (MQ=255)                   tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                                     >  1:3018705/3‑71 (MQ=9)                   tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                                     >  1:2223873/3‑71 (MQ=9)                   tGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                                     >  1:1956194/1‑71 (MQ=255)                    cgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGaaa                                                             >  1:852263/2‑62 (MQ=14)                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                                                    >  1:1606738/1‑70 (MQ=255)                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                                     <  1:58901/69‑1 (MQ=255)                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                                     <  1:3059465/69‑1 (MQ=255)                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                   <  1:2097776/71‑1 (MQ=255)                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                   >  1:1593923/1‑71 (MQ=255)                      aTAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggt                                                  >  1:3418903/1‑71 (MQ=255)                       tAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATCAACGCCAAATACCGCGAAACCTCACGCggtg                                                 >  1:1263373/1‑71 (MQ=255)                       tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg                                                 >  1:263893/1‑69 (MQ=9)                          ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                   <  1:1367076/66‑1 (MQ=255)                          ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTc                                              >  1:935991/1‑71 (MQ=255)                          ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTc                                              <  1:626062/71‑1 (MQ=255)                           gttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCt                                             <  1:2906174/68‑1 (MQ=12)                           gTCATCGAAACTATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                                             <  1:2756947/71‑1 (MQ=255)                           gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                                             <  1:1826045/71‑1 (MQ=255)                           gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                                             <  1:3562811/71‑1 (MQ=255)                            tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                                     >  1:2438567/1‑46 (MQ=255)                            tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                                                   >  1:3059730/1‑48 (MQ=255)                             taTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                                                                 >  1:1111565/2‑49 (MQ=17)                              aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                   >  1:909735/1‑62 (MQ=255)                              aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                           >  1:2350951/1‑70 (MQ=255)                              aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                           <  1:1054546/70‑1 (MQ=255)                               tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                                    >  1:3040285/1‑44 (MQ=255)                               tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                                    >  1:3445904/1‑44 (MQ=255)                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                        >  1:3831990/1‑70 (MQ=255)                                  aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                       >  1:1703542/1‑70 (MQ=255)                                   aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTc                                              <  1:2609792/62‑1 (MQ=255)                                   aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                                     <  1:2277316/71‑1 (MQ=255)                                     cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                                     >  1:121230/1‑69 (MQ=255)                                      caTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt                                       <  1:2067284/65‑1 (MQ=17)                                       aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                                                                 >  1:1203697/1‑39 (MQ=17)                                       aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                                      >  1:2686122/1‑65 (MQ=17)                                       aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                                      >  1:2598076/1‑65 (MQ=17)                                        tGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgcg                                                                >  1:1984711/1‑39 (MQ=17)                                         gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc                                >  1:718376/1‑70 (MQ=255)                                          tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATc                                      <  1:3846907/63‑1 (MQ=255)                                          tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                                     <  1:2684182/64‑1 (MQ=255)                                          tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                              <  1:491799/71‑1 (MQ=255)                                           aCGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtg                                                 <  1:1528441/51‑1 (MQ=255)                                            cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                            >  1:2844133/1‑71 (MQ=255)                                             gAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtg                           >  1:2327289/1‑71 (MQ=255)                                              aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                          <  1:45108/71‑1 (MQ=255)                                              aaaCAGGTAAAGATATGAACGCCAAGTACCGCGaa                                                              >  1:2588943/1‑35 (MQ=9)                                               aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                          <  1:3421238/70‑1 (MQ=255)                                                    gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                              >  1:3270597/1‑61 (MQ=255)                                                    gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                              >  1:198561/1‑61 (MQ=255)                                                    gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACt                      <  1:412999/69‑1 (MQ=255)                                                     tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTTGTCTGGCAATCAAAGTCCAGTGTGACt                      <  1:830853/68‑1 (MQ=255)                                                     tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                    >  1:457045/1‑70 (MQ=255)                                                       aGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAc                 <  1:2996347/71‑1 (MQ=255)                                                       aGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAc                 >  1:1262430/1‑71 (MQ=255)                                                         gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt               <  1:1824390/71‑1 (MQ=255)                                                         gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt               >  1:3215150/1‑71 (MQ=255)                                                           cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                          <  1:1774333/58‑1 (MQ=255)                                                           cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct             >  1:1631395/1‑71 (MQ=255)                                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                             <  1:1213692/54‑1 (MQ=255)                                                             tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                  >  1:2101311/1‑64 (MQ=255)                                                              gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAc                 <  1:2774320/64‑1 (MQ=255)                                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                             <  1:3383363/51‑1 (MQ=255)                                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc          <  1:2159856/70‑1 (MQ=255)                                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt         <  1:3650946/71‑1 (MQ=255)                                                                aCGCCAAATACCGCGAACCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTc        >  1:3297575/1‑71 (MQ=255)                                                                aCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                             <  1:124844/50‑1 (MQ=255)                                                                aCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt         >  1:2964595/1‑70 (MQ=255)                                                                  gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt         >  1:269641/1‑68 (MQ=255)                                                                  gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGc      >  1:1159513/1‑71 (MQ=255)                                                                    cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc          <  1:2135566/65‑1 (MQ=255)                                                                     aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                    >  1:2556333/1‑54 (MQ=255)                                                                     aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa   >  1:2252232/1‑71 (MQ=255)                                                                      aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct             <  1:768172/60‑1 (MQ=255)                                                                      aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct             <  1:772815/60‑1 (MQ=255)                                                                      aaTACCGCGAAACCTCACGCGGTGGTATGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAt  <  1:2852716/71‑1 (MQ=255)                                                                       |                                                                      GTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAT  >  minE/1202091‑1202230

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