Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,213,335	A→G	34.0%	N301S (AAT→AGT)	yehQ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,213,335	0	A	G	34.0%	48.8 / 42.8	47	N301S (AAT→AGT)	yehQ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base G (8/8);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTGGCGCACG  >  W3110S.gb/2213273‑2213399                                                               |                                                                 cgcCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGcac                                                          <  1:2828044/71‑1 (MQ=255)   ccagccaCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACaa                                                        <  1:440696/71‑1 (MQ=255)   ccagccaCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACaa                                                        <  1:772551/71‑1 (MQ=255)     agccaCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAaga                                                      <  1:1953767/71‑1 (MQ=255)     agccaCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAaga                                                      <  1:1986222/71‑1 (MQ=255)      gccaCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAagag                                                     <  1:1524869/71‑1 (MQ=255)      gccaCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAagag                                                     <  1:458016/71‑1 (MQ=255)        caCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTaaa                                                                  >  1:1596735/1‑56 (MQ=255)        caCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTaaa                                                                  >  1:287977/1‑56 (MQ=255)        caCTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTaaa                                                                  >  1:1509942/1‑56 (MQ=255)           tataATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAaga                                                      <  1:3158903/65‑1 (MQ=255)           tataATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAagag                                                     <  1:3290667/66‑1 (MQ=255)            ataATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAAGAGATGGcc                                               >  1:45396/1‑71 (MQ=255)            ataATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGcc                                               >  1:678037/1‑71 (MQ=255)            ataATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGcc                                               >  1:3276823/1‑71 (MQ=255)             taATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAAGAGATGGccc                                              >  1:655537/1‑71 (MQ=255)                  ccGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTGGATTAAGTTGCGCACAAGAGATGGccc                                              >  1:363668/1‑66 (MQ=255)                  ccGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAAGAGATGGccc                                              >  1:3655756/1‑66 (MQ=255)                   cGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGcac                                                          <  1:909424/53‑1 (MQ=255)                        aaTGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAg                                    <  1:2979014/70‑1 (MQ=255)                            cTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTg                               >  1:3012111/1‑71 (MQ=255)                             ttACATCAGCTTGCGGCATTAAACAGTCGATTaaa                                                                  >  1:2523548/1‑35 (MQ=255)                               aCATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAGTa                                  >  1:2503018/1‑65 (MQ=255)                               aCATCAGCTTGCGGCATTAAACAGTCGATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAGTa                                  >  1:3279901/1‑65 (MQ=255)                                 aTCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAaga                                                      <  1:2559674/43‑1 (MQ=255)                                 aTCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAaga                                                      <  1:3473190/43‑1 (MQ=255)                                 aTCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAaga                                                      <  1:751383/43‑1 (MQ=255)                                  tCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATgg                                                 <  1:30831/47‑1 (MQ=255)                                    aGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGcc                                               >  1:1724069/1‑47 (MQ=255)                                    aGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGcc                                               >  1:834612/1‑47 (MQ=255)                                         gCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTtcctc                   <  1:196102/70‑1 (MQ=255)                                          cGGCATTAAACAGTCGATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTtcctcc                  <  1:1878839/70‑1 (MQ=255)                                          cGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTtcctc                   <  1:2189447/69‑1 (MQ=255)                                              aTTAAACAGTCCATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGt                        <  1:3692033/60‑1 (MQ=255)                                               ttAAACAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGt                        <  1:3207107/59‑1 (MQ=255)                                                 aaaCAGTCGATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCg           >  1:184413/1‑70 (MQ=255)                                                 aaaCAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTgg                              >  1:2825135/1‑51 (MQ=255)                                                 aaaCAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGt          >  1:185906/1‑71 (MQ=255)                                                  aaCAGTCGATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTg         <  1:1936141/71‑1 (MQ=255)                                                  aaCAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTAt                                 <  1:46227/47‑1 (MQ=255)                                                    cAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTg                               >  1:655925/1‑47 (MQ=255)                                                      gTCGATTAAGTTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAg                         >  1:1107570/1‑51 (MQ=255)                                                      gTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTtcctcc                  <  1:1175745/58‑1 (MQ=255)                                                      gTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTtcctcc                  <  1:1981207/58‑1 (MQ=255)                                                         gATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTa                                  >  1:74079/1‑39 (MQ=255)                                                         gATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTGGCGCACg  >  1:405865/1‑71 (MQ=255)                                                         gATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTGGCGCACg  >  1:182943/1‑71 (MQ=255)                                                         gATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTGGCGCACg  >  1:3692173/1‑71 (MQ=255)                                                         gATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTGGCGCACg  >  1:3399752/1‑71 (MQ=255)                                                         gATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTGGCGCACg  >  1:2786940/1‑71 (MQ=255)                                                         gATTAAATTGCGCACAAGAGATGGCACGGCGCGACAGTa                                  >  1:365711/1‑39 (MQ=255)                                                             aaaTTGCGCACAAGAGATGGCCCGGCGCGACAGTATTgg                              >  1:683210/1‑39 (MQ=255)                                                               |                                                                 CGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGATGGCCCGGCGCGACAGTATTGGTGAAGTTCCTCCTGTGCCGTGGCGCACG  >  W3110S.gb/2213273‑2213399

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