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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,202,161	0	A	G	13.4%	215.5 / 19.7	82	K438K (AAA→AAG)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base G (5/6);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAT  >  Exported/1202091‑1202230                                                                       |                                                                      gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                                       <  1:2388045/71‑1 (MQ=255) gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                                       <  1:5386880/71‑1 (MQ=255)  tGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAAt                                                                      <  1:3202003/71‑1 (MQ=255)   gCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                                     >  1:1085806/1‑71 (MQ=255)    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                                     <  1:588897/70‑1 (MQ=255)         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcgcg                                                                >  1:3409604/1‑70 (MQ=255)         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcgcg                                                                >  1:2217578/1‑70 (MQ=255)         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                                               <  1:3744168/71‑1 (MQ=255)         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                                               <  1:4411502/71‑1 (MQ=255)         cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                                               <  1:4482551/71‑1 (MQ=255)             tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcgc                                                                 <  1:5520954/65‑1 (MQ=255)               tcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                                       >  1:3928558/1‑57 (MQ=255)                cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCa                                                        >  1:4829874/1‑71 (MQ=255)                  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                                       >  1:2878054/1‑54 (MQ=255)                  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                                      >  1:5574699/1‑71 (MQ=255)                  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                                      <  1:4143883/71‑1 (MQ=255)                  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                                      <  1:3343172/71‑1 (MQ=255)                  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                                      <  1:466880/71‑1 (MQ=255)                  cTGGATAAGGGCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                                      <  1:2006494/71‑1 (MQ=255)                   tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                                     >  1:4221071/3‑71 (MQ=9)                   tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                                     >  1:1473167/3‑71 (MQ=9)                    ggATAAGGTCATCGAAACGATGTACGAAACCGTTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                                     >  1:5586504/1‑70 (MQ=255)                       tAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtg                                                 <  1:2025406/71‑1 (MQ=255)                         aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                                                    <  1:3738303/66‑1 (MQ=255)                         aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg                                                <  1:5193472/70‑1 (MQ=255)                           gttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCt                                             <  1:4544821/68‑1 (MQ=12)                            tCATCGAAACGATGTAAGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                                      <  1:5067371/61‑1 (MQ=255)                             cATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                           >  1:1888669/1‑71 (MQ=255)                              aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                                          <  1:3744487/71‑1 (MQ=255)                                cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaa                                                              <  1:5697285/49‑1 (MQ=255)                                cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                        >  1:2265199/1‑71 (MQ=255)                                cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                        >  1:759792/1‑71 (MQ=255)                                cGAAACGATGTACGAAACCGGTAAGGACATGAACGACAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                        >  1:5452808/1‑71 (MQ=255)                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                                     >  1:4248574/1‑41 (MQ=255)                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                        >  1:160257/1‑70 (MQ=255)                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                       >  1:3277154/1‑71 (MQ=255)                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                       >  1:3688775/1‑71 (MQ=255)                                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                       >  1:5074494/1‑71 (MQ=255)                                 gAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt                                       <  1:4977722/71‑1 (MQ=17)                                  aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                   >  1:2875460/1‑58 (MQ=255)                                   aaCCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtga                                     <  1:3170059/71‑3 (MQ=12)                                     ccaTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                   <  1:2610426/69‑2 (MQ=12)                                      caTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaaga                                  <  1:761970/70‑3 (MQ=12)                                        tGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaa                                                              <  1:1006210/41‑1 (MQ=255)                                        tGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc                                <  1:3174331/71‑1 (MQ=255)                                          tACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCgg                                                   >  1:499273/1‑50 (MQ=17)                                          tACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                                      <  1:1496302/63‑2 (MQ=17)                                            cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaaa                                   >  1:1115203/1‑64 (MQ=255)                                            cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaaa                                   >  1:2100008/1‑64 (MQ=255)                                            cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                            >  1:48996/1‑71 (MQ=255)                                            cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                            <  1:5020958/71‑1 (MQ=255)                                            cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                            >  1:1861590/1‑71 (MQ=255)                                            cGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCaa                                        >  1:240513/1‑59 (MQ=17)                                            cGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCaa                                        >  1:5320277/1‑59 (MQ=17)                                              aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                          <  1:3417817/71‑1 (MQ=255)                                              aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                          >  1:2959377/1‑71 (MQ=255)                                              aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                          <  1:2329887/71‑1 (MQ=255)                                              aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                          <  1:2330294/71‑1 (MQ=255)                                                    gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                    <  1:1662673/71‑1 (MQ=255)                                                     tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                    >  1:2895223/1‑70 (MQ=255)                                                     tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                    >  1:4780250/1‑70 (MQ=255)                                                      aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGAc                       >  1:1363211/1‑66 (MQ=255)                                                      aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                  >  1:4117496/1‑71 (MQ=255)                                                      aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                  >  1:721408/1‑71 (MQ=255)                                                         gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt               <  1:5143428/71‑1 (MQ=255)                                                           cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGTCAATCAAAGTCCAGTGTGACTaa                    <  1:2364631/64‑1 (MQ=255)                                                           cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct             <  1:4912416/71‑1 (MQ=255)                                                           cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct             <  1:3857742/71‑1 (MQ=255)                                                           cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtc              <  1:1245795/70‑1 (MQ=255)                                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTgaga                        >  1:1917519/1‑59 (MQ=255)                                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                        >  1:4426291/1‑59 (MQ=255)                                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                        >  1:850636/1‑59 (MQ=255)                                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                        >  1:1910741/1‑59 (MQ=255)                                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                    >  1:600639/1‑63 (MQ=255)                                                             tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtctt            <  1:131133/70‑1 (MQ=255)                                                             tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTa           >  1:907824/1‑71 (MQ=255)                                                              gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTa           >  1:2953457/1‑70 (MQ=255)                                                              gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc          >  1:2256805/1‑71 (MQ=255)                                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtctt            <  1:134453/68‑1 (MQ=255)                                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc          <  1:1087350/70‑1 (MQ=255)                                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt         >  1:5120370/1‑71 (MQ=255)                                                                  gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTcc                               <  1:3840604/46‑1 (MQ=255)                                                                   ccAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                              <  1:3666635/46‑1 (MQ=255)                                                                    cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                             >  1:4025890/1‑46 (MQ=255)                                                                    cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                             >  1:2886856/1‑46 (MQ=255)                                                                    cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGccc    <  1:1561643/71‑1 (MQ=255)                                                                     aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa   >  1:201241/1‑71 (MQ=255)                                                                      aaTACCGCGAAACCTCACGCGGTTGTCTGGCAATCAAAGTCCAgt                            <  1:1155662/45‑1 (MQ=255)                                                                      aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa   >  1:1569990/1‑70 (MQ=255)                                                                      aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAt  >  1:514951/1‑71 (MQ=255)                                                                       aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                                     <  1:463645/35‑1 (MQ=255)                                                                       aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                  >  1:2129067/1‑54 (MQ=255)                                                                       aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                  >  1:2004569/1‑54 (MQ=255)                                                                       |                                                                      GTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAT  >  Exported/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: ‑