Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,362,459	C→T	S116S (AGC→AGT)	yhcA →	putative fimbrial chaperone YhcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,362,459	0	C	T	100.0%	110.7 / NA	33	S116S (AGC→AGT)	yhcA	putative fimbrial chaperone YhcA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/13);  total (20/13)

GTGGTGAAGCTAGTATTAATATTAAGAATACAGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGCAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAAAACGGGGGCTATAAATAAGACCTATTTTCTACCACACACCAGTACAACTACTGCATTAACAAA  >  NC_000913/3362215‑3362704                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      gtggtgAAGCTAGTATTAATATTAAGAATACAGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGtt                                                                                                                                                                                                                                                 >  1:383132/1‑251 (MQ=255)       aaGCTAGTATTAATATTAAGAATACAGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGATCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACACGAGCAGTAGAGTTACTGtt                                                                                                                                                                                                                                           >  1:278633/1‑251 (MQ=255)        aGCTAGTATTAATATTAAGAATACAGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGttt                                                                                                                                                                                                                                          <  2:85466/251‑1 (MQ=255)        aGCTAGTATTAATATTAAGAATACAGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGttt                                                                                                                                                                                                                                          <  1:327859/251‑1 (MQ=255)             gTATTAATATTAAGAATACAGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAAtt                                                                                                                                                                                                                                     <  1:285901/251‑1 (MQ=255)                       taggaagttagatgaacacacTTCATTGCTTTACACCACCATAGTTGATTTACCTGTAAGCAATAAATAATTTCGATTAATCCCGACTCAGCCCGTTACTCCTGTAGAACAAGGTCAAGCACAACAAGTTCGATATTTACTTCCGGCCACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGCTGACAGGAGCAGTAGAGTTACTGTTTAAATTAGTCAGGATc                                                                                                                                                                                                                           >  2:255517/21‑251 (MQ=255)                                aGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTAc                                                                                                                                                                                                                  >  1:221076/1‑251 (MQ=255)                                 gatgatCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACt                                                                                                                                                                                                                 >  1:184817/1‑251 (MQ=255)                                    gatCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGCGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGAt                                                                                                                                                                                                              >  2:213629/1‑251 (MQ=255)                                    gatCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATATACCTGTACTGAt                                                                                                                                                                                                              >  2:319157/1‑251 (MQ=255)                                                     ttACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGc                                                                                                                                                                                             >  1:210405/1‑251 (MQ=255)                                                            accaTTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGa                                                                                                                                                                                         >  1:135558/1‑248 (MQ=255)                                                                 tgttgATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCg                                                                                                                                                                                   >  1:348521/1‑249 (MQ=255)                                                                  gttgACTTGCATGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGCTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAAAGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTAACTTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCg                                                                                                                                                                                   >  2:255300/1‑248 (MQ=255)                                                                           ccTGAAAGCAATAAATCAATTCGTTTAATCCCGGCTCAGCCCGTTAGTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCTGAAGAGCGAGAGACCTGGaa                                                                                                                                                                       <  2:2385/251‑1 (MQ=255)                                                                                           cAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCACGTACAACAGGTTCGATTTTTACTTCAGGCAACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACATGTACTGATCCACCCGGCTTCCCTGCCGTCAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGt                                                                                                                                                       >  2:338741/1‑251 (MQ=255)                                                                                             aTTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTaa                                                                                                                                                     <  2:221076/251‑1 (MQ=255)                                                                                                         ccGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGAcc                                                                                                                                          >  1:390527/1‑250 (MQ=255)                                                                                                                                  aaGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGAGACCAAATTGAGATTAGCAAACCCAGTAAt                                                                                                                >  1:261703/1‑251 (MQ=255)                                                                                                                                     ccGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTAt                                                                                                             >  1:328470/1‑251 (MQ=255)                                                                                                                                       ggTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGt                                                                                                           <  1:36871/251‑1 (MQ=255)                                                                                                                                            aGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGtcgttc                                                                                                      <  2:348521/251‑1 (MQ=255)                                                                                                                                             gTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGGCCAAATTGAGATTAGCAACCCCAGTAATTATGtcgttcg                                                                                                     <  1:23823/251‑1 (MQ=255)                                                                                                                                                                         cGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCt                                                                         <  2:392835/251‑1 (MQ=255)                                                                                                                                                                          gACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCtt                                                                        >  1:284441/1‑251 (MQ=255)                                                                                                                                                                                   tctGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAAtt                                                                                                               <  2:314759/203‑1 (MQ=255)                                                                                                                                                                                   tctGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAAtt                                                                                                               >  1:314759/1‑203 (MQ=255)                                                                                                                                                                                      gCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCgaga                                                                                                                                                                                >  1:370640/1‑135 (MQ=255)                                                                                                                                                                                      gCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCgaga                                                                                                                                                                                <  2:370640/135‑1 (MQ=255)                                                                                                                                                                                                                             tCCACCAAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAAAACGGGGGCTATAAATAAGACCTATTTTCTACCACACACCAGt                      >  1:192120/1‑250 (MQ=255)                                                                                                                                                                                                                              ccaccaAAAGATGACAAGAGCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAAAACGGGGGCTATAAATAAGACCTATTTTCTACCACACACCAGTa                     <  1:83298/250‑1 (MQ=255)                                                                                                                                                                                                                                             aaGNGCAGTAGAGTTACTGTTTCAAATCGTCAGGATCTTCCTGTTCTGAACCCCCCCGCTTTCCTGGCGGAAGAGCGAGAGACCTTGAAAATCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAAAACGGGGGCTATAAATAAGACCTATTTTCTACCACACACCAGTACAACTACTGCATTAAc     <  2:383132/251‑1 (MQ=255)                                                                                                                                                                                                                                              agagCAGTAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAAAACGGGGGCTATAAATAAGACCTATTTTCTACCACACACCAGTACAACTACTGCATTAAc     >  1:10406/1‑250 (MQ=255)                                                                                                                                                                                                                                                agcagtAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAAGTTAAAACGCTCCCTTCAGGAAAAACGGGGGCTATAAATAAGACCTATTTTCTACCACACACCAGTACAACTACTGCATTAACaaa  <  1:255300/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      GTGGTGAAGCTAGTATTAATATTAAGAATACAGATGATCATCCTTCATTGCTTTACACCACCATTGTTGATTTGCCTGAAAGCAATAAATCAATTCGTTTAATCCCGACTCAGCCCGTTATTCGTGTTGAAGCCGGTCAAGTACAACAGGTTCGATTTTTACTTCAGGCGACTGTTCCTCTGCAATCAGAAGAGCTGAAACGCGTAACGTTCGAAGGTATTCCACCAAAAGATGACAAGAGCAGCAGAGTTACTGTTTCAATTCGTCAGGATCTACCTGTACTGATCCACCCGGCTTCCCTGCCGGAAGAGCGAGAGACCTGGAAATTCCTGGAATGGCGTAAAAATGGCGACCAAATTGAGATTAGCAACCCCAGTAATTATGTCGTTCGTATGACCTTACAATTTAAAACGCTCCCTTCAGGAAAAACGGGGGCTATAAATAAGACCTATTTTCTACCACACACCAGTACAACTACTGCATTAACAAA  >  NC_000913/3362215‑3362704

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