Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,586,637	G→A	P365L (CCG→CTG)	hsdM ←	DNA methylase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,586,637	0	G	A	100.0%	211.4 / NA	65	P365L (CCG→CTG)	hsdM	DNA methylase M
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (60/5);  total (60/5)

GCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCG  >  W3110S.gb/4586590‑4586693                                                |                                                         gCAGGTCATACACCCACACATCATCGGTAGAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:2786113/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGt                                             >  1:2284278/1‑61 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:3482840/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:3646514/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:2980164/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:2971599/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:2908450/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:3678325/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:3824342/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:4348815/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:2197347/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:4374301/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:1900201/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:1874119/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:4397687/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:575673/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:678759/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:1548974/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:1313470/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:1453126/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTc                                            >  1:1392674/1‑62 (MQ=255) gCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCAACCGTc                                            >  1:2292098/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:1737296/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:3298400/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:3596395/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:3016285/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:3646962/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:833924/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:4248936/1‑62 (MQ=255)   aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:4484323/1‑62 (MQ=255)       cATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTg                                      <  1:1077817/62‑1 (MQ=255)       cATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTg                                      <  1:3539190/62‑1 (MQ=255)       cATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTg                                      <  1:1979685/62‑1 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCAcc                                               >  1:1316427/1‑48 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTaa                                  >  1:773542/1‑61 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTaa                                  >  1:3834081/1‑61 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTaa                                  >  1:4020863/1‑61 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTaa                                  >  1:951183/1‑61 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTaa                                  >  1:1117594/1‑61 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTaa                                  >  1:1325780/1‑61 (MQ=255)            acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTa                                   >  1:263168/1‑60 (MQ=255)             cccACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAg                                >  1:931355/1‑62 (MQ=255)              ccACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAg                                >  1:2776812/1‑61 (MQ=255)                  acaTCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                           >  1:3053832/1‑62 (MQ=255)                  acaTCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                           >  1:2052648/1‑62 (MQ=255)                  acaTCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                           >  1:4104922/1‑62 (MQ=255)                  acaTCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                           >  1:4206350/1‑62 (MQ=255)                  acaTCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                           >  1:4242817/1‑62 (MQ=255)                      catcGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:4252747/1‑43 (MQ=255)                      catcGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:1324872/1‑43 (MQ=255)                      catcGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:4330582/1‑43 (MQ=255)                      catcGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTccc                                          >  1:4259129/1‑43 (MQ=255)                            tACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGt                  >  1:2794284/1‑61 (MQ=255)                            tACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTc                 >  1:2696073/1‑62 (MQ=255)                            tACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTc                 <  1:1692331/62‑1 (MQ=255)                            tACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTc                 >  1:1635606/1‑62 (MQ=255)                            tACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTc                 >  1:317397/1‑62 (MQ=255)                                 ttCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCAc            >  1:2963032/1‑62 (MQ=255)                                         cTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGaa                              >  1:1651899/1‑36 (MQ=255)                                         cTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAg    >  1:3325769/1‑62 (MQ=255)                                         cTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAg    >  1:2376368/1‑62 (MQ=255)                                         cTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAg    >  1:1859608/1‑62 (MQ=255)                                          tGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAg    >  1:2620117/1‑61 (MQ=255)                                          tGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAg    >  1:3277472/1‑61 (MQ=255)                                           gATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCg  <  1:4342965/62‑1 (MQ=255)                                                |                                                         GCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCG  >  W3110S.gb/4586590‑4586693

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