Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,549,754	G→A	intergenic (+43/‑24)	fimC → / → fimD	chaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,549,754	0	G	A	100.0%	125.6 / NA	42	intergenic (+43/‑24)	fimC/fimD	chaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/24);  total (18/24)

AAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACA  >  W3110S.gb/4549688‑4549818                                                                   |                                                                 aaaaTGACGGGCGTAATGGAATAACGCATGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTg                                                              <  1:1797866/71‑1 (MQ=255) aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTg                                                              <  1:3509680/71‑1 (MQ=255)    aTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGAtt                                                           >  1:1835546/1‑71 (MQ=255)                aTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCggag                                                >  1:3461297/1‑70 (MQ=255)                     aTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaataa                                          <  1:2620056/71‑1 (MQ=255)                      tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaataat                                         <  1:4498845/71‑1 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       >  1:675103/1‑71 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       >  1:4281095/1‑71 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       <  1:4149667/71‑1 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       >  1:3867432/1‑71 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       >  1:3536298/1‑71 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       <  1:2379623/71‑1 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       <  1:1764058/71‑1 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       >  1:1289825/1‑71 (MQ=255)                        aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       <  1:1698149/71‑1 (MQ=255)                            aGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCat                                    <  1:3406252/70‑1 (MQ=255)                             gggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaa                                             <  1:3082020/60‑1 (MQ=255)                              ggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                            <  1:1169585/60‑1 (MQ=255)                               gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATc                                 >  1:1294179/1‑70 (MQ=255)                                 gAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCt                                >  1:487216/1‑69 (MQ=255)                                 gAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTg                               >  1:173827/1‑70 (MQ=255)                                  aaTTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTaa                                                       >  1:2862366/1‑45 (MQ=255)                                  aaTTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaataa                                          >  1:4164497/1‑58 (MQ=255)                                        tCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                        <  1:1477209/70‑1 (MQ=255)                                         cGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAg                       <  1:1763788/70‑1 (MQ=255)                                          gCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGAc                     <  1:1421079/71‑1 (MQ=255)                                             tGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                            <  1:495268/45‑1 (MQ=255)                                             tGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                            <  1:2048137/45‑1 (MQ=255)                                              gAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                       <  1:4058546/49‑1 (MQ=255)                                               aaTAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaataa                                          >  1:1587686/1‑45 (MQ=255)                                                   aaaaGAATTGACTGCCAGGTGATTTTAAGCCggagg                                               >  1:715365/1‑36 (MQ=255)                                                   aaaaGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGc            >  1:809309/1‑71 (MQ=255)                                                    aaaGAATTGACTGCCAGGTGATTTTAAGCCGGAGGaat                                            <  1:1990810/38‑1 (MQ=255)                                                      aGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAg                       <  1:943948/57‑1 (MQ=255)                                                          ttGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAaca      >  1:2508133/1‑70 (MQ=255)                                                            gACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaaa        >  1:981784/1‑66 (MQ=255)                                                            gACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca    >  1:1114676/1‑70 (MQ=255)                                                             aCTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTTAATTTAAGACTTTACCAGCGAAacacac   <  1:884894/70‑1 (MQ=255)                                                             aCTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacacac   <  1:4113999/70‑1 (MQ=255)                                                             aCTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacacac   <  1:874971/70‑1 (MQ=255)                                                             aCTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacacac   <  1:2338686/70‑1 (MQ=255)                                                             aCTGCCAGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAAcacaca  <  1:3355478/71‑1 (MQ=255)                                                                   |                                                                 AAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACA  >  W3110S.gb/4549688‑4549818

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