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%	48.8 / NA	19	I502I (ATC→ATT)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/8);  total (11/8)

CATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTG  >  NC_000913/4546385‑4546819                                                                                                                                                                                                                         |                                                                                                                                                                                                                           cATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGATAATTCCACCCTTCCCGATGACAGTCAGCATGAAGGACAATAGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACTGTTATTCGACCCGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                                                                                                                                                                           >  1:445586/1‑250 (MQ=255)                  ctctGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTGAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTaan                                                                                                                                                                                              <  1:80360/230‑2 (MQ=255)                                                         ccGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAacac                                                                                                                                                                                                                     >  1:236899/1‑168 (MQ=255)                                                                              gACAATCTGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCTGATATTTTAATTTCGCTGATACAACATACAGGCGAATGAAGGTCTACAACATTGAAACACAGGGCGGAGTGATTCAAGTTAAGCCGAAATTCACCGACGATTACAACCTCGCTTATAACAAACGCGGGAAATTACAtctc                                                                                                                                     <  1:451267/227‑1 (MQ=255)                                                                                               ctctATAACAACTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGGGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata                                                                                                                                                            >  1:149366/1‑187 (MQ=255)                                                                                                tctATAACAAATCGATCAATGAATCCGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCTAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTAACCGACTATTACAAGCTCGCTTATAACAAACGCGTGAAATTACAACTCACCGTTACTAAGAAACTCGGGCGCACATCAACCCTGTATTTg                                                                                           >  1:102096/1‑251 (MQ=255)                                                                                                             gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCggg                                                                                                                                                 >  1:169297/1‑184 (MQ=255)                                                                                                                          gnCACGACTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACa                                                                                                                                         >  1:349464/1‑179 (MQ=255)                                                                                                                              cGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACg                                                             <  1:342117/251‑1 (MQ=255)                                                                                                                                              gggTTACCGTTATTCGACCAGCGGATATTTTACTTTCGCTGATACAACATACAGTCGAATGCATGGCTACAACATTGACACACAGGACGGAGTTATTCAGGTTAAGCCGCAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCACCACTGTATTTGAGTGGTAGCAATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                                             >  1:321846/1‑251 (MQ=255)                                                                                                                                              gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTTAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACTCTGGAAATTACAACGCACCGTTACTCAGCAACCCGGACCCGCATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAg                                              <  1:202584/250‑1 (MQ=255)                                                                                                                                              gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATAGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCAGAAATTCCCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACAGTTACTCAGCAACTCGGGCGCACATCAACACTGAATTTGAGTGGGAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                                             >  1:286015/1‑251 (MQ=255)                                                                                                                                                                    ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATAc                       >  1:369167/1‑251 (MQ=255)                                                                                                                                                                    ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCGACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCGCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTc                                                                                                                   <  1:369278/159‑1 (MQ=255)                                                                                                                                                                          tnnAATTTCNNTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCTGAAATTCACCGACAATTAc                                                                                                                                                                         <  2:5330/99‑1 (MQ=255)                                                                                                                                                                                   gCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGatat        >  1:16581/1‑251 (MQ=255)                                                                                                                                                                                       ntcnaacnTNNANNCGANTGAATGGCTACAACATTGAAANNCAGGACGGANTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAANGCGGGAAATNNCANCTCACNGNNNCNNNGCNACTCGGGCGCACATCAACANTGNNTTTGAGNGGTAGCCATCAAACTTATTGGGGAAc                                                              <  2:102096/189‑1 (MQ=255)                                                                                                                                                                                         acaacaTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGATCGAGTAATGGCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACtn  <  1:89905/251‑2 (MQ=255)                                                                                                                                                                                                   gTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                                                    >  1:139948/1‑159 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                           CATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTG  >  NC_000913/4546385‑4546819

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