Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,658,871	G→A	12.7%	P357P (CCG→CCA)	ppx →	exopolyphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,658,871	0	G	A	12.7%	167.0 / 12.6	63	P357P (CCG→CCA)	ppx	exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (6/49);  new base A (3/5);  total (9/54)
Fisher's exact test for biased strand distribution p-value = 7.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAC  >  minE/1658802‑1658931                                                                      |                                                             aCGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCTGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:3420278/71‑1 (MQ=255)      tGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTg                                                        <  1:1571405/71‑1 (MQ=255)           gATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgc                                                   <  1:3824303/71‑1 (MQ=255)           gATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgc                                                   <  1:2192599/71‑1 (MQ=255)             tACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgcgc                                                 >  1:5181152/1‑71 (MQ=255)                cACTATGCAAATGTACGAACATTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:1998725/62‑1 (MQ=255)                cACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:4721775/62‑1 (MQ=255)                cACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:3641685/62‑1 (MQ=255)                cACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:4440943/62‑1 (MQ=255)                 aCTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTgg                                                       <  1:2112406/61‑1 (MQ=255)                   tATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTGGAGGCGCTACTGc                                           >  1:475909/1‑71 (MQ=255)                   tATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTGGAGGCGCTACTGc                                           >  1:2633753/1‑71 (MQ=255)                   tATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCACAACTGGAGGCGCTACTGc                                           >  1:3755705/1‑71 (MQ=255)                     tGCAAATGTACGAACAGTGGCGGGAACAGCAACTGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGa                                         >  1:2379050/1‑71 (MQ=255)                     tGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTg                                            <  1:3480246/68‑1 (MQ=255)                       cAAATGTACGAACAGTTGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGAt                                        <  1:518715/70‑1 (MQ=255)                       cAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGAt                                        <  1:4203805/70‑1 (MQ=255)                       cAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:1900269/71‑1 (MQ=255)                       cAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:4420442/71‑1 (MQ=255)                        aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGAt                                        <  1:251663/69‑1 (MQ=255)                        aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATgg                                      <  1:4786114/71‑1 (MQ=255)                        aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:1403103/70‑1 (MQ=255)                        aaaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:2884352/70‑1 (MQ=255)                         aaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATggg                                     <  1:1588765/71‑1 (MQ=255)                         aaTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATggc                                     <  1:2485885/71‑2 (MQ=255)                         aaTGTACGAACAGTGGCGGGAACAGCAACCGAAGATGGCGCATCCGCAACTGGAGGCGCTACTGCGATggg                                     <  1:959902/71‑1 (MQ=255)                          aTGTACGAACAGTGGCGGGAACAGCAACCTAAGCTGGCGCATCCGc                                                             <  1:327614/46‑1 (MQ=255)                          aTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:3133695/46‑1 (MQ=255)                          aTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:1833610/46‑1 (MQ=255)                          aTGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGc                                                             <  1:544182/46‑1 (MQ=255)                                    aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:4892124/49‑1 (MQ=255)                                    aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:5007386/49‑1 (MQ=255)                                    aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:565121/49‑1 (MQ=255)                                    aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:1355999/49‑1 (MQ=255)                                    aGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCt                                                <  1:844200/49‑1 (MQ=255)                                     gTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCAtgctg                          >  1:3739563/1‑70 (MQ=255)                                     gTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCAtgctg                          >  1:3863871/1‑70 (MQ=255)                                        gCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGgc                                                   <  1:4149296/42‑1 (MQ=255)                                          gggAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTATTGCGAt                                        <  1:9836/51‑1 (MQ=255)                                          gggAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:135283/52‑1 (MQ=255)                                          gggAACAGCAACCGAAGCGGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:3409147/52‑1 (MQ=255)                                            gAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt                  <  1:796742/71‑1 (MQ=255)                                             aaCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATg                                       <  1:5013665/49‑1 (MQ=255)                                              aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:4610877/71‑1 (MQ=255)                                              aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:753733/71‑1 (MQ=255)                                              aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:1521391/71‑1 (MQ=255)                                              aCAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:1652537/71‑1 (MQ=255)                                               cAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTg                                            <  1:4785590/42‑1 (MQ=255)                                                 gCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTg                                            <  1:1262017/40‑1 (MQ=255)                                                  cAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt                  <  1:3886468/65‑1 (MQ=255)                                                  cAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCg                <  1:4744732/67‑1 (MQ=255)                                                   aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:1251660/71‑1 (MQ=255)                                                   aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:1515770/71‑1 (MQ=255)                                                   aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:4691287/71‑1 (MQ=255)                                                   aaCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:3589504/71‑1 (MQ=255)                                                    aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGttag          <  1:2150800/71‑3 (MQ=255)                                                    aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGttag          <  1:3098259/71‑3 (MQ=255)                                                    aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGttag          <  1:4432055/71‑3 (MQ=255)                                                    aCCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGa          <  1:4571143/71‑1 (MQ=255)                                                        taGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt                  <  1:3279671/58‑1 (MQ=255)                                                        aaGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTg           <  1:3921552/66‑1 (MQ=255)                                                                     cGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAc  >  1:2924265/1‑62 (MQ=255)                                                                     cGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAc  >  1:1417490/1‑62 (MQ=255)                                                                      |                                                             ACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAATATCAAC  >  minE/1658802‑1658931

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