Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,191	A→G	K155R (AAA→AGA)	fimI →	fimbrial protein involved in type 1 pilus biosynthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,191	0	A	G	100.0%	89.1 / NA	24	K155R (AAA→AGA)	fimI	fimbrial protein involved in type 1 pilus biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/15);  total (9/15)

AGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAAT  >  NC_000913/4544059‑4544324                                                                                                                                     |                                                                                                                                      agGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCg                                                                                                                                 >  1:127257/1‑139 (MQ=255) agGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCg                                                                                                                                 >  2:166066/1‑139 (MQ=255)   gggCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGATACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTg                                                                                                                               >  2:317775/1‑139 (MQ=255)           gATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGCCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCggg                                                                                                                       >  1:329457/1‑139 (MQ=255)                                    ttgttTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCa                                                                                              <  2:83977/139‑1 (MQ=255)                                         tgatgatgaAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGcc                                                                                         <  1:498982/139‑1 (MQ=255)                                          gatgatgaAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGccc                                                                                        <  1:231223/139‑1 (MQ=255)                                              atgaAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGc                                                                                    <  2:460984/139‑1 (MQ=255)                                              atgaAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGc                                                                                    <  1:317775/139‑1 (MQ=255)                                                  aGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTgg                                                                                <  2:92329/139‑1 (MQ=255)                                                   ggAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGt                                                                               <  1:166066/139‑1 (MQ=255)                                                            gTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTaa                                                                      <  1:361210/139‑1 (MQ=255)                                                             tACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAAc                                                                     <  1:372913/139‑1 (MQ=255)                                                             tACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAAc                                                                     <  2:127257/139‑1 (MQ=255)                                                                cGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTa                                                                  <  1:179949/139‑1 (MQ=255)                                                                            ctccAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGcc                                                                                   >  1:46276/1‑110 (MQ=255)                                                                            ctccAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGcc                                                                                   <  2:46276/110‑1 (MQ=255)                                                                               cAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTcag                                                   >  1:245832/1‑139 (MQ=255)                                                                                 gCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTcagca                                                 >  1:14165/1‑139 (MQ=255)                                                                                               gCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATaac                                   <  2:245832/139‑1 (MQ=255)                                                                                                 tttACTCAGGCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATaacag                                 <  2:329457/139‑1 (MQ=255)                                                                                                         ggCTCTACTTCGCTACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAacaggaacagga                         <  2:144310/139‑1 (MQ=255)                                                                                                                      tACATTTCATCGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGtaataa            >  2:280960/1‑139 (MQ=255)                                                                                                                                cGCCAGATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAAt  >  1:307466/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                      AGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAAT  >  NC_000913/4544059‑4544324

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: ‑