Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	381,851	A→C	G88G (GGT→GGG)	yaiF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	381,851	0	A	C	97.9%	166.8 / ‑4.9	47	G88G (GGT→GGG)	yaiF	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (21/25);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

AGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGC  >  W3110S.gb/381786‑381910                                                                  |                                                            aGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTg                                                        <  1:967292/71‑1 (MQ=255)    tGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAAt                                                     <  1:3436041/71‑1 (MQ=255)     gTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAAt                                                     >  1:2054315/1‑70 (MQ=255)     gTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAAt                                                     >  1:311970/1‑70 (MQ=255)          gAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGc                                               <  1:1120546/71‑1 (MQ=255)          gAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGc                                               <  1:1635671/71‑1 (MQ=255)                    ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaa                                     <  1:2736464/71‑1 (MQ=255)                    ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaa                                     <  1:2846603/71‑1 (MQ=255)                    ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaa                                     <  1:3991627/71‑1 (MQ=255)                    ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaa                                     <  1:3807325/71‑1 (MQ=255)                            cAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACggg                              <  1:1149817/70‑1 (MQ=255)                            cAATTACCTATCAGGCAGTTTGCGCCAATCACTACCAGCCCCTGAATAACGGCACCATGACAAATACggg                              <  1:2457422/70‑1 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGaa                                                      >  1:4129529/1‑44 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGaa                                                      >  1:2334424/1‑44 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGaa                                                      >  1:2284235/1‑44 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTa                            >  1:1921708/1‑70 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTa                            >  1:2572630/1‑70 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTAt                           >  1:2055309/1‑71 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTAt                           >  1:59458/1‑71 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTAt                           >  1:1361011/1‑71 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTAt                           >  1:628554/1‑71 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTAt                           >  1:837801/1‑71 (MQ=255)                              aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTAt                           >  1:1276155/1‑71 (MQ=255)                                 aCCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTcg                        >  1:3235192/1‑71 (MQ=255)                                  ccTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTcgc                       <  1:2682084/71‑1 (MQ=255)                                        aGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCa                                           >  1:3270086/1‑45 (MQ=255)                                        aGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGa                  >  1:329159/1‑70 (MQ=255)                                        aGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGa                  >  1:1377315/1‑70 (MQ=255)                                          gCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACgg                                                >  1:112930/1‑38 (MQ=255)                                          gCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaa               >  1:3617151/1‑71 (MQ=255)                                          gCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaa               >  1:147564/1‑71 (MQ=255)                                           cAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                             >  1:1660292/1‑40 (MQ=255)                                             gTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaa            <  1:555725/71‑1 (MQ=255)                                              tttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaaa                                    <  1:2789675/46‑1 (MQ=255)                                              tttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaaa                                    <  1:686906/46‑1 (MQ=255)                                              tttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaaa                                    <  1:326597/46‑1 (MQ=255)                                              tttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaac              <  1:2526009/68‑1 (MQ=255)                                              tttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaac              <  1:1920372/68‑1 (MQ=255)                                              tttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaac              <  1:2087957/68‑1 (MQ=255)                                                       aTCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc  <  1:1436473/71‑1 (MQ=255)                                                       aTCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc  <  1:1296961/71‑1 (MQ=255)                                                       aTCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc  <  1:1292118/71‑1 (MQ=255)                                                       aTCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc  <  1:2109693/71‑1 (MQ=255)                                                        tCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaac           <  1:3822555/61‑1 (MQ=255)                                                        tCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaac           <  1:3817758/61‑1 (MQ=255)                                                        tCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaac           <  1:727502/61‑1 (MQ=255)                                                              ccGGCCCCTGAATAACGGCACCATGACAAATACGGGTa                            >  1:955372/1‑38 (MQ=255)                                                                  |                                                            AGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGC  >  W3110S.gb/381786‑381910

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