Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,547,129	G→A	M137I (ATG→ATA)	fimE →	tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,547,129	0	G	A	100.0%	161.4 / NA	49	M137I (ATG→ATA)	fimE	tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/30);  total (19/30)

TATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGG  >  W3110S.gb/4547067‑4547190                                                               |                                                              tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGATCCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:1969264/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:407645/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:2141657/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:221446/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:2511602/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:2761688/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:511719/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:2853082/65‑1 (MQ=255) tATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:1124098/65‑1 (MQ=255)  aTCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATAtt                                                             <  1:2549357/64‑1 (MQ=255)             cgcgATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTg                                                 <  1:146056/65‑1 (MQ=255)              gcgATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATATTAAg                                                          >  1:39333/1‑55 (MQ=255)              gcgATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGc                                                >  1:1981209/1‑65 (MQ=255)                      ggTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGa                                        >  1:2792877/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCAt                                                      >  1:1081917/1‑41 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATg                                                     >  1:2100000/1‑42 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGa                               >  1:562460/1‑64 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:2985583/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:3024639/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:571358/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:691811/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:2100877/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:823442/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:999094/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:951301/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:1223899/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:12637/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:1633888/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:1894732/1‑65 (MQ=255)                                cTGGAACCGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAg                              >  1:1869365/1‑65 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:706755/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:1867922/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:1114046/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:46170/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:318989/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:3053240/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:115355/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:2951085/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:2902944/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:287194/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:2034767/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:514675/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:2685290/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:1427857/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:909316/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:884389/65‑1 (MQ=255)                                       cGTAACGCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:613599/65‑1 (MQ=255)                                             gCAGACTCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:2625310/59‑1 (MQ=255)                                                   tCATCCTCATATATTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                       <  1:141464/53‑1 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAAt       >  1:1768217/1‑60 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:941527/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:2633952/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:488389/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:828478/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:1917964/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:1108864/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:1169852/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:1364974/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:14619/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:1781197/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:1837560/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:1914432/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:431090/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:2013831/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:2306870/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:2403136/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:2581844/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:270992/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:2835489/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:3056575/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAgg  >  1:35928/1‑65 (MQ=255)                                                            tataTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAg   >  1:787798/1‑64 (MQ=255)                                                               |                                                              TATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGG  >  W3110S.gb/4547067‑4547190

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