Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,591,162	C→A	28.4%	V112V (GTG→GTT)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,591,162	0	C	A	28.4%	88.9 / 20.6	42	V112V (GTG→GTT)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand):  ref base C (20/10);  new base A (4/8);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 8.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGC  >  W3110S.gb/4591093‑4591227                                                                      |                                                                  ccAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCa                                                                  >  1:488740/1‑71 (MQ=255) ccAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCa                                                                  >  1:2503812/1‑71 (MQ=255) ccAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCa                                                                  >  1:3437582/1‑71 (MQ=255)   aGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCAc                                                                 >  1:2270758/1‑70 (MQ=255)     aCAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCAcc                                                                <  1:1696464/69‑1 (MQ=255)              cGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                     <  1:3025217/71‑1 (MQ=255)              cGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                     <  1:1515604/71‑1 (MQ=255)              cGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                     <  1:4080638/71‑1 (MQ=255)               ggAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                     <  1:2997018/70‑1 (MQ=255)                  aCCCGAGTCGC‑‑GCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATgcg                                               <  1:2371105/71‑1 (MQ=255)                    ccGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATgcg                                               <  1:938359/71‑1 (MQ=255)                     cGAGTCGCAGGCACATCTGGGCATCGTTGAGATGGTTATGATATTCGTGCACCGCCTGGTTACCAATgcgg                                              >  1:3654361/1‑71 (MQ=255)                      gAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACAAATgcg                                               >  1:3298995/1‑69 (MQ=255)                           gCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAAt                                        <  1:1664533/71‑1 (MQ=255)                                 acaTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTtgtg                                   >  1:239063/1‑70 (MQ=255)                                 acaTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGa                                  <  1:429191/71‑1 (MQ=255)                                   aTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCTCCGCCTGGTTACCAATGCGGCGTAATTTGTCaaa                                <  1:1383360/71‑1 (MQ=255)                                    tCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGaaa                                <  1:3076129/70‑1 (MQ=255)                                    tCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGaaa                                >  1:1286594/1‑70 (MQ=255)                                      tGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCaa                                                   >  1:1376282/1‑49 (MQ=255)                                      tGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCaa                                                   >  1:3132352/1‑49 (MQ=255)                                        ggCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTTAAATa                              <  1:3795082/68‑1 (MQ=255)                                            tCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACa                            >  1:2119467/1‑66 (MQ=255)                                            tCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACagagag                       <  1:3185262/71‑1 (MQ=255)                                              gTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgc                                                >  1:2413182/1‑44 (MQ=255)                                              gTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGaaa                                >  1:2899106/1‑60 (MQ=255)                                              gTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGaaa                                >  1:1552763/1‑60 (MQ=255)                                              gTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGaaa                                >  1:239704/1‑60 (MQ=255)                                                tGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCaa                                                   >  1:2137066/1‑39 (MQ=255)                                                         tATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACagagag                       <  1:3323712/58‑1 (MQ=255)                                                         tATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACagagag                       <  1:3153344/58‑1 (MQ=255)                                                            gATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATa                              <  1:3081283/48‑1 (MQ=255)                                                             atatTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATcaaca       <  1:2253538/70‑1 (MQ=255)                                                              tatTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACaga                          >  1:1409983/1‑50 (MQ=255)                                                              tatTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACaaa     <  1:615062/71‑1 (MQ=255)                                                               atTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATcaa         >  1:3779138/1‑66 (MQ=255)                                                               atTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATcaa         >  1:522514/1‑66 (MQ=255)                                                                  cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGc  >  1:2966928/1‑70 (MQ=255)                                                                  cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGc  >  1:98815/1‑70 (MQ=255)                                                                  cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGc  >  1:3272313/1‑70 (MQ=255)                                                                  cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGc  >  1:2472591/1‑70 (MQ=255)                                                                  cGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGc  >  1:3484238/1‑70 (MQ=255)                                                                      |                                                                  CCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGC  >  W3110S.gb/4591093‑4591227

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