Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,179,239	C→T	16.2%	L55L (CTG→TTG)	ygiC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,179,239	0	C	T	16.2%	171.6 / 23.0	68	L55L (CTG→TTG)	ygiC	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (6/5);  total (36/32)
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.68e-01

CTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAC  >  W3110S.gb/3179169‑3179297                                                                       |                                                           cTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAAc                                                            >  1:3813477/1‑71 (MQ=255) cTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAAc                                                            >  1:841772/1‑71 (MQ=255)    gtgtGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                          >  1:2749436/1‑70 (MQ=255)         aaGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                       >  1:350313/1‑68 (MQ=255)         aaGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                       >  1:3365914/1‑68 (MQ=255)         aaGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                       >  1:3019535/1‑68 (MQ=255)         aaGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCa                                                     <  1:1079236/70‑1 (MQ=255)         aaGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCa                                                     <  1:1146035/70‑1 (MQ=255)         aaGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCa                                                     <  1:1410403/70‑1 (MQ=255)          aGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCaca                                                   <  1:2183682/71‑3 (MQ=255)          aGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                   <  1:2683596/71‑1 (MQ=255)            aTGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                   <  1:258487/69‑1 (MQ=255)            aTGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg                                                 >  1:3264055/1‑71 (MQ=255)             tGCTTACTACAAGTTGACTCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgt                                                >  1:4053584/1‑71 (MQ=255)             tGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgt                                                >  1:3205340/1‑71 (MQ=255)             tGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgt                                                >  1:2934708/1‑71 (MQ=255)             tGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgt                                                <  1:3215803/71‑1 (MQ=255)             tGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgt                                                >  1:716593/1‑71 (MQ=255)                    tacAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                   <  1:1876781/61‑1 (MQ=255)                      cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAg                                       >  1:1832936/1‑71 (MQ=255)                      cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg                                                 >  1:327060/1‑61 (MQ=255)                      cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                       >  1:2646171/1‑71 (MQ=255)                         gTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGc                                                         <  1:2975492/50‑1 (MQ=255)                         gTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGc                                                         <  1:126639/50‑1 (MQ=255)                         gTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                   <  1:3049432/56‑1 (MQ=255)                           tGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                  >  1:1213521/1‑71 (MQ=255)                            gACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAgtggtg                                  >  1:638314/1‑70 (MQ=255)                                cTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaa                             <  1:2523536/71‑1 (MQ=255)                                  cGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                              >  1:1982628/1‑52 (MQ=255)                                  cGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                              >  1:680378/1‑52 (MQ=255)                                  cGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                           >  1:1334381/1‑71 (MQ=255)                                   gCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                          <  1:1356714/71‑1 (MQ=255)                                       aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                              >  1:486768/1‑47 (MQ=255)                                       aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                              >  1:3956052/1‑47 (MQ=255)                                       aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                              >  1:2195393/1‑47 (MQ=255)                                        ggTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                     >  1:1348496/1‑71 (MQ=255)                                        ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                       >  1:1993703/1‑69 (MQ=255)                                         gTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                         <  1:309107/66‑1 (MQ=255)                                         gTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                         <  1:571462/66‑1 (MQ=255)                                         gTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                    <  1:1231759/55‑1 (MQ=255)                                           tGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGc                    >  1:104261/1‑69 (MQ=255)                                            gAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                 >  1:3225664/1‑71 (MQ=255)                                             aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                <  1:1004955/71‑1 (MQ=255)                                             aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                <  1:834523/71‑1 (MQ=255)                                             aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                <  1:2997577/71‑1 (MQ=255)                                             aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                <  1:773340/71‑1 (MQ=255)                                             aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGATGGAAAAAGTGATCGCCAGc                <  1:2381153/71‑1 (MQ=255)                                              aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                >  1:2791165/1‑70 (MQ=255)                                              aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                <  1:656025/70‑1 (MQ=255)                                              aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCg               <  1:3355675/71‑1 (MQ=255)                                                aGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGa              >  1:3447331/1‑70 (MQ=255)                                                aGCTGGAAGAAGACACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGa              >  1:937071/1‑70 (MQ=255)                                                 gCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg            <  1:1367231/71‑1 (MQ=255)                                                 gCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCg               <  1:539081/68‑1 (MQ=255)                                                 gCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt             >  1:3603272/1‑70 (MQ=255)                                                  cTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                         <  1:3810089/57‑1 (MQ=255)                                                     gaagaaGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt        <  1:1042212/71‑1 (MQ=255)                                                     gaagaaGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt        <  1:222371/71‑1 (MQ=255)                                                     gaagaaGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGa           >  1:4009180/1‑68 (MQ=255)                                                     gaagaaGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGa           >  1:3386418/1‑68 (MQ=255)                                                     gaagaaGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt        >  1:2059737/1‑71 (MQ=255)                                                     gaagaaGTCACAGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt        <  1:3200519/71‑1 (MQ=255)                                                      aagaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                    <  1:3984274/42‑1 (MQ=255)                                                       agaagTCACCGCCGAACTGCACCAGATGTGCCTGa                                          >  1:3138423/1‑35 (MQ=255)                                                        gaagTCACCGCCTAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga      >  1:1321288/1‑70 (MQ=255)                                                          agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga   >  1:3770243/1‑71 (MQ=255)                                                            tCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGa           <  1:4079062/61‑1 (MQ=255)                                                            tCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc  >  1:1401748/1‑70 (MQ=255)                                                              aCCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                      >  1:1658741/1‑48 (MQ=255)                                                                      aCTGCACCAGATGTGCCTGAAAGTGGTGGAATAAGTGATCGCCAGCg               <  1:3964046/47‑1 (MQ=255)                                                                       |                                                           CTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAC  >  W3110S.gb/3179169‑3179297

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