Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,546,601	C→T	I502I (ATC→ATT)	fimD →	type I fimbriae usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,546,601	0	C	T	100.0%	33.2 / NA	14	I502I (ATC→ATT)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/9);  total (5/9)

GGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAG  >  NC_000913/4546371‑4546784                                                                                                                                                                                                                                       |                                                                                                                                                                                        ggTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTAt                                                                                                                                                                     >  1:467417/1‑251 (MQ=255)                                         gggATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATGCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTAGTCAGGTTAAGCCGAAATTCTCCGACTATTACAACCTCGCt                                                                                                                             <  1:184683/250‑1 (MQ=255)                                                       ncTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAac                                                                                                                                                                                    >  1:54846/2‑182 (MQ=255)                                                                                                                           gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGc                                                                                     <  1:154378/209‑1 (MQ=255)                                                                                                                               aaTGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGAAATTCTAATTTCGCTGATACAACATCCAGTCGAATCAATGGCTACCACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACAGAATATTACAACCTCGCTTATAACAACCGCGGGAACTTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAAc                                       >  1:217657/1‑251 (MQ=255)                                                                                                                                   aaTCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGCCCGGCGGATATTTTAATTTCGCTGATACAACATCCAGTCAAATGAATGGCTACAACATTGAAACACAGGACGGATTTATTCAGTTTAAGCCGAAATTCCCCGACTATTACACCCTCGCTTATAACAAACGCgg                                                                                                               <  1:193294/175‑1 (MQ=255)                                                                                                                                                                  ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAGTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgg                                <  1:348705/223‑1 (MQ=255)                                                                                                                                                                    gTTATTTGACCAGCGGATATTTTAATTTCGCTTATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGATTTATTCAGGTTAAGCCGAAATTCTCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCGTCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATGTGAGTGGTAGCCAGCAACCTTATTGGGGAACGAGTAATGTCGATGATCAATTCCAg  <  1:87551/251‑1 (MQ=255)                                                                                                                                                                     ttATTCGACCAGCGGATATTTTAATTTCGCTGATACATCATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATGCACCGACTATTGCAACCGCGCTTATAACAAACGCGGGGAATTACAACTCGCCGTTTCGCGGCAACTCTCTCGCTCATTAACACTGTATTTGAGTGTTAGCCATCAAACTTATTTGGGAACGAGTAATGTCGATGAGCAATTCCAg  <  1:335573/250‑1 (MQ=255)                                                                                                                                                                                ncGGATANTTTAATTTCGCTGATACAANATACAGTCGNATGAANNNCNACAANATTGAAac                                                                                                                                                                                    <  2:54846/60‑1 (MQ=255)                                                                                                                                                                                                  cTGATACAACATACAGTCGAATNAATNGCTANAANATTGAAACACAGGACNGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCNNTNNNANCAAAcgcg                                                                                                                <  2:348961/111‑1 (MQ=255)                                                                                                                                                                                                  cTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                                                                                                                >  1:348961/1‑111 (MQ=255)                                                                                                                                                                                                       acaacaTTCAGTCGAATGAATGGCTACAACATTGAAACACAGGACTGAGTTATTCACGTTGAGCCGAAAGTCACCGACTATTACAACCTCGCGTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAg                                                                                      <  1:88152/132‑1 (MQ=255)                                                                                                                                                                                                        caacaTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGATATTCAGGTTAAGCCGACATTCACCGACTATTACAACCTCGCTTATAACAACCGCGGGCCATTACAACTCACCGTTACTc                                                                                        >  1:66661/1‑129 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                        GGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAG  >  NC_000913/4546371‑4546784

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