Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,736	C→G	16.4%	A372A (GCG→GCC)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,736	0	C	G	16.4%	160.9 / 20.7	67	A372A (GCG→GCC)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (21/35);  new base G (4/7);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCC  >  minE/2093669‑2093805                                                                    |                                                                      gTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgcg                                                                    <  1:1594202/71‑1 (MQ=255) gTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                                      <  1:1090986/69‑1 (MQ=255) gTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                                      <  1:2607773/69‑1 (MQ=255)     aaTGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                                                 <  1:1763668/70‑1 (MQ=255)     aaTGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                                                 <  1:1713354/70‑1 (MQ=255)       tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                                              <  1:3704535/71‑1 (MQ=255)         cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                                            <  1:327095/71‑1 (MQ=255)         cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                                            <  1:482353/71‑1 (MQ=255)         cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccg                                                             <  1:385379/70‑1 (MQ=255)           cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTc                                                               <  1:2370356/66‑1 (MQ=255)             ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                                      >  1:3059344/1‑57 (MQ=255)               aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                                      >  1:1312345/1‑71 (MQ=255)               aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                                      <  1:755610/71‑1 (MQ=255)               aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTACGCCGCAtt                                                      <  1:3689283/71‑1 (MQ=255)                tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                                              <  1:188381/62‑1 (MQ=255)                 ttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                                     <  1:2937035/70‑1 (MQ=255)                  tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                                     <  1:2619489/69‑1 (MQ=255)                   gTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgc                                                                     >  1:1220318/1‑52 (MQ=255)                    tACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCac                                                  <  1:2748418/70‑2 (MQ=255)                    tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                  <  1:1846827/70‑1 (MQ=255)                    tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                  <  1:660514/70‑1 (MQ=255)                    tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                  <  1:665244/70‑1 (MQ=255)                      ccTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACAcc                                               <  1:621551/71‑1 (MQ=255)                      ccTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACAcc                                               <  1:3773520/71‑1 (MQ=255)                       cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCt                                              <  1:3419973/71‑1 (MQ=255)                       cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCt                                              <  1:3625457/71‑1 (MQ=255)                          tCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTg                                             <  1:3655475/69‑1 (MQ=9)                            cGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                          <  1:2887763/70‑1 (MQ=255)                              gCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgc                                                                     <  1:212616/41‑1 (MQ=255)                              gCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgc                                                                     <  1:1494542/41‑1 (MQ=255)                              gCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCAATCGCGATTTCCGCCGCATTGCATACCTGCGCCg                                        >  1:3141530/1‑70 (MQ=255)                               cgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCggg                                      <  1:1670782/69‑2 (MQ=255)                               cgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCggg                                      <  1:2687725/69‑2 (MQ=255)                               cAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                                      <  1:1213022/39‑1 (MQ=255)                                aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                     <  1:1732695/71‑1 (MQ=255)                                aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                     >  1:513882/1‑71 (MQ=255)                                aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                     <  1:3014455/71‑1 (MQ=255)                                 aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                    >  1:2419385/1‑71 (MQ=255)                                 aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                    >  1:1686890/1‑71 (MQ=255)                                  ccGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCGGGCtt                                   >  1:2433961/1‑70 (MQ=255)                                  ccGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAt                                                                  <  1:610524/40‑1 (MQ=255)                                   cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTgc                                  >  1:613875/1‑71 (MQ=255)                                     gATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                                                 >  1:556891/1‑38 (MQ=255)                                      aTCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                                              <  1:3048837/40‑1 (MQ=255)                                         gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGccc                             <  1:838594/70‑1 (MQ=255)                                          cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag                           >  1:2013484/1‑71 (MQ=255)                                             gTCAGACCGAGGTTGTGCTCCATGGCGATTTccgcc                                                           >  1:179941/1‑36 (MQ=19)                                             gTCAGACCGAGGTTGTGCTCCATGGCGATTTccgcc                                                           >  1:748225/1‑36 (MQ=19)                                             gTCAGACCGAGGTTGTGCTCCATGGCGATTTccgcc                                                           >  1:2021829/1‑36 (MQ=19)                                              tCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt                       <  1:103575/71‑1 (MQ=255)                                                aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTTCATACCTGCGCCGGACTGCCGCCCAGTAGTTc                     <  1:1713055/71‑1 (MQ=255)                                                aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                      >  1:328614/1‑70 (MQ=255)                                                 gCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCa                    <  1:1277156/71‑1 (MQ=255)                                                  cccAAGGTTATGGTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                          >  1:2891604/1‑48 (MQ=255)                                                    cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCa                    <  1:3106494/68‑1 (MQ=255)                                                    cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                 <  1:2743762/71‑1 (MQ=255)                                                    cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGt                  >  1:1230071/1‑70 (MQ=255)                                                     aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagt                          <  1:3016703/61‑1 (MQ=255)                                                      aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                          >  1:833009/1‑44 (MQ=255)                                                      aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa               >  1:143060/1‑71 (MQ=255)                                                      aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                >  1:76089/1‑70 (MQ=255)                                                       ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAc              <  1:265030/71‑1 (MQ=255)                                                        gTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGtt                      <  1:3850049/62‑1 (MQ=255)                                                        gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTg                                   >  1:1920682/1‑49 (MQ=255)                                                         ttATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAccc            <  1:306517/71‑1 (MQ=255)                                                           aTGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCt           <  1:2570331/70‑1 (MQ=255)                                                            tGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTg          <  1:995383/70‑1 (MQ=255)                                                             gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGcc                                         <  1:2275269/38‑1 (MQ=255)                                                             gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                   >  1:2008961/1‑60 (MQ=255)                                                               tCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgc      >  1:2338377/1‑71 (MQ=255)                                                                ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc     >  1:3727431/1‑71 (MQ=255)                                                                ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc     >  1:902015/1‑71 (MQ=255)                                                                 cATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccg    >  1:115773/1‑71 (MQ=255)                                                                  aTCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgcc        >  1:1887974/1‑66 (MQ=255)                                                                   tCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgcc  <  1:2749942/71‑1 (MQ=255)                                                                    |                                                                      GTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCC  >  minE/2093669‑2093805

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