Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,087	G→T	G94W (GGG→TGG)	fimH →	type 1 fimbriae D‑mannose specific adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,087	0	G	T	100.0%	98.6 / NA	26	G94W (GGG→TGG)	fimH	type 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/16);  total (10/16)

ATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTC  >  NC_000913/4549013‑4549160                                                                           |                                                                          aTCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGAcc                                                                      <  1:863914/80‑1 (MQ=255)    cGGAAACCATTACAGACTATGTCACACTGCAGCGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTaa                                                                  <  2:863916/81‑1 (MQ=255)    cGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTaa                                                                  <  1:863917/81‑1 (MQ=255)    cGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTaa                                                                  <  2:863915/81‑1 (MQ=255)        aaCCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAata                                                              <  1:863918/81‑1 (MQ=255)          ccATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAata                                                              >  1:863957/1‑79 (MQ=255)            aTTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGt                                                          <  2:863919/81‑1 (MQ=255)                cAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGc                                                       <  1:863921/80‑1 (MQ=255)                cAGACTATGTCACACTGCAACGAGACTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCa                                                      <  2:863920/81‑1 (MQ=255)                 aGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAg                                                     <  1:863922/81‑1 (MQ=255)                  gACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGt                                                    <  2:863923/81‑1 (MQ=255)                             acTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCAtt                                         >  1:863948/1‑81 (MQ=255)                              cTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCAttt                                        <  1:863924/81‑1 (MQ=255)                               tGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTc                                       >  1:863953/1‑81 (MQ=255)                                 cAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTcc                                      <  2:863925/80‑1 (MQ=255)                                  aaCGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTa                                    <  1:863926/81‑1 (MQ=255)                                    cGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc                                  <  1:863927/81‑1 (MQ=255)                                     gAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacca                                 >  1:863963/1‑81 (MQ=255)                                     gAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacca                                 >  1:863965/1‑81 (MQ=255)                                           cGGCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGc                            >  2:863968/1‑80 (MQ=255)                                            ggCTTATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCg                           <  2:863928/80‑1 (MQ=255)                                                tATGGCGGCGTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcg                      >  1:863959/1‑81 (MQ=255)                                                    gcggcgTGTTATCTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcgccgc                  <  1:863929/81‑1 (MQ=255)                                                                cTAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtat       >  2:863964/1‑80 (MQ=255)                                                                 tAATTTTTCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATaa     >  1:863954/1‑81 (MQ=255)                                                                       ttCCTGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTc  >  2:863960/1‑78 (MQ=255)                                                                           |                                                                          ATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTC  >  NC_000913/4549013‑4549160

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