Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,061	A→G	Y85C (TAT→TGT)	fimH →	minor component of type 1 fimbriae

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,061	0	A	G	100.0%	177.9 / NA	51	Y85C (TAT→TGT)	fimH	minor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/27);  total (24/27)

GTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTT  >  NC_000913/4548928‑4549191                                                                                                                                      |                                                                                                                                   gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                               <  1:860586/139‑1 (MQ=255)      tGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgtg                                                                                                                          <  1:910056/139‑1 (MQ=255)          aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTAt                                                                                                                      >  1:128362/1‑139 (MQ=255)          aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTAt                                                                                                                      >  2:501158/1‑139 (MQ=255)           aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATc                                                                                                                     >  2:400736/1‑139 (MQ=255)               tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAAt                                                                                                                 >  1:872183/1‑139 (MQ=255)                       tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                         >  2:458163/1‑139 (MQ=255)                         gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  1:1068476/139‑1 (MQ=255)                         gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  1:458163/139‑1 (MQ=255)                         gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  2:128362/139‑1 (MQ=255)                         gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  2:872183/139‑1 (MQ=255)                           gAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                                                                                     <  2:927915/139‑1 (MQ=255)                                  gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                               >  1:255560/1‑106 (MQ=255)                                  gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggcg                                                                                                                               <  2:255560/106‑1 (MQ=255)                                         aCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTgg                                                                                       <  2:995492/139‑1 (MQ=255)                                         aCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTgg                                                                                       <  1:890336/139‑1 (MQ=255)                                         aCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTgg                                                                                       >  2:890336/1‑139 (MQ=255)                                             ggTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       >  2:752621/1‑119 (MQ=255)                                             ggTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  1:752621/119‑1 (MQ=255)                                                      tCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                          <  1:560708/139‑1 (MQ=255)                                                       cTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                         <  1:486915/139‑1 (MQ=255)                                                       cTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAt                                                                         <  1:868341/139‑1 (MQ=255)                                                                        ttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGaa                                                        <  1:501158/139‑1 (MQ=255)                                                                                        ccGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgtg                                                                                                                          <  2:99196/57‑1 (MQ=255)                                                                                        ccGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCgtg                                                                                                                          >  1:99196/1‑57 (MQ=255)                                                                                                 aTTACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       >  2:1017738/1‑67 (MQ=255)                                                                                                 aTTACAGACAATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  1:1017738/67‑1 (MQ=255)                                                                                                  ttACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgag                              >  2:514655/1‑139 (MQ=255)                                                                                                   tACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggc                                                                                                                                >  2:611932/1‑40 (MQ=255)                                                                                                   tACAGACTATGTCACACTGCAACGAGGCTCGGCTTGTggc                                                                                                                                <  1:611932/40‑1 (MQ=255)                                                                                                            tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       >  2:457478/1‑56 (MQ=255)                                                                                                            tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  1:457478/56‑1 (MQ=255)                                                                                                            tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGAt                       >  2:726826/1‑136 (MQ=255)                                                                                                            tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGAt                       <  1:726826/136‑1 (MQ=255)                                                                                                            tGTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg                    >  1:845480/1‑139 (MQ=255)                                                                                                             gTCACACTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGc                   >  2:1017108/1‑139 (MQ=255)                                                                                                                   cTGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGc             <  1:514655/139‑1 (MQ=255)                                                                                                                    tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACTACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc            >  2:749046/1‑139 (MQ=255)                                                                                                                    tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc            >  2:381917/1‑139 (MQ=255)                                                                                                                    tGCAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc            >  2:730005/1‑139 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       >  1:66916/1‑46 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                       <  2:66916/46‑1 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAatat                                                                                            >  2:4269/1‑57 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAatat                                                                                            <  1:4269/57‑1 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                                 >  2:1008482/1‑68 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                                 <  1:1008482/68‑1 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                              <  1:594241/103‑1 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                              >  2:594241/1‑103 (MQ=255)                                                                                                                      cAACGAGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg          >  2:350974/1‑139 (MQ=255)                                                                                                                           aGGCTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGgcg     <  1:749046/139‑1 (MQ=255)                                                                                                                              cTCGGCTTGTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCtt  <  1:350974/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   GTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTT  >  NC_000913/4548928‑4549191

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