Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,639,730	T→C	K222R (AAA→AGA)	xerC ←	site‑specific tyrosine recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,639,730	0	T	C	100.0%	133.3 / NA	36	K222R (AAA→AGA)	xerC	site‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/2);  total (34/2)

ATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGT  >  W3110S.gb/3639681‑3639787                                                  |                                                          aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTAc                                               <  1:2173870/62‑1 (MQ=255) aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTAc                                               <  1:3141092/62‑1 (MQ=255)                 aGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCtt                               >  1:322295/1‑62 (MQ=255)                    ccTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGga                                        >  1:967612/1‑50 (MQ=255)                    ccTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGga                                        >  1:159458/1‑50 (MQ=255)                    ccTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTg                              >  1:3487431/1‑60 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:395976/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:387426/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:1177548/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:3297307/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:4086025/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:320006/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:4430508/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:4535591/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:279937/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:677524/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:768370/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:2432341/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg                                          >  1:1843476/1‑38 (MQ=255)                              tGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgc                                           >  1:2438481/1‑37 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGt                        >  1:3567063/1‑42 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc    >  1:4017828/1‑62 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc    >  1:4133888/1‑62 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc    >  1:4554426/1‑62 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc    >  1:48729/1‑62 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc    >  1:2678087/1‑62 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc    >  1:1400161/1‑62 (MQ=255)                                            aaGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgt     >  1:3212362/1‑61 (MQ=255)                                            aaGCGTCCCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc    >  1:3069797/1‑62 (MQ=255)                                              gCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt  >  1:4157470/1‑62 (MQ=255)                                              gCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt  >  1:4203471/1‑62 (MQ=255)                                              gCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt  >  1:2828594/1‑62 (MQ=255)                                              gCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt  >  1:611311/1‑62 (MQ=255)                                              gCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt  >  1:661869/1‑62 (MQ=255)                                              gCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt  >  1:2718474/1‑62 (MQ=255)                                              gCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt  >  1:836260/1‑62 (MQ=255)                                                  |                                                          ATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGT  >  W3110S.gb/3639681‑3639787

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