Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,064	C→T	P113S (CCA→TCA)	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,064	0	C	T	100.0%	147.7 / NA	47	P113S (CCA→TCA)	fimI	fimbrial protein involved in type 1 pilus biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/23);  total (24/23)

GAAGATAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAA  >  NC_000913/4543932‑4544192                                                                                                                                     |                                                                                                                                 gAAGATAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGa                                                                                                                            <  2:535243/139‑1 (MQ=255)   aGATAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATa                                                                                                                          <  2:325065/139‑1 (MQ=255)      tAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGcc                                                                                                                       <  1:325403/139‑1 (MQ=255)         cgcACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTACTTTCCGTGGGAGAGGGGTCAGGGATAGccacc                                                                                                                    >  2:571112/1‑139 (MQ=255)           cACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCaa                                                                                                                  >  2:513538/1‑139 (MQ=255)                       ttGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTACTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGc                                                                                                      <  1:571112/139‑1 (MQ=255)                        tGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCg                                                                                                     <  1:914071/139‑1 (MQ=255)                          ttATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCgtt                                                                                                   >  1:668138/1‑139 (MQ=255)                                                 aCGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCg                                                                            <  2:974799/139‑1 (MQ=255)                                                     tggtgAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTAcc                                                                        >  2:132402/1‑139 (MQ=255)                                                     tggtgAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTAcc                                                                        >  1:156963/1‑139 (MQ=255)                                                     tggtgAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTAcc                                                                        >  1:779914/1‑139 (MQ=255)                                                      ggtgAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGCTGATGAAGGAAACCTCGTACCg                                                                       >  1:545673/1‑139 (MQ=255)                                                       gtgagtgaACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGa                                                                                                                            >  1:387601/1‑85 (MQ=255)                                                       gtgagtgaACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGa                                                                                                                            <  2:387601/85‑1 (MQ=255)                                                       gtgagtgaACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGACACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGa                                                                      >  1:670021/1‑139 (MQ=255)                                                       gtgagtgaACGTGTAGGAGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGa                                                                      >  1:660273/1‑139 (MQ=255)                                                             gaACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTATAAATCAGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGCCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATc                                                                >  2:963826/1‑139 (MQ=255)                                                              aaCGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCg                                                               >  2:506580/1‑139 (MQ=255)                                                              aaCGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCg                                                               <  1:820056/139‑1 (MQ=255)                                                               aCGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGccacc                                                                                                                    >  2:465161/1‑85 (MQ=255)                                                               aCGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGccacc                                                                                                                    <  1:465161/85‑1 (MQ=255)                                                               aCGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATACGGATGTGCTTGCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGt                                                              >  2:447014/1‑139 (MQ=255)                                                                         gtgGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGCATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCaa                                                    >  2:718677/1‑139 (MQ=255)                                                                          tgGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCaaa                                                   <  1:506580/139‑1 (MQ=255)                                                                           ggCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAAc                                                  >  1:137688/1‑139 (MQ=255)                                                                            gCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACt                                                 >  2:1038117/1‑139 (MQ=255)                                                                            gCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACt                                                 >  1:913436/1‑139 (MQ=255)                                                                                 tCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGaaa                                            >  1:1006858/1‑139 (MQ=255)                                                                                 tCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGaaa                                            <  2:668138/139‑1 (MQ=255)                                                                                  cACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTtgatgat                                                                                         >  1:338538/1‑93 (MQ=255)                                                                                  cACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTtgatgat                                                                                         <  2:338538/93‑1 (MQ=255)                                                                                     ggTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGc                                        <  2:267042/139‑1 (MQ=255)                                                                                          cgcgGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAata                                                                                                                >  2:384814/1‑62 (MQ=255)                                                                                          cgcgGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAata                                                                                                                <  1:384814/62‑1 (MQ=255)                                                                                            cgGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTc                                 <  1:513538/139‑1 (MQ=255)                                                                                                 ggTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGct                            >  1:697809/1‑139 (MQ=255)                                                                                                  gTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGTGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGctc                           <  2:545673/139‑1 (MQ=255)                                                                                                   tAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGctct                          >  1:224948/1‑139 (MQ=255)                                                                                                                   ggctCCGTGGGCGAGGGGTCAGGGATAGCCACCAATATTTGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCa          <  2:268571/136‑1 (MQ=255)                                                                                                                   cTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTcg          <  2:660273/139‑2 (MQ=255)                                                                                                                   cTTTCCGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCa          <  2:913436/139‑1 (MQ=255)                                                                                                                        cGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGcc     <  1:447014/139‑1 (MQ=255)                                                                                                                        cGTGGGAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGcc     <  2:429418/139‑1 (MQ=255)                                                                                                                           gggAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTa                         <  1:9303/116‑1 (MQ=255)                                                                                                                           gggAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTa                         >  2:9303/1‑116 (MQ=255)                                                                                                                           gggAGAGGGGTCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCaaa  <  2:156963/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                 GAAGATAGCGCACCGGTGCCTTTTGTTATTCATTTACGGGAATGTAGCACGGTGGTGAGTGAACGTGTAGGTGTGGCGTTTCACGGTGTCGCGGATGGTAAAAATCCGGATGTGCTTTCCGTGGGAGAGGGGCCAGGGATAGCCACCAATATTGGCGTAGCGTTGTTTGATGATGAAGGAAACCTCGTACCGATTAATCGTCCTCCAGCAAACTGGAAACGGCTTTATTCAGGCTCTACTTCGCTACATTTCATCGCCAAA  >  NC_000913/4543932‑4544192

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