Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,544,435	A→G	intergenic (‑252/‑1204)	yjhA ← / → fimB	N‑acetylnuraminic acid outer membrane channel protein/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,544,435	0	A	G	100.0%	274.7 / NA	76	intergenic (‑252/‑1204)	yjhA/fimB	N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (44/32);  total (44/32)

CTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGATCTC  >  W3110S.gb/4544384‑4544492                                                    |                                                          cTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACTTa                                                  >  1:1519415/1‑61 (MQ=255) cTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACTTa                                                  >  1:2307336/1‑61 (MQ=255)            aTCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCAttt                                      <  1:1632860/62‑1 (MQ=255)            aTCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCAttt                                      <  1:711724/62‑1 (MQ=255)            aTCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCAttt                                      <  1:2095235/62‑1 (MQ=255)            aTCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCAttt                                      <  1:139670/62‑1 (MQ=255)             tctcTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:255548/1‑47 (MQ=255)             tctcTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:1903647/1‑47 (MQ=255)             tctcTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:3066074/1‑47 (MQ=255)             tctcTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:2404490/1‑47 (MQ=255)             tctcTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:2962000/1‑47 (MQ=255)             tctcTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:252055/1‑47 (MQ=255)              ctctTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:1573986/1‑46 (MQ=255)              ctctTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:347945/1‑46 (MQ=255)                 ttctGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACt                                                    >  1:1773414/1‑43 (MQ=255)                 ttctGATGTTTAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTa                                    >  1:757157/1‑59 (MQ=255)                         tttAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCAt                         <  1:3015652/62‑1 (MQ=255)                         tttAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCAt                         <  1:53253/62‑1 (MQ=255)                         tttAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCAt                         <  1:3109879/62‑1 (MQ=255)                           tAATTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       <  1:2411981/62‑1 (MQ=255)                             aTTGATTTGAATGTTGGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:2773585/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:2597874/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:1665001/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:1657463/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:72734/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:2574763/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:98081/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:630479/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:2475398/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:3025944/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:1454976/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:1345657/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:1112786/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGCTGTATAACTTACTCAATCCATTTTACCCAGAGTCATaa                       >  1:1472497/1‑60 (MQ=255)                             aTTGATTTGAATGTTCGTAAGATGTATCACTTACTCAATCCATTTTACCCAGAGTCATaaa                      >  1:2036932/1‑61 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:1053294/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:770618/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:2464383/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:1050905/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:916337/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:2438646/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:3053190/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:2186698/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:2155759/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:505280/1‑62 (MQ=255)                                  tttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                >  1:362478/1‑62 (MQ=255)                                   ttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAAt               <  1:417508/62‑1 (MQ=255)                                   ttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAAt               <  1:576834/62‑1 (MQ=255)                                   ttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAAt               <  1:715633/62‑1 (MQ=255)                                   ttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAAt               <  1:2379203/62‑1 (MQ=255)                                   ttGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAAt               <  1:1722526/62‑1 (MQ=255)                                    tGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATg              >  1:606635/1‑62 (MQ=255)                                    tGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATg              >  1:1087504/1‑62 (MQ=255)                                    tGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATg              >  1:1822144/1‑62 (MQ=255)                                    tGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATg              >  1:3106206/1‑62 (MQ=255)                                    tGAATGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATg              >  1:2101760/1‑62 (MQ=255)                                      aaTGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGt             >  1:1572383/1‑61 (MQ=255)                                       aagTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTc            <  1:302348/59‑1 (MQ=255)                                       aTGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTc            <  1:2271937/61‑1 (MQ=255)                                       aTGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTc            <  1:2275267/61‑1 (MQ=255)                                       aTGTTCGTAAGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTc            <  1:2650861/61‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCTTTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:448096/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:85860/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:86323/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:244773/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:2909649/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:2521807/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:2460523/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1005224/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1983521/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1933173/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1843540/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1739016/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1546616/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1054681/62‑1 (MQ=255)                                                aGCTGTATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc  <  1:1039395/62‑1 (MQ=255)                                                    |                                                          CTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGATCTC  >  W3110S.gb/4544384‑4544492

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