Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,175,236	C→A	Q102H (CAG→CAT)	mfd ←	transcription‑repair coupling factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,175,236	0	C	A	100.0%	112.0 / NA	36	Q102H (CAG→CAT)	mfd	transcription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/25);  total (11/25)

GAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAA  >  W3110S.gb/1175217‑1175301                    |                                                                  gAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                      <  1:2527018/49‑1 (MQ=255) gAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                      <  1:3453214/49‑1 (MQ=255) gAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                      <  1:3400794/49‑1 (MQ=255) gAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                      <  1:1945398/49‑1 (MQ=255) gAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                      <  1:2836764/49‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:3969864/56‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:2813556/56‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:2215471/56‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:402048/56‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:1867313/56‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:1722368/56‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:3264622/56‑1 (MQ=255)         aGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:385814/56‑1 (MQ=255)         aGTACGCCACGATCCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:4029145/56‑1 (MQ=255)         aGTACGCCACGATCCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:601462/56‑1 (MQ=255)          gTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                       <  1:319913/55‑1 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:676093/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:2355364/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:542555/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:180389/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:1754099/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:1703813/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:164856/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:1570460/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATaa                     >  1:150101/1‑56 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGc       >  1:2429811/1‑70 (MQ=255)           tACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGc       >  1:3548171/1‑70 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:1312194/71‑1 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:3503743/71‑1 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:2900902/71‑1 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:2826628/71‑1 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:2350464/71‑1 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:1954967/71‑1 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:1314173/71‑1 (MQ=255)               ccACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:1028397/71‑1 (MQ=255)                cACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaaa  <  1:806228/70‑1 (MQ=255)                    |                                                                  GAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAA  >  W3110S.gb/1175217‑1175301

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