Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,910,338	C→T	8.3%	*48* (TGA→TAA)	yobG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,910,338	0	C	T	8.3%	145.8 / 2.9	48	*48* (TGA→TAA)	yobG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base T (2/2);  total (27/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.85e-01

CACCTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCT  >  W3110S.gb/1910268‑1910392                                                                       |                                                       caccTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATc                                                        <  1:2646661/71‑1 (MQ=255) caccTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATc                                                        <  1:1516310/71‑1 (MQ=255) caccTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATc                                                        <  1:1944764/71‑1 (MQ=255)      ctGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcaccac                                                   >  1:69839/1‑71 (MQ=255)      ctGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcaccac                                                   >  1:684058/1‑71 (MQ=255)      ctGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcaccac                                                   >  1:578831/1‑71 (MQ=255)                 cgcTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTg                                        >  1:733668/1‑71 (MQ=255)                  gcTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTgg                                       >  1:644541/1‑71 (MQ=255)                     aTTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGtt                                     >  1:728226/1‑70 (MQ=255)                      ttCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATcac                                                      <  1:2105281/52‑1 (MQ=255)                      ttCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:1570573/1‑71 (MQ=255)                        cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                 >  1:2997901/1‑71 (MQ=255)                        cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                 >  1:2564313/1‑71 (MQ=255)                        cTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATg                                 >  1:1757896/1‑71 (MQ=255)                             aCTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGcaca                            <  1:737596/71‑1 (MQ=255)                               tgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                          >  1:46304/1‑71 (MQ=255)                               tgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                          >  1:4059665/1‑71 (MQ=255)                               tgctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACaaa                          >  1:2579260/1‑71 (MQ=255)                               agctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATCGCACaaa                          >  1:2583904/2‑71 (MQ=255)                                gctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAAt                         >  1:5239/1‑71 (MQ=255)                                 ctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATaaa                                           >  1:347112/1‑52 (MQ=255)                                 ctgGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATaaa                                           >  1:948052/1‑52 (MQ=255)                                           agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTa                                    >  1:925949/1‑49 (MQ=255)                                           agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTa                                    >  1:115167/1‑49 (MQ=255)                                           agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTa                                    >  1:2188871/1‑49 (MQ=255)                                           agaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa              <  1:859078/71‑1 (MQ=255)                                            gaaAATCTAGTGCTGAAAAAATGATATCACCACGGGATa                                             <  1:3030000/39‑1 (MQ=255)                                             aaaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaa              <  1:3758526/69‑1 (MQ=255)                                              aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATa                                             <  1:1616271/37‑1 (MQ=255)                                              aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:3953363/1‑47 (MQ=255)                                              aaaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:3108795/1‑47 (MQ=255)                                               aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaa               <  1:2713563/66‑1 (MQ=255)                                               aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:3706595/71‑1 (MQ=255)                                               aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:2976381/71‑1 (MQ=255)                                               aaTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          >  1:779246/1‑71 (MQ=255)                                                aTCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:1564016/70‑1 (MQ=255)                                                   tAGTGCTGAAAAAATGATATCACCACGGGATAAACTgg                                       >  1:3388313/1‑38 (MQ=255)                                                     gTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTaa                                   >  1:1491659/1‑40 (MQ=255)                                                     gTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTAc       <  1:2083158/68‑1 (MQ=255)                                                      tGCTGAAAAAATGATATTACCACGGGATAAACTGGtt                                     <  1:97111/37‑1 (MQ=255)                                                      tGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTa        <  1:2520392/66‑1 (MQ=255)                                                        cTGAAAAAATGATATTACCACGGGATAAACTGGTTaa                                   >  1:3187843/1‑37 (MQ=255)                                                        cTGAAAAAATGATATTACCACGGGATAAACTGGTTaa                                   >  1:125420/1‑37 (MQ=255)                                                        cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCt  <  1:1896717/70‑1 (MQ=255)                                                        cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCt  >  1:577635/1‑70 (MQ=255)                                                        cTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCt  <  1:746556/70‑1 (MQ=255)                                                                aTGATATTACCACGGGATAAACTGGTTAATGGCACaaa                          <  1:2404575/38‑1 (MQ=255)                                                                aTGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTAc       >  1:1164366/1‑57 (MQ=255)                                                                aTGATAGCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTAc       >  1:522408/1‑57 (MQ=255)                                                                  gATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa             <  1:3152796/49‑1 (MQ=255)                                                                  gATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGaaaaa             <  1:2241729/49‑1 (MQ=255)                                                                      tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:607262/48‑1 (MQ=255)                                                                      tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:882785/48‑1 (MQ=255)                                                                      tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:887647/48‑1 (MQ=255)                                                                      tCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTg          <  1:2681780/48‑1 (MQ=255)                                                                       |                                                       CACCTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCT  >  W3110S.gb/1910268‑1910392

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