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

GTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAC  >  NC_000913/4546372‑4546847                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        gTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGACTGAATGGCTACAACATTGAACCACAGGACGGAGTTAtt                                                                                                                                                                                                                                    >  2:194307/1‑251 (MQ=255)          aaaaCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTa                                                                                                                                                                                                                             >  1:260413/1‑249 (MQ=255)           aaaCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGc                                                                                                                                                                                                                          >  1:130252/1‑251 (MQ=255)                        cTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                                                                                                                                                                              >  1:6178/1‑250 (MQ=255)                             cgcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                        <  1:3132/251‑1 (MQ=255)                               ctctGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                                                                                                                                                                                      >  1:303102/1‑251 (MQ=255)                                             atGACGCAGGCTAATTCCACACGTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCGCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                                                                                                                                                                                        <  2:149008/251‑1 (MQ=255)                                                                         aTGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCCCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTAGACAAGCGGATATTTTAATTTCGCTGCTACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACAGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTA‑AACATCACCGtt                                                                                                                                                           >  2:275220/1‑251 (MQ=255)                                                                             cAGTCAGCATGACGGACAATCGGTGCGTTTCCTGTATTACAAATCACTCAATGAATAAGGCACGAATATTCAGTTCGAGGGTTACCGTTATTCGACCAGCGGATCTTTTAATTTCGCTGATACAACATCCAGTCGAATGAATGGCTACAACATTGAAACACAGGAGGGAAATATTCAGGTTAAGCCGAAATTCACCGACTATTACCACCTAGCTTAAAACAAACGCAGGAAATTACAAC‑TCACCGTTACt                                                                                                                                                        >  1:293257/1‑250 (MQ=255)                                                                                     aTGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACAGACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGTGAAATTACAAC‑TCACCGTTACTCAGCAACTc                                                                                                                                               <  1:38085/251‑1 (MQ=255)                                                                                                                          gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTGACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                                                                                                          <  1:267127/251‑1 (MQ=255)                                                                                                                               aTGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATGAACACTGTATTTGAGTGGTAGCCATCAAACt                                                                                                     <  1:265465/251‑1 (MQ=255)                                                                                                                                     cAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGt                                                                                                                            <  2:51832/222‑1 (MQ=255)                                                                                                                                     cAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGt                                                                                                                            >  1:51832/1‑222 (MQ=255)                                                                                                                                                atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAAGTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCGTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACGGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                                                                                    <  1:140783/251‑1 (MQ=255)                                                                                                                                                              ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACACCATACAGTCGAATGCATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTACGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGTGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAAt                                                                      >  2:203380/1‑251 (MQ=255)                                                                                                                                                                            ccAGCGGATATTTTACTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTACGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAtt                                                         >  1:268033/1‑250 (MQ=255)                                                                                                                                                                                 ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCGCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATAc                                                   <  2:303102/251‑1 (MQ=255)                                                                                                                                                                                         tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAAGACTGCGTTc                                            <  2:268033/250‑1 (MQ=255)                                                                                                                                                                                                       caacaTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTgg                             >  2:123612/1‑251 (MQ=255)                                                                                                                                                                                                                                  acaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCGCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAc  <  1:302222/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        GTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAC‑TCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAC  >  NC_000913/4546372‑4546847

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