Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	494,477	C→T	17.0%	L272L (CTG→CTA)	galM ←	galactose‑1‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	494,477	0	C	T	17.0%	117.6 / 15.7	47	L272L (CTG→CTA)	galM	galactose‑1‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base T (4/4);  total (26/21)
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.61e-01

GCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGT  >  minE/494410‑494543                                                                    |                                                                   gcgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGc                                                                 >  1:663473/1‑71 (MQ=255) gcgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGc                                                                 >  1:2230171/1‑71 (MQ=255)  cgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCt                                                                >  1:3565559/1‑71 (MQ=255)   gACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTg                                                               >  1:3945651/1‑71 (MQ=255)    aCGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGc                                                              <  1:2376437/71‑1 (MQ=255)         gtgCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAAttt                                                         >  1:3537249/1‑71 (MQ=255)         gtgCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAAtt                                                          >  1:1848858/1‑70 (MQ=255)         gtgCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAAtt                                                          >  1:2286259/1‑70 (MQ=255)            ccgccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAAtt                                                          <  1:3419475/67‑1 (MQ=255)            ccgccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTc                                                      <  1:2601509/71‑1 (MQ=255)              gccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCAt                                                    <  1:2807820/71‑1 (MQ=255)                cgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGTAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATc                                                   >  1:1323859/1‑70 (MQ=255)                cgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTTAGCTGCAATTTTTCATct                                                  >  1:2499339/1‑71 (MQ=255)                 gAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATctg                                                 >  1:1108520/1‑71 (MQ=255)                          gCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAg                                         <  1:1731389/70‑1 (MQ=255)                           ccGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTTAGCTGCAATTTTTCATCTGCTGACCAGa                                        <  1:4115551/70‑1 (MQ=255)                           ccGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTTAGCTGCAATTTTTCATCTGCTGACCAGa                                        <  1:253561/70‑1 (MQ=255)                           ccGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGAc                                       <  1:3573604/71‑1 (MQ=255)                           ccGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGAc                                       <  1:228394/71‑1 (MQ=255)                            cGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCa                                          >  1:1720356/1‑67 (MQ=255)                              gAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATg                                    >  1:3818684/1‑71 (MQ=255)                                gTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATgcg                                  >  1:679336/1‑71 (MQ=255)                                 tAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCAt                                                    <  1:2074995/52‑1 (MQ=255)                                   gAATTGCAGAGCCGGAGCGGTGGTGTATACCTTCAGCTGCAATTTTTCATCTGCTGACCa                                          >  1:210227/1‑60 (MQ=255)                                   gAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTTAGCTGCAATTTTTCATCTGCTGACCAGACATGcgccg                               >  1:3594088/1‑71 (MQ=255)                                   gAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTc                                                      <  1:1125178/48‑1 (MQ=255)                                   gAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGcgcc                                >  1:1763043/1‑70 (MQ=255)                                   gAATTGCAGAGCCGAAGCGGTGGTGTAGACCTTTAGCTGCAATTTTTCATCTGCTGACCAGACATGcgccg                               >  1:1855022/1‑71 (MQ=255)                                    aaTTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGcgccgc                              <  1:2588569/71‑1 (MQ=255)                                    aaTTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGcgccgc                              <  1:2416001/71‑1 (MQ=255)                                     aTTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGcgccgc                              >  1:3189601/1‑70 (MQ=255)                                      ttGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATct                                                  >  1:2410548/1‑49 (MQ=255)                                        gCAGAGCCGGAGCGGTGGTGTAGACCTTTAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACt                          >  1:2271856/1‑71 (MQ=255)                                           gagCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTc                       >  1:836995/1‑71 (MQ=255)                                           gagCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTc                       >  1:502638/1‑71 (MQ=255)                                            agCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTc                       <  1:3001490/70‑1 (MQ=255)                                              ccGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGAc                                            <  1:1198838/47‑1 (MQ=255)                                                ggAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGcc                  <  1:2634074/71‑1 (MQ=255)                                                  aGCGGTGGTGTAGACCTTTAGCTGCAATTTTTCATCTGCTGACCAGa                                        <  1:245648/47‑1 (MQ=255)                                                      gtggtgTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGcc            <  1:3984593/71‑1 (MQ=255)                                                            tAGACCTTTAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTc                       <  1:3845236/54‑1 (MQ=255)                                                            tAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACtt                         <  1:1667525/52‑1 (MQ=255)                                                            tAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCt           <  1:2181060/66‑1 (MQ=255)                                                             aGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTg        >  1:2182058/1‑68 (MQ=255)                                                                 cTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGt  >  1:2543853/1‑70 (MQ=255)                                                                   tCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGc             >  1:423265/1‑57 (MQ=255)                                                                   tCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGc             >  1:676970/1‑57 (MQ=255)                                                                    |                                                                   GCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGT  >  minE/494410‑494543

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