Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,202,182	T→C	11.9%	G445G (GGT→GGC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,202,182	0	T	C	11.9%	247.5 / 18.9	92	G445G (GGT→GGC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base T (44/37);  new base C (2/9);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

TAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTA  >  Exported/1202113‑1202252                                                                      |                                                                       tAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtg                                                                       <  1:2025406/71‑1 (MQ=255)   aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg                                                                      <  1:5193472/70‑1 (MQ=255)     gttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCt                                                                   <  1:4544821/68‑1 (MQ=12)       cATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                                                 >  1:1888669/1‑71 (MQ=255)        aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                                                                <  1:3744487/71‑1 (MQ=255)          cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                                              >  1:2265199/1‑71 (MQ=255)          cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                                              >  1:759792/1‑71 (MQ=255)          cGAAACGATGTACGAAACCGGTAAGGACATGAACGACAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                                              >  1:5452808/1‑71 (MQ=255)           gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                                                              >  1:160257/1‑70 (MQ=255)           gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                                             >  1:3277154/1‑71 (MQ=255)           gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                                             >  1:5074494/1‑71 (MQ=255)           gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                                             >  1:3688775/1‑71 (MQ=255)           gAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt                                                             <  1:4977722/71‑1 (MQ=17)             aaCCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtga                                                           <  1:3170059/71‑3 (MQ=12)               ccaTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                                         <  1:2610426/69‑2 (MQ=12)                caTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaaga                                                        <  1:761970/70‑3 (MQ=12)                  tGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc                                                      <  1:3174331/71‑1 (MQ=255)                    tACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                                                            <  1:1496302/63‑2 (MQ=17)                      cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaaa                                                         >  1:1115203/1‑64 (MQ=255)                      cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaaa                                                         >  1:2100008/1‑64 (MQ=255)                      cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                                                  >  1:1861590/1‑71 (MQ=255)                      cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                                                  <  1:5020958/71‑1 (MQ=255)                      cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                                                  >  1:48996/1‑71 (MQ=255)                      cGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCaa                                                              >  1:240513/1‑59 (MQ=17)                      cGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCaa                                                              >  1:5320277/1‑59 (MQ=17)                        aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                                <  1:2330294/71‑1 (MQ=255)                        aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                                <  1:2329887/71‑1 (MQ=255)                        aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                                <  1:3417817/71‑1 (MQ=255)                        aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                                >  1:2959377/1‑71 (MQ=255)                              gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                          <  1:1662673/71‑1 (MQ=255)                               tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                          >  1:4780250/1‑70 (MQ=255)                               tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                          >  1:2895223/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:2886856/1‑46 (MQ=255)                                              cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                                                   >  1:4025890/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)                                                  tACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                                         <  1:1328273/36‑2 (MQ=12)                                                  tACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                                         <  1:5181447/36‑2 (MQ=12)                                                  tACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                                         <  1:4040918/36‑2 (MQ=12)                                                  tACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCt                      >  1:2090788/1‑71 (MQ=255)                                                   aCCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                          >  1:3560348/1‑50 (MQ=255)                                                   aCCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                          >  1:4849932/1‑50 (MQ=255)                                                    ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCTATCTGc                    <  1:739819/71‑1 (MQ=255)                                                    ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa                         <  1:1091376/66‑1 (MQ=255)                                                    ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa                         <  1:2010905/66‑1 (MQ=255)                                                    ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGc                    <  1:2620639/71‑1 (MQ=255)                                                      gcgAAACCTCACGCGGTGTTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCaa                  <  1:5013659/71‑1 (MQ=255)                                                      gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                                        >  1:5338672/1‑49 (MQ=255)                                                      gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCaa                  <  1:2207420/71‑1 (MQ=255)                                                       cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                                                    >  1:2993534/1‑36 (MQ=255)                                                       cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                                <  1:5617285/40‑1 (MQ=255)                                                       cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                                <  1:2064428/40‑1 (MQ=255)                                                       cgAAACCTAACGCGGTGGTCTGGCAATCAAAGTCCa                                                    >  1:2164211/1‑36 (MQ=255)                                                         aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGccc                          <  1:2483240/60‑1 (MQ=255)                                                         aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACgg               <  1:1692655/71‑1 (MQ=255)                                                           aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:2737675/71‑1 (MQ=255)                                                           aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:4266285/71‑1 (MQ=255)                                                           aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:3426651/71‑1 (MQ=255)                                                            ccTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATg            <  1:2294783/71‑1 (MQ=255)                                                            ccTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATg            <  1:719396/71‑1 (MQ=255)                                                              tCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                              >  1:5195265/1‑35 (MQ=255)                                                              tCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                              >  1:2376557/1‑35 (MQ=255)                                                                 cgcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                                        >  1:874680/1‑38 (MQ=255)                                                                 cgcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                                        >  1:3197837/1‑38 (MQ=255)                                                                 cgcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCg        <  1:4435621/70‑1 (MQ=255)                                                                   cgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAAt     <  1:4796274/71‑1 (MQ=255)                                                                     gtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:1051568/61‑1 (MQ=255)                                                                     gtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:4335264/61‑1 (MQ=255)                                                                      tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:5036189/60‑1 (MQ=255)                                                                      tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta  <  1:5513497/71‑1 (MQ=255)                                                                      tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta  <  1:1874762/71‑1 (MQ=255)                                                                      |                                                                       TAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTA  >  Exported/1202113‑1202252

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