Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,088	G→A	G94E (GGG→GAG)	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,088	0	G	A	100.0%	104.5 / NA	33	G94E (GGG→GAG)	fimH	minor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/15);  total (18/15)

GAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGT  >  NC_000913/4548954‑4549223                                                                                                                                       |                                                                                                                                        gAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGAcc                                                                                                                                     <  1:749235/139‑1 (MQ=255)  aaTGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCg                                                                                                                                    <  2:762838/139‑1 (MQ=255)        gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGGAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTACTTTTTCCGAGCCCGTAAAat                                                                                                                              >  2:548706/1‑139 (MQ=255)             aaaCCTGGTCGTGCATCTTTCGACGCAAATCTTTTGCCATAACGCTTGTCAGGAAACCATGACAAACTATGTCACACTGCAACGAGGCTCGGCTTATTGCGGCGCGTTAACTAATTTTTCCGAGACCGTAAAATATAGt                                                                                                                         <  2:261287/139‑1 (MQ=255)               aCCTGGTCGTGGATCTTTCGACGCAAATCTTTTCCCATAACGATTATCCGTAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTgg                                                                                                                       <  2:524080/139‑1 (MQ=255)                ccTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGc                                                                                                                      <  1:101496/139‑1 (MQ=255)                         ggATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTAt                                                                                                             <  1:15326/139‑1 (MQ=255)                         ggATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTAt                                                                                                             >  2:321273/1‑139 (MQ=255)                         ggATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTAGCAGTAGCTAt                                                                                                             >  1:501191/1‑139 (MQ=255)                            tCTTTCTAAGCAAATCTTATGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGGTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                                                          <  2:228474/139‑1 (MQ=255)                                                gCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAc                                                                                      <  2:461782/139‑1 (MQ=255)                                                       cGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgc                                                                               <  1:45963/139‑1 (MQ=255)                                                           tATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgttg                                                                           >  1:360695/1‑139 (MQ=255)                                                                  aaaCCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCCCGCGTTGTTTATaa                                                                    >  1:475923/1‑139 (MQ=255)                                                                        ttACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACAACCAGCGAAACGCCGCGCGTTGTTTATAATTCgag                                                              >  1:425499/1‑139 (MQ=255)                                                                          aCAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgttgtt                                                                         <  1:649883/126‑1 (MQ=255)                                                                          aCAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgttgtt                                                                         >  2:649883/1‑126 (MQ=255)                                                                          aCAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACAACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa                                                            <  2:425499/139‑1 (MQ=255)                                                                                  tGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg                                                    >  1:112287/1‑139 (MQ=255)                                                                                     cacaCTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCg                                                 >  1:158579/1‑139 (MQ=255)                                                                                        acTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTgg                                              >  2:739747/1‑139 (MQ=255)                                                                                          tGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTa                                                                                                                                  <  1:760316/53‑1 (MQ=255)                                                                                          tGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTa                                                                                                                                  >  2:760316/1‑53 (MQ=255)                                                                                          tGCAACGAGGCTCGGCTGATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc                                            >  1:483146/1‑139 (MQ=255)                                                                                               cGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTgg                                       >  2:299887/1‑139 (MQ=255)                                                                                                  ggCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg                                          <  2:174025/133‑1 (MQ=255)                                                                                                  ggCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg                                          >  1:174025/1‑133 (MQ=255)                                                                                                    cTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTa                                                                                                                                  <  1:704072/43‑1 (MQ=255)                                                                                                    cTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTa                                                                                                                                  >  2:704072/1‑43 (MQ=255)                                                                                                    cTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCtt                                  <  2:483146/139‑1 (MQ=255)                                                                                                    cTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCtt                                  >  1:402345/1‑139 (MQ=255)                                                                                                              ggcggcGTGTTATCTAATTTTTCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGc                        >  2:532064/1‑139 (MQ=255)                                                                                                                               tttttCCGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTgcgggc       >  1:606488/1‑139 (MQ=255)                                                                                                                                    ccGAGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGt  >  2:477114/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                        GAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGT  >  NC_000913/4548954‑4549223

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