Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,188,957	(C)5→4	coding (1954/2481 nt)	yehB ←	putative fimbrial usher protein YehB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,188,953	0	C	.	100.0%	214.5 / NA	51	coding (1958/2481 nt)	yehB	putative fimbrial usher protein YehB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/30);  total (21/30)

TATCACCACTCCATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACT  >  NC_000913/2188814‑2189083                                                                                                                                            |                                                                                                                                   tatCACCACTCCATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGc                                                                                                                                    <  2:578318/140‑1 (MQ=255)       cacTCCATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGcccctga                                                                                                                              <  1:819938/140‑4 (MQ=255)         cTCCATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCCCCCAACCAGGCACCAAT‑CCCCCTaaaa                                                                                                                           >  1:507226/1‑140 (MQ=255)           ccATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAaca                                                                                                                         >  1:196670/1‑140 (MQ=255)            cATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGCTTATCGCTACCCGACCTGGCGGCAATG‑CCCCTGAAacac                                                                                                                        >  2:323584/1‑140 (MQ=255)            cATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAacac                                                                                                                        <  2:575519/140‑1 (MQ=255)            cATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGACAGATTAACGCAACCCGACCGGGCGTCAATG‑CCCATGAAacac                                                                                                                        >  2:382109/1‑140 (MQ=255)                  ggTTTGTTGTGCGATCTTTTTGGCCGTTGGCATAAGCTTCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTc                                                                                                                  <  1:347611/140‑1 (MQ=255)                    tttGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCa                                                                                                                >  1:85343/1‑140 (MQ=255)                      tgttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGc                                                                                                              <  2:184332/140‑1 (MQ=255)                      tgttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGc                                                                                                              <  1:639311/140‑1 (MQ=255)                      tgttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGc                                                                                                              >  2:203754/1‑140 (MQ=255)                      tgttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCAGACCAGGGGACAATG‑CCCCTGAAACACTGGCTCCAGc                                                                                                              >  1:397623/1‑140 (MQ=255)                      tgttgtGCGATATTTATGCCCATTGTCATAAGCATCTCTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGc                                                                                                              <  1:258482/140‑1 (MQ=255)                       gttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGcc                                                                                                             <  1:419301/140‑1 (MQ=255)                         tgttccataaTTTTGCCCCTTGGCATAAGCAACTTTAATTCCTTGCGCCTCCCACACCGCAAACGGTTCCGAAAGGCGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTg                                                                                                           <  2:507226/134‑1 (MQ=255)                          gtgCGATATTTTTGCCCATTGATATAAGCATCTTTAACTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGt                                                                                                          <  1:323584/140‑1 (MQ=255)                            gCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGc                                                                                                              <  2:806833/134‑1 (MQ=255)                            gCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGc                                                                                                              >  1:806833/1‑134 (MQ=255)                            gCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCg                                                                                                        >  2:824345/1‑140 (MQ=255)                                  tttttGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGt                                                                                                  >  1:626477/1‑140 (MQ=255)                                    tttGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTAc                                                                                                <  2:85343/140‑1 (MQ=255)                                     ttGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACt                                                                                               <  1:688061/140‑1 (MQ=255)                                     ttGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACt                                                                                               <  1:824345/140‑1 (MQ=255)                                        cccATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCga                                                                                            <  1:660126/140‑1 (MQ=255)                                        cccATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCga                                                                                            <  2:184956/140‑1 (MQ=255)                                         ccATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCgac                                                                                           <  1:637345/140‑1 (MQ=255)                                          cATTGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCgact                                                                                          <  2:650541/140‑1 (MQ=255)                                            ttGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCgactga                                                                                        >  2:548742/1‑140 (MQ=255)                                            ttGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCga                                                                                            >  1:665765/1‑136 (MQ=255)                                            ttGACATAAGCATCTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCga                                                                                            <  2:665765/136‑1 (MQ=255)                                                   aaGCATCTGTAATTCCGGGCGCGTCCATCACCGCAAACGATTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAAc                                                                                 <  2:397623/140‑1 (MQ=255)                                                         cTTTAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCAt                                                                           >  2:743826/1‑140 (MQ=255)                                                           ttAATTCCTGGCGCATCCATCACAGTAAAGGTTTCGGCAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTc                                                                         <  1:320588/140‑1 (MQ=255)                                                            tAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCa                                                                        <  2:775470/140‑1 (MQ=255)                                                            tAATTCCTGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCa                                                                        <  2:666838/140‑1 (MQ=255)                                                                   tGGCGCATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGATTGACTATAACTGCCATTCACTGTCGc                                                                 >  2:504402/1‑140 (MQ=255)                                                                      cgcATCCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAAc                                                              <  1:743826/140‑1 (MQ=255)                                                                                acaGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTc                                                    <  2:626477/140‑1 (MQ=255)                                                                                     aaaCGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGt                                               <  2:17604/140‑1 (MQ=255)                                                                                       aCGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCa                                             >  2:41861/1‑140 (MQ=255)                                                                                         gTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCaaa                                           <  1:41861/140‑1 (MQ=255)                                                                                                   gACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAg                                 >  2:477108/1‑140 (MQ=255)                                                                                                     cGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCGGCt                               >  1:675151/1‑140 (MQ=255)                                                                                                           gCCAGATTAACGCCAACCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGtt                         >  1:491607/1‑140 (MQ=255)                                                                                                             cAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTTCAGGTCAAATTCGCCCCAGCTGTCGTTTc                       >  1:3610/1‑140 (MQ=255)                                                                                                             cAGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTAACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTc                       >  1:92886/1‑140 (MQ=255)                                                                                                              aGATTAACGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCa                      <  1:422896/140‑1 (MQ=255)                                                                                                                   aaCGCCACCCGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGATTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTcc                 <  1:504402/140‑1 (MQ=255)                                                                                                                      gCCACCCGTCCAGGCGATAATG‑ACCCTGAAACACTTGCTCCAGCCTGTAGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTg              <  2:393033/140‑1 (MQ=255)                                                                                                                          cccGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGtt                         >  1:46067/1‑125 (MQ=255)                                                                                                                          cccGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGtt                         <  2:46067/125‑1 (MQ=255)                                                                                                                          cccGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGGCAAATTCGCCCTAGCTGTCGTTTCATTTCCCTGATGt          >  1:352843/1‑140 (MQ=255)                                                                                                                            cGACCAGGCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTtg        <  2:196670/140‑1 (MQ=255)                                                                                                                                  ggCGACAATG‑CCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTTCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACt  <  2:3610/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                   TATCACCACTCCATTACGGTTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACT  >  NC_000913/2188814‑2189083

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