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

CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAA  >  NC_000913/4546370‑4546845                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     cGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCCCACGTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTa                                                                                                                                                                                                                                   <  2:93400/251‑1 (MQ=255)  ggTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGTCTACAACATTGAAAAACAGGACGGAGTTAt                                                                                                                                                                                                                                  >  1:300606/1‑251 (MQ=255)   gTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTAtt                                                                                                                                                                                                                                 <  1:191160/251‑1 (MQ=255)             aaaCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGc                                                                                                                                                                                                                       >  1:409481/1‑251 (MQ=255)              aaCATGGGGGCACTGGGCGATCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATAAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATTTTTTAATTTAGCTGATACAACATACAGACGAATGAATAGCTACAACATTGAAACACAGGACGGGTTTATTAAGGTTAGGcc                                                                                                                                                                                                                      >  2:421693/1‑251 (MQ=255)                catggaagCACTGGGCGCTCTGTCTGTCGATCTCACGCAGGCTAATTCCACACTTCCCGCTGACAGTCAGCATGACGGACAATCGGTGCGTTTACTCTATAACAAATCGCTCACTGAAACAGGCACGAATATTAAGTTAGTGGGTTACCGTTATTCGACAAGCGGATATTTTAATTTAGCTGATCCAACATACAGTCGCATGAAAGGCTACAACATTTAAACACAGAACAGCATTATTCAGGCTAATCCg                                                                                                                                                                                                                     >  2:206303/8‑250 (MQ=255)                                         gtgGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCTGTGCGTTTTCTCTATAACAAATCGCCCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTGCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCg                                                                                                                                                                                            <  1:441882/250‑1 (MQ=255)                                               atGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                                                                                                                                                                                     <  1:368424/251‑1 (MQ=255)                                                             ttCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAt                                                                                                                                                                       >  1:418898/1‑251 (MQ=255)                                                              cccactctcctcgaagACAGTCAGCATGACGGACAATCGGTGTGTTTTTGCTATAACAAATCGCTCAATGAATCAGGAACGAATATTCAGTTAGTGGGTTACCGGTATTCGACCAGCTGATATTTTAAGTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCGTATACCAAACGCGGGAAAtt                                                                                                                                                                      <  2:460268/236‑1 (MQ=255)                                                                  cacTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGCATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGACATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACa                                                                                                                                                                   >  1:437393/1‑250 (MQ=255)                                                                   acTTCCCGATGACCGTCAGCATGACGGACAATCAGTGCGTTTTCTCCATAACAAATCGCTCAATGAATAAGGCACGAATATTCAGTTAGTGGGTTACCGATATTCGACCAGCGGATATTTTAATTTCGCTGATACAACCTACACTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGCAATTCACCGACTATTACAACCTCGCTTATAACAAAAGCGGTAAATTACaa                                                                                                                                                                  >  2:126662/1‑250 (MQ=255)                                                                                   cAGCATGACGGACAATCGGTGCGTTGTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGc                                                                                                                                                  <  1:172493/250‑1 (MQ=255)                                                                                       aTGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGATATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAAAACTCACCGTTACTCCACAACTc                                                                                                                                             >  2:65569/1‑251 (MQ=255)                                                                                           cGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGTAATTACAACTCACCGTTACTCAGCAACTCGGgc                                                                                                                                         >  1:174674/1‑251 (MQ=255)                                                                                                                           cGCTCAAGGAATCAGGCACTAATATTTAGTTAGTGGGTTACCGTTATTCGACCAGCGGATAATATAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGGTATTCAGGTTAAGCTGAAATTCACCGACTATTACACCCTCGCTTATAACAAACGCGTGAAATTACAATTCTCCGTTACTCAGCAACTCGGGCTCACATCAACACTGTATTTGAGTGGTAGCCAt                                                                                                         <  2:453559/251‑1 (MQ=255)                                                                                                                             cTCAATGAATCAGGCACGAATATTCAGTTCGTGGGTTACCCTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTAGTCAGGTTAATCCGAAATTCACCGACTATTACAACCTAGCTTATAACAAACGCGGCAAATTACAACTCACAGTTACTCCGCAACTCGGACGCACATCAACACTGAATTTGAGTGGTAGCCATc                                                                                                        >  2:26164/1‑250 (MQ=255)                                                                                                                                      acAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                                                                              <  1:240077/250‑1 (MQ=255)                                                                                                                                       cAGGCACGAATATTCAGTTAGTGTGTTACCTTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATAAATGGCTACAACATTGAAACACATTACGGAGt                                                                                                                                                                                                                                     <  2:208300/115‑1 (MQ=255)                                                                                                                                       cAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGt                                                                                                                                                                                                                                     >  1:208300/1‑115 (MQ=255)                                                                                                                                         ggCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                                                                           <  1:122192/251‑1 (MQ=255)                                                                                                                                            aCGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCGCGCTTATAACAAACGCGGGAAATTACAACTCTCCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTTGTAGCCATCAAACTTATTGGGGAAc                                                                                        <  1:65569/251‑1 (MQ=255)                                                                                                                                                             gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                                                                       <  1:376073/251‑1 (MQ=255)                                                                                                                                                               gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCa                                                                      >  2:347762/1‑250 (MQ=255)                                                                                                                                                                gtACCGATATTCGACCAGCGGATAATTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACAGTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACTCGGGAAATTACAACTCACCGTTACTCAGCAACTCGTGCGCACATCAACACTGTATTTGAGTGGTAGACATCAAACTTATTGGGGAACGAGTAATGGCGTTGAGCAAt                                                                    <  2:170117/250‑1 (MQ=255)                                                                                                                                                                  aCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTTAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTCTTTGAGTGGTAGCCATCAAACTTATTAGGGAACGCGTAATGTCGATGAGCAATTc                                                                  >  2:133307/1‑251 (MQ=255)                                                                                                                                                                                 gCGGATATTTTAATTTCGCTGATACAACATACAGACGAATGAATGGCTATAACATTTAAACACAGGACGGAGTTATTCAGGTTAAACCGAAATTCACCGGCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCTTTACTCAGCAACTCCGGCGCAAATCAACACAATATTTGCGTAGTAGCCATCAAAATTATAGGGGAACGAGGAATGTCGATGAGCACTTACAGGATGGATTAAAt                                                   >  2:394875/1‑251 (MQ=255)                                                                                                                                                                                     atatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTg                                               <  1:347762/251‑1 (MQ=255)                                                                                                                                                                                        ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTTAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAACAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAACCATCAAAATTATTGGTGAACGAGTAATGTCTATGATCAATTCCAGACTGGTTTaagtaatgggt                                            >  2:257345/1‑242 (MQ=255)                                                                                                                                                                                           taatttatcTGATACAACATACAGTCGAATTAAAGGCTACAACATTGAAACACAGGACTGAGTTATTCATGATAATCCGAAAATCACCGACTATTACACCCTCGCTTATAACAATCGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACTCTGTATTGGAGTGGTAGCCATCAAACTTATTGGGGAACGAGAAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTc                                          <  2:411430/242‑1 (MQ=255)                                                                                                                                                                                                               aCAGTCGAATGAATGGCTAAAACATTTAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCCCCGTCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGTGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATGCTGCGTTCGAAGCTATCAACTGGACGCt                      <  2:382295/250‑1 (MQ=255)                                                                                                                                                                                                                aattCGAATGAAGGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCa                    <  2:418898/248‑1 (MQ=255)                                                                                                                                                                                                                  gTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGc                  <  1:46965/251‑1 (MQ=255)                                                                                                                                                                                                                    cGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCt                 >  1:252074/1‑250 (MQ=255)                                                                                                                                                                                                                             ggCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGAc       >  1:179283/1‑251 (MQ=255)                                                                                                                                                                                                                              gCTACAACATTGAAACACAGTACGGAGTTATGCAGGTTAAGCCGAAATCCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACg      <  2:3004/251‑1 (MQ=255)                                                                                                                                                                                                                                 acaacaTTGAAACACAGGACGGAGTTATTCAGGGTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaaa   <  2:323618/251‑1 (MQ=255)                                                                                                                                                                                                                                   aacaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTCCAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTCCTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaaaa  >  2:287464/1‑250 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAA  >  NC_000913/4546370‑4546845

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