Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,748,488	A→G	11.1%	P25P (CCA→CCG)	ydhR →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,748,488	0	A	G	11.1%	112.3 / 5.8	45	P25P (CCA→CCG)	ydhR	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/19);  new base G (3/2);  total (24/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.94e-01

TTGTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGAAAAGAA  >  W3110S.gb/1748423‑1748553                                                                  |                                                                  ttGTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACtt                                                                >  1:1158393/1‑69 (MQ=255) ttGTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACtt                                                                >  1:926271/1‑69 (MQ=255) ttGTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACtt                                                                >  1:2433970/1‑69 (MQ=255)  tGTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCt                                                             >  1:161597/1‑71 (MQ=255)   tttACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCt                                                             <  1:2514530/69‑1 (MQ=255)   gTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCt                                                             >  1:3590232/1‑70 (MQ=255)   gTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCt                                                             <  1:2454148/70‑1 (MQ=255)          cTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGAt                                                     <  1:385528/71‑1 (MQ=255)          cTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGAt                                                     <  1:2678052/71‑1 (MQ=255)             cATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGAtt                                                    <  1:3140461/69‑1 (MQ=255)             cATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGAtt                                                    <  1:1186439/69‑1 (MQ=255)              aTTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAAt                                                 <  1:84957/71‑1 (MQ=255)              aTTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAAt                                                 <  1:386483/71‑1 (MQ=255)                         aaTGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAAt                                                 >  1:3451802/1‑60 (MQ=255)                             gCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGtttt                                  <  1:1844052/71‑1 (MQ=255)                              cccATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTc                                 >  1:1731364/1‑71 (MQ=255)                                  tttGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGa                                            <  1:527911/56‑1 (MQ=255)                                  tttGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCtgtg                             >  1:1525000/1‑71 (MQ=255)                                    tGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCa                                               >  1:899274/1‑51 (MQ=255)                                    tGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTgg                            >  1:749374/1‑70 (MQ=255)                                       cGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAgg                        <  1:1363834/71‑1 (MQ=255)                                       cGACGCAATGACTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAgg                        <  1:671915/71‑1 (MQ=255)                                          cGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTAt                     >  1:78188/1‑71 (MQ=255)                                          cGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTAt                     >  1:1121326/1‑71 (MQ=255)                                           gCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCtg                               >  1:2781418/1‑60 (MQ=255)                                           gCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATg                    <  1:1696874/71‑1 (MQ=255)                                             aaTGGCTGAGCAGCTTAAACCGCTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGt                       <  1:80148/66‑1 (MQ=255)                                             aaTGGCTGAGCAGCTTAAACCGCTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGt                       <  1:448331/66‑1 (MQ=255)                                              aTGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGa                  >  1:3584684/1‑70 (MQ=255)                                               tGGCTGAGCAGCTTAAACCACTTGCTAAGTCGATTAAt                                                 <  1:3561029/38‑1 (MQ=255)                                                ggCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAAc                                          >  1:2299231/1‑44 (MQ=255)                                                ggCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGaa                          >  1:209959/1‑60 (MQ=255)                                                 gCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTAt                     >  1:2728583/1‑64 (MQ=255)                                                 gCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAg               >  1:1302341/1‑70 (MQ=255)                                                       cagcTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGa        <  1:2440982/71‑1 (MQ=255)                                                         gcTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAgg                        <  1:1017327/53‑1 (MQ=255)                                                          cTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGaa       <  1:94351/69‑1 (MQ=255)                                                           ttAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTa                      >  1:1886719/1‑53 (MQ=255)                                                            tAAACCGCTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCt                                >  1:2459537/1‑42 (MQ=255)                                                            tAAACCGCTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCt                                >  1:3350180/1‑42 (MQ=255)                                                            tAAACCGCTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCt                                >  1:3134349/1‑42 (MQ=255)                                                            tAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGaaaa     >  1:2199439/1‑69 (MQ=255)                                                            tAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGAAAAg    >  1:796410/1‑70 (MQ=255)                                                             aaaCCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGAAAAg    <  1:289299/69‑1 (MQ=255)                                                              aaCCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGAAAAGaa  <  1:3330873/70‑1 (MQ=255)                                                                  |                                                                  TTGTTACAACTTCATTTTGCTTTTAATGGCCCATTTGGCGACGCAATGGCTGAGCAGCTTAAACCACTTGCTGAGTCGATTAATCAGGAACCTGGTTTTCTGTGGAAGGTATGGACAGAAAGTGAAAAGAA  >  W3110S.gb/1748423‑1748553

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑