Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	610,451	G→A	9.5%	P239S (CCG→TCG)	focA ←	formate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	610,451	0	G	A	9.5%	120.0 / 4.0	42	P239S (CCG→TCG)	focA	formate transporter
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base A (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATA  >  minE/610385‑610521                                                                   |                                                                       cGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCa                                                                    <  1:2514293/71‑1 (MQ=255)  gATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAg                                                                   <  1:2836437/71‑1 (MQ=255)         aCCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGt                                                                       >  1:1596267/1‑60 (MQ=255)         aCCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGATGCAGAACCGa                                                             <  1:3091229/70‑1 (MQ=255)             gAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCg                                                        >  1:1778155/1‑71 (MQ=255)             gAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCg                                                        <  1:1840975/71‑1 (MQ=255)             gAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCg                                                        >  1:1703829/1‑71 (MQ=255)             gAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCg                                                        >  1:849380/1‑71 (MQ=255)             gAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTACCGGTGCAGAACCGACTGc                                                         >  1:473834/1‑70 (MQ=255)                   ggTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCaa                                                  <  1:3613616/71‑1 (MQ=255)                     ttATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGAc                                                            <  1:522507/59‑1 (MQ=255)                     ttATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGAc                                                            <  1:2583263/59‑1 (MQ=255)                     ttATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCaaaa                                                >  1:87919/1‑71 (MQ=255)                     ttATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCaaaa                                                >  1:251490/1‑71 (MQ=255)                        tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGATGCAGAACCGACTGCGGTCCAAAATTc                                             <  1:33909/71‑1 (MQ=255)                          aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg                                           >  1:3594554/1‑71 (MQ=255)                             gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCggg                                         >  1:1127845/1‑70 (MQ=255)                             gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCggg                                         >  1:3020070/1‑70 (MQ=255)                              atgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGa                                       >  1:3581976/1‑71 (MQ=255)                              atgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGa                                       <  1:2973041/71‑1 (MQ=255)                                   aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGAt                                      <  1:3306387/67‑1 (MQ=255)                                   aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGAt                                      <  1:1805868/67‑1 (MQ=255)                                   aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGAt                                      <  1:1766727/67‑1 (MQ=255)                                   aTTCATCACGGTCAGGTGAGAAAAATTTTCCGATGCAGaa                                                                 >  1:289015/1‑40 (MQ=255)                                       atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTc                              <  1:3686906/71‑1 (MQ=255)                                       atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTc                              <  1:3669105/71‑1 (MQ=255)                                         caCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCgggg                                        <  1:253641/59‑1 (MQ=255)                                                    gagaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTcc                                            >  1:2379043/1‑44 (MQ=255)                                                    gagaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCgg                                          >  1:471062/1‑46 (MQ=255)                                                         aaaTTTTCCGGTGCAGAACCGCCTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCAt            >  1:912749/1‑71 (MQ=255)                                                         aaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTc                              >  1:2511289/1‑53 (MQ=255)                                                         aaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCAt            >  1:1839230/1‑71 (MQ=255)                                                         aaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCAt            >  1:3408455/1‑71 (MQ=255)                                                         aaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCAt            >  1:630034/1‑71 (MQ=255)                                                         aaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCAt            >  1:693499/1‑71 (MQ=255)                                                            ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTAc                     <  1:2032242/59‑1 (MQ=255)                                                            ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTAc                     <  1:1760838/59‑1 (MQ=255)                                                            ttttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTAc                     <  1:1462292/59‑1 (MQ=255)                                                            ttttCCGATGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCgg         >  1:1768428/1‑71 (MQ=255)                                                             tttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCgg         >  1:3653287/1‑70 (MQ=255)                                                              ttCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCggg        <  1:1018649/70‑1 (MQ=255)                                                                ccGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTc                              <  1:508919/46‑1 (MQ=255)                                                                   gTGCAGAACCGACTGCGGTCCAAAATTCCTGTGATGc                                    <  1:2994301/37‑1 (MQ=255)                                                                   gTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATa  <  1:199757/71‑1 (MQ=255)                                                                   gTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATa  <  1:1189383/71‑1 (MQ=255)                                                                   |                                                                       CGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATA  >  minE/610385‑610521

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