Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	807,506	A→C	12.4%	A66A (GCT→GCG)	ybhC ←	predicted pectinesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	807,506	0	A	C	12.4%	77.6 / 6.4	40	A66A (GCT→GCG)	ybhC	predicted pectinesterase
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base C (2/3);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAA  >  W3110S.gb/807441‑807573                                                                  |                                                                    cccGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCAcc                                                                >  1:1704736/1‑71 (MQ=255)   cgccgGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTg                                                              <  1:1104534/71‑1 (MQ=255)             acaacaAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGa                                                    <  1:82608/71‑1 (MQ=255)             acaacaAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGa                                                    <  1:3630204/71‑1 (MQ=255)              caacaaAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTgcagca                                                                  >  1:3206109/1‑56 (MQ=255)              aaacaaAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTgcagca                                                                  >  1:4285917/2‑56 (MQ=255)                acaaAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCTCGCTGCAGCACCTGGTGTCAGGGATg                                                  >  1:2038450/1‑70 (MQ=255)                  aaaGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATg                                                  >  1:2081060/1‑68 (MQ=255)                   aaGTCAGGTTTCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCCGCACCTGGTGTCAGGGATGCaa                                               >  1:1476816/1‑70 (MQ=255)                   aaGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCCGCACCTGGTGTCAGGGATGCaa                                               >  1:697991/1‑70 (MQ=255)                   aaGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCaa                                               >  1:565715/1‑70 (MQ=255)                    aGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCaaaa                                             <  1:639868/71‑1 (MQ=255)                    aGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCaaaa                                             <  1:1049727/71‑1 (MQ=255)                          gTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAgtgt                                       >  1:1674827/1‑71 (MQ=255)                          gTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAgtgt                                       <  1:3725009/71‑1 (MQ=255)                          gTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAgtgt                                       <  1:2694256/71‑1 (MQ=255)                           ttGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGtt                                      >  1:2029034/1‑71 (MQ=255)                           ttGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGtt                                      <  1:344632/71‑1 (MQ=255)                             gcgcAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCa                                                <  1:3676533/59‑1 (MQ=255)                                 aggcaggGTAATCGGGTAAGGATTCCACGCTGCCGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCAt                                <  1:3393986/71‑1 (MQ=255)                                 aggcaggGTAATCGGGGAAGGATTCCACGCTGCCGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCAt                                <  1:3130671/71‑1 (MQ=255)                                 aggcaggGTAATCGGGGAAGGATTCCACGCTGCCGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCAt                                <  1:4065642/71‑1 (MQ=255)                                     aggGTAATCGGGGAAGGATTCCACGCTGCAGCAcc                                                                >  1:1730846/1‑35 (MQ=255)                                      gggTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCaa                                               <  1:2186359/51‑1 (MQ=255)                                       ggTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAAt                          >  1:992884/1‑71 (MQ=255)                                       ggTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAAt                          >  1:987830/1‑71 (MQ=255)                                         tAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTc                        >  1:2876225/1‑71 (MQ=255)                                         tAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTc                        >  1:352079/1‑71 (MQ=255)                                           aTCGGTGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCtt                      <  1:4322004/71‑2 (MQ=255)                                              gtggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCaaa                                              >  1:2134322/3‑44 (MQ=255)                                              ggggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGtt                                      >  1:3406912/1‑52 (MQ=255)                                              ggggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGa                                    >  1:3180042/1‑54 (MQ=255)                                               gggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCt                       >  1:847468/1‑66 (MQ=255)                                               gggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCt                       >  1:758821/1‑66 (MQ=255)                                                ggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGctt                 >  1:421703/1‑71 (MQ=255)                                                   aGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGcttctt              <  1:2895289/71‑1 (MQ=255)                                                   aGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGcttctt              <  1:1417602/71‑1 (MQ=255)                                                    ggATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCttt             <  1:1771089/71‑1 (MQ=255)                                                     gATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCtttt            <  1:2108407/71‑1 (MQ=255)                                                        tCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGcc         >  1:369501/1‑71 (MQ=255)                                                             gCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAg    <  1:492771/71‑1 (MQ=255)                                                               tGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGaa  >  1:4076088/1‑71 (MQ=255)                                                                  agcaCCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGtaa                           <  1:2384156/43‑1 (MQ=255)                                                                  |                                                                    CCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAA  >  W3110S.gb/807441‑807573

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