Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,363,653	C→T	L283L (CTG→TTG)	yhcD →	putative fimbrial usher protein YhcD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,363,653	0	C	T	100.0%	115.0 / NA	32	L283L (CTG→TTG)	yhcD	putative fimbrial usher protein YhcD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/14);  total (18/14)

ACATTCAGTAACCAGAACTCGTCAACCTACCTTCAACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTCTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACAC  >  NC_000913/3363407‑3363890                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              acaTTCAGTAACCAGAACTCGTCAACCTACCTTCAACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGt                                                                                                                                                                                                                                           >  2:204114/1‑251 (MQ=255)     ccAGTAACCAGAACTCGTCAACCTACCTTCAACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATGCCATTTTGGTTCCt                                                                                                                                                                                                                                       <  1:315167/250‑1 (MQ=255)      cAGTAACCAGAACTCGTCAACCTACCTTCAACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGGCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTg                                                                                                                                                                                                                                      <  1:254258/251‑1 (MQ=255)          aaCCAGAACACGTCAACCTTCCTTCAACGTACCTTTACAGATCTTAAAACACTCATGCGTGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCCGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGCCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCgg                                                                                                                                                                                                                                   <  2:144539/250‑1 (MQ=255)                  cTCGTCAACCTACATTCAACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGCGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTc                                                                                                                                                                                                                           >  2:226320/1‑250 (MQ=255)                                     cGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGt                                                                                                                                                                                                       >  2:101148/1‑251 (MQ=255)                                      gTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAATGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGtt                                                                                                                                                                                                      >  1:212034/1‑251 (MQ=255)                                         ccTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGAATTGCAAACCAGCGGCAGTGGTGTGCACGTTTCTGGTATAGACAAAACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCATGATGTCCCTGTTCg                                                                                                                                                                                                    >  1:2533/1‑250 (MQ=255)                                            ttACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAAt                                                                                                                                                                                                >  2:190870/1‑251 (MQ=255)                                                     ttAAAACACTCATGCGAGCAGGTGAAGTTACCCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAANAGt                                                                                                                                                                                       >  2:380664/1‑251 (MQ=255)                                                                     aGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCg                                                                                                                                                                       <  1:204114/251‑1 (MQ=255)                                                                      gCAGGTGAAGTTAACCTCACTAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGCAATCGAACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGCGTCTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCACTGGTACTAGTGATCTTAATGTCACCGt                                                                                                                                                                      >  1:374156/1‑251 (MQ=255)                                                                         ggTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGt                                                                                                                                                                   <  1:1469/251‑1 (MQ=255)                                                                           tGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCg                                                                                                                                                                 >  1:358051/1‑251 (MQ=255)                                                                              aGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGaa                                                                                                                                                               >  1:290307/1‑250 (MQ=255)                                                                                   aCCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGAc                                                                                                                                                         >  2:210712/1‑251 (MQ=255)                                                                                            ataGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCg                                                                                                                                                >  1:317972/1‑251 (MQ=255)                                                                                                             gCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAATGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTcc                                                                                                                               <  2:212034/251‑1 (MQ=255)                                                                                                                 gTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTcc                                                                                                                           <  2:317972/251‑1 (MQ=255)                                                                                                                          gTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTNttt                                                                                                                  >  2:379407/1‑251 (MQ=255)                                                                                                                            aTCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTGCCCTCCACCCTGTTTa                                                                                                                 <  1:152438/250‑1 (MQ=255)                                                                                                                                                         aaCCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATAtt                                                                                                                                    >  1:317678/1‑202 (MQ=255)                                                                                                                                                         aaCCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATAtt                                                                                                                                    <  2:317678/202‑1 (MQ=255)                                                                                                                                                             aGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCtt                                                                               <  1:46191/251‑1 (MQ=255)                                                                                                                                                              gCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTc                                                                              >  2:197239/1‑251 (MQ=255)                                                                                                                                                                    atggtgTGCAAGGTACTGGTATAGCCAACACCTCTCAGGCGCTTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATACACTATCCCTGATGTACCTGTGCGCAAGGGTAATAGTTATCTTAATGTCACCGATGTCGAAACAGACGGTAGTGCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTATAAAGCTTCCAGGTTTATCGCTTCGCGATa                                                                        <  1:249981/250‑1 (MQ=255)                                                                                                                                                                                                            tCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAg                                                                                                            <  2:93058/175‑1 (MQ=255)                                                                                                                                                                                                            tCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAg                                                                                                            >  1:93058/1‑175 (MQ=255)                                                                                                                                                                                                              gtgtCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAAtt                              <  1:190870/251‑1 (MQ=255)                                                                                                                                                                                                                  cGAGATTCGTCAACAACGAGTTTTAATTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCTCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTg                          >  2:273304/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaTTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGacac  >  2:94181/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaTTCATTCCATTTTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGaca   <  1:98436/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              ACATTCAGTAACCAGAACTCGTCAACCTACCTTCAACGTACCTTTACAGATCTTAAAACACTCATGCGAGCAGGTGAAGTTAACCTCAATAATAGCGTGTTGGAAGGAGCCAGTATTTACGGTATCGAAATCGCACCGGACAACGCATTGCAAACCAGCGGCAGTGGTGTGCAAGTTACTGGTATAGCCAACACCTCTCAGGCTCGTGTCGAGATTCGTCAACAAGGAGTTTTAATTCATTCCATTCTGGTTCCTGCGGGCGCATTCACTATCCCTGATGTACCTGTTCGCAATGGTAATAGTGATCTTAATGTCACCGTTGTCGAAACAGACGGTAGTTCGCACAACTATATTGTTCCCTCCACCCTGTTTAATCAGCATGTAGAAAGCTTCCAGGGTTATCGCTTCGCGATAGGGCGGGTAGACGATGACTATGACGAATCACCTTGGGTAATTAGTGCATCGAGCGGATGGAATCTGACAC  >  NC_000913/3363407‑3363890

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