Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,202,095	A→G	13.1%	A416A (GCA→GCG) ‡	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,202,095	0	A	G	13.1%	198.5 / 16.9	76	A416A (GCA→GCG) ‡	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base G (3/7);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATA  >  Exported/1202027‑1202163                                                                     |                                                                     gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                                      >  1:1055199/1‑69 (MQ=255) gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                                      >  1:5573887/1‑69 (MQ=255) gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCAcc                                                                    >  1:580349/1‑71 (MQ=255)     gTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgc                                                                <  1:58668/71‑1 (MQ=255)     gTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgc                                                                <  1:5164087/71‑1 (MQ=255)      tAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgcg                                                               >  1:880853/1‑71 (MQ=255)          gCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                           >  1:3739672/1‑71 (MQ=255)          gCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                           >  1:3333967/1‑71 (MQ=255)           ccATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc                                                           >  1:4229551/1‑70 (MQ=255)           ccATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg                                                          >  1:2211402/1‑71 (MQ=255)             aTTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg                                                          <  1:845554/69‑1 (MQ=255)               tGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTg                                                       >  1:3480006/1‑70 (MQ=255)               tGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTg                                                       >  1:611908/1‑70 (MQ=255)                    ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                                      >  1:4869360/1‑50 (MQ=255)                    ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAg                                                 >  1:4865052/1‑71 (MQ=255)                    ctGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAg                                                 >  1:4823271/1‑71 (MQ=255)                     tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                                >  1:1298988/1‑71 (MQ=255)                     tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                                >  1:4900608/1‑71 (MQ=255)                     tgtgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                                >  1:3969631/1‑71 (MQ=255)                       tgAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTc                                              <  1:3333494/71‑1 (MQ=255)                         aaGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCa                                             >  1:1598500/1‑70 (MQ=255)                         aaGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCAt                                            <  1:1483214/71‑1 (MQ=255)                          aGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATc                                           <  1:1486707/71‑1 (MQ=255)                             cGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa                                        <  1:1944295/71‑1 (MQ=255)                             cGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaa                                        <  1:1098307/71‑1 (MQ=255)                                ttAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAAc                                      >  1:3092243/1‑70 (MQ=255)                                 tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                    >  1:573100/1‑71 (MQ=255)                                  aaCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                    <  1:2635090/70‑1 (MQ=255)                                  aaCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                   >  1:4061821/1‑71 (MQ=255)                                  aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                                    <  1:2338470/70‑3 (MQ=9)                                  aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                                    <  1:2331851/70‑3 (MQ=9)                                  aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAcca                                    >  1:1675178/1‑68 (MQ=9)                                   aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                   <  1:5636342/70‑1 (MQ=255)                                   aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                   <  1:1416027/70‑1 (MQ=255)                                   aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                  <  1:214850/71‑1 (MQ=255)                                   aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                  >  1:1591409/1‑71 (MQ=255)                                    cgccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAAc                                      <  1:1318944/66‑1 (MQ=255)                                      ccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                    >  1:2896731/1‑66 (MQ=255)                                       cgcacgTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACg                              <  1:2094673/71‑1 (MQ=12)                                        gcgcgGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                             <  1:3272937/71‑1 (MQ=255)                                           cgtaTGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAAc                          <  1:3611663/68‑1 (MQ=12)                                             gATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCg                        >  1:3287383/1‑71 (MQ=255)                                                ggCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGt                      >  1:1703545/1‑70 (MQ=255)                                                 gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                     <  1:5028452/70‑1 (MQ=255)                                                 gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                     <  1:374468/70‑1 (MQ=255)                                                 gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                     <  1:602788/70‑1 (MQ=255)                                                  ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCg                                          <  1:1429092/48‑1 (MQ=255)                                                   gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag                  <  1:155195/70‑3 (MQ=12)                                                   gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag                  <  1:3752269/70‑3 (MQ=12)                                                    ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaa                    >  1:4658923/1‑68 (MQ=16)                                                    ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaa                    >  1:1338565/1‑68 (MQ=16)                                                    ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                 >  1:5130055/1‑71 (MQ=255)                                                     tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                               <  1:5138410/56‑1 (MQ=255)                                                     tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                               <  1:279462/56‑1 (MQ=255)                                                       cgccgtACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata               <  1:1741280/70‑3 (MQ=9)                                                       cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                             >  1:4681059/1‑56 (MQ=255)                                                       cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt              <  1:5233609/71‑1 (MQ=255)                                                        gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:3386637/71‑1 (MQ=255)                                                        gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:1053086/71‑1 (MQ=255)                                                        gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:2945933/71‑1 (MQ=255)                                                         ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:485031/71‑1 (MQ=255)                                                         ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:740483/71‑1 (MQ=255)                                                         ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:3133687/71‑1 (MQ=255)                                                         ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa            <  1:2237289/71‑1 (MQ=255)                                                          cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                       >  1:4772866/1‑43 (MQ=255)                                                          cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                       >  1:1721249/1‑43 (MQ=255)                                                          cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg             <  1:122387/69‑1 (MQ=255)                                                          cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa           >  1:2482802/1‑71 (MQ=255)                                                            cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg         >  1:3152718/1‑71 (MQ=255)                                                            cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg         >  1:2128450/1‑71 (MQ=255)                                                              ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc       <  1:2282120/71‑1 (MQ=255)                                                              ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc       <  1:1873247/71‑1 (MQ=255)                                                               cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc       >  1:4645873/1‑70 (MQ=255)                                                                aGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTa                                >  1:3994514/1‑44 (MQ=255)                                                                 gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa    <  1:5386880/71‑1 (MQ=255)                                                                 gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa    <  1:2388045/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)                                                                     |                                                                     GAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATA  >  Exported/1202027‑1202163

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑