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

ATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAA  >  NC_000913/4546354‑4546845                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     aTCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGaa                                                                                                                                                                                                                                                   <  2:162767/251‑1 (MQ=255)    gTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAaca                                                                                                                                                                                                                                                <  1:65486/251‑1 (MQ=255)          ttAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGa                                                                                                                                                                                                                                           <  1:56250/250‑1 (MQ=255)                 cGGTATCGGGAAAAACATGGGGGCACTGGGCGCGCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTGATTCGACTAGCGTATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTa                                                                                                                                                                                                                                   <  2:365187/251‑1 (MQ=255)                    tATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATTACGCAGGCTCATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCCCGACTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGCTACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTc                                                                                                                                                                                                                                >  2:167838/1‑251 (MQ=255)                          gAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGtt                                                                                                                                                                                                                           >  1:127797/1‑250 (MQ=255)                               aCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                                                                                                                                                                                                     >  1:221312/1‑251 (MQ=255)                                     gggCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                                                                                                                                                                               >  2:233278/1‑251 (MQ=255)                                           tGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                                                                                                                                                                                         >  1:404658/1‑251 (MQ=255)                                                gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                     >  1:163180/1‑250 (MQ=255)                                                              tatGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                                                                                                                                                                                         >  1:120504/1‑232 (MQ=255)                                                              tatGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                                                                                                                                                                                         <  2:120504/232‑1 (MQ=255)                                                                                       cccGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGTGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa                                                                                                                                                              <  1:171911/250‑1 (MQ=255)                                                                                       cccGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAc                                                                                                                                                             <  2:258538/251‑1 (MQ=255)                                                                                           aTGACAGTCAGAATTAAGGACAATTGGTGCGTTTTCACTAAAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTGAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACGATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAATCACCGtt                                                                                                                                                         <  2:238972/251‑1 (MQ=255)                                                                                                aGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTc                                                                                                                                                     >  1:237889/1‑250 (MQ=255)                                                                                                                       gCGTTTTCTCGATAACAAATCGCGCAATGAATCAGGCACGAATATTCAGTGAGTGGGTTACCGTTATTCGACCAGCGGATATCTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAGTCACCGACTAGTACAGCCTCTCTTATAACAACCGCGGGAAATTACAACTCACCGTTACGCAGCAGCTCGGTCTCACATCGaca                                                                                                                              <  2:379782/250‑1 (MQ=255)                                                                                                                          ttttCTCTATAACAAATCGCTCTATGAATCAGGCACGCATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTTCAACCCAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCCACACTGt                                                                                                                          >  1:101050/1‑251 (MQ=255)                                                                                                                                    aaCAAAGCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCTCCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTgg                                                                                                                <  2:397651/251‑1 (MQ=255)                                                                                                                                       aaaTCGCTCAATGAATCAGGCACGAATATTCAGCTAGTGGGCCACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATAAAACCTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACCGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACACCTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGt                                                                                                               >  1:421394/1‑249 (MQ=255)                                                                                                                                          tcgcttaattaatccgggtcGAATAATCAGTTAGAGGGTTACCGTGATTCGACCAGCGGATATTTAAATTTCGCTGATACAACATACAGTCGAATGAATGGATACAACATTGAAACACAGGACGGAGTGTTCCAGGTTACGCCGACAGTCAACTACTATTACAACCTCTCTTATAACCAGCGCGGGAAAGTACAACTATCCGTGAGTCACCAAATCGGGGGCACATCAACCCTGGATTTGAGTGGTAggc                                                                                                           <  2:258996/231‑3 (MQ=255)                                                                                                                                                                     cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTGGGGCGGACATCAACACTGTATTTGAGTGTTAGTCATCAAACTTATTGGGTATCGAGTAATg                                                                                <  2:221312/250‑1 (MQ=255)                                                                                                                                                                            tGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCa                                                                                                                                 >  1:428494/1‑194 (MQ=255)                                                                                                                                                                            tGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACGCACCGTTACTCAGCAACTCGGGCGCACATCa                                                                                                                                 <  2:428494/194‑1 (MQ=255)                                                                                                                                                                                   ccGTTATTCGACCAGCGGATATTTTAATTTCGCTTATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGGCGGAGTTAGTCAGGTTAAGCCGAACTTCACCGTCTATTACACCCTCGACTATAACAAACGCGGTACATTACACCTCACCGTTACGCAGCAACTCGGGCGCATATCAAACCCGTACTTGAGTGGTAGCCAGCAAGAGTATGGGGGAACGAGTATTGTCGATGTGCAAtgcc                                                                 <  2:345476/251‑4 (MQ=255)                                                                                                                                                                                   ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGCCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGCCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAATGGCGCACGTGTAATGTCGAGTAGCAATTcc                                                                 >  2:133750/1‑251 (MQ=255)                                                                                                                                                                                    cGTTATTCGACCAGCGGATATTTTAATTTCGCGGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGTACGGAGTTATTCAGGGTAAGCCGAAAGTCACCTACTATTACAACCTCGCTTAGAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCATCTCGGGCGCAGATCAACACTGTATTTTAGTGGTAGCCATCAAACTTAGTGGGGAACGAGTAATGCCGATGAGCAATTCCa                                                                <  2:428822/251‑1 (MQ=255)                                                                                                                                                                                     gTTATTCGACGCGCGGATAATTTAATTTCGCTGATACAACATAAAGTCGAATGAATGGATACAAAATTGAATCACAGGCAGGAGTTATACAGGATAAGCCGAAAATCACCGACTATTAAAACCCCGCTTATAACAAACGCGGGAAAATACAAAGCACAGCTACTCAGCAACACGGACGCACAACAAAAATGCATTTGAGTGGTAGCCATAAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTcc                                                                 >  2:256534/1‑249 (MQ=255)                                                                                                                                                                                                           tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGGTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAGTTGAGTGGTAGCCATCAAACTTATTGGGGTACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTc                                          <  1:39149/250‑1 (MQ=255)                                                                                                                                                                                                              tttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGa                                        >  1:55907/1‑249 (MQ=255)                                                                                                                                                                                                                 cGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGat                                    >  2:53675/1‑250 (MQ=255)                                                                                                                                                                                                                 cGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTACTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGata                                   <  2:250978/251‑1 (MQ=255)                                                                                                                                                                                                                           acaTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGAc                         <  1:310711/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGa        <  2:312616/250‑1 (MQ=255)                                                                                                                                                                                                                                             ggCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAGTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGAc       <  2:128492/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aacaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGAAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaaaa  >  2:268553/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     ATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAA  >  NC_000913/4546354‑4546845

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