Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,978,961	C→T	10.6%	A111T (GCA→ACA)	ygeA ←	predicted racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,978,961	0	C	T	10.6%	132.7 / 6.5	47	A111T (GCA→ACA)	ygeA	predicted racemase
Reads supporting (aligned to +/- strand):  ref base C (14/28);  new base T (2/3);  total (16/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

CTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTA  >  W3110S.gb/2978899‑2979031                                                               |                                                                       cTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCgg                                                                >  1:2573186/1‑71 (MQ=255)    ttCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGTACGTCCGGTgg                                                             <  1:1627824/71‑1 (MQ=255)      ccATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCa                                                           >  1:567371/1‑71 (MQ=255)           gtgtAACGCGTACCCAGCAGCGCCACACTAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGc                                                      <  1:3288653/71‑1 (MQ=255)                aCGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAAtgt                                                 >  1:1051214/1‑71 (MQ=255)                     taaccAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGTACGTCCGGTGGCATCCGCAATgtgt                                               >  1:2397610/4‑68 (MQ=255)                     tACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGTACGTCCGGTGGCATCCGCAATgtgt                                               >  1:2781100/1‑68 (MQ=255)                     tACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAg                                            <  1:2993880/71‑1 (MQ=255)                       cccAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGa                                           <  1:2610030/70‑1 (MQ=255)                       cccAGCAGCGCCACACGAGTAATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGa                                           <  1:1394089/70‑1 (MQ=255)                         cagcagCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGGATCCGCAATGTGTAAGAAAg                                        >  1:1684781/1‑71 (MQ=255)                          agcagcGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAgg                                       <  1:2765463/71‑1 (MQ=255)                          agcagcGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAgg                                       <  1:2469819/71‑1 (MQ=255)                           gcagcGCCACACGAGTCATTCCGGCCCCGGTAATTGCAGGTCCGGTGGCATCCGCAATGTGTAAGAAAgg                                       <  1:4014999/70‑1 (MQ=255)                            cagcGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCa                                     <  1:2045456/71‑1 (MQ=255)                            cagcGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCa                                     <  1:159031/71‑1 (MQ=255)                             agcGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAgg                                       >  1:490427/1‑68 (MQ=255)                             agcGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCa                                     <  1:287961/70‑1 (MQ=255)                             agcGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTCTAAGAAAGGCag                                    <  1:3817750/71‑1 (MQ=255)                                gcCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCagag                                  >  1:2723790/1‑70 (MQ=255)                                    cacGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCa                                     <  1:3772292/63‑1 (MQ=255)                                    cacGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCa                                     <  1:1850949/63‑1 (MQ=255)                                      cGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACg                           >  1:18322/1‑71 (MQ=255)                                       gAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGt                          <  1:1855295/71‑1 (MQ=255)                                       gAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGt                          <  1:3252336/71‑1 (MQ=255)                                         gTCATTCCGGCCCCGGTAATTGTACGTCCGGTGGCATCCGCAATGTGTAAGa                                           <  1:3442770/52‑1 (MQ=255)                                         gTCATTCCGGCCCCGGTAATTGTACGTCCGGTGGCATCCGCAATGTGTAAGa                                           <  1:3347943/52‑1 (MQ=255)                                          taaTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCag                                    <  1:3028810/56‑1 (MQ=255)                                            aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACt                      >  1:2817419/1‑70 (MQ=255)                                            aTTCAGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACt                      >  1:767567/1‑70 (MQ=255)                                             ttCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCag                                    <  1:381774/55‑1 (MQ=255)                                               ccGGCCCCGGTAATTGCACGTCCGCTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAAt                  >  1:3884796/1‑71 (MQ=255)                                                       ggTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCa              <  1:3215479/67‑1 (MQ=255)                                                       ggTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCa              <  1:3428064/67‑1 (MQ=255)                                                       ggTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCa              <  1:3222697/67‑1 (MQ=255)                                                        gTAATTGCACGTCCTGTGGCATCCGCAATGTGTAAGa                                           <  1:463279/37‑1 (MQ=255)                                                        gTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTc                     <  1:70534/59‑1 (MQ=255)                                                        gTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGc         <  1:2813342/71‑1 (MQ=255)                                                           aTTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCa       <  1:133222/70‑1 (MQ=255)                                                            ttGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCa                                     <  1:2287498/39‑1 (MQ=255)                                                            ttGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGt                          <  1:4060880/50‑1 (MQ=255)                                                            ttGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGc               <  1:221351/61‑1 (MQ=255)                                                             tGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCaga                                   >  1:368463/1‑40 (MQ=255)                                                             tGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCaga                                   >  1:3839153/1‑40 (MQ=255)                                                             tGCACGT‑CGGTGGCATCCGCAATGTGTAAGAAAGGCaga                                   >  1:147828/1‑39 (MQ=38)                                                              gCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCa       >  1:576845/1‑67 (MQ=255)                                                              gCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt    <  1:2500294/70‑1 (MQ=255)                                                               cACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTa  >  1:315856/1‑71 (MQ=255)                                                               |                                                                       CTGTTCCATGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTA  >  W3110S.gb/2978899‑2979031

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