Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,187,468	A→C	I316R (ATA→AGA)	yehA ←	putative fimbrial‑like adhesin protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,187,468	0	A	C	100.0%	196.4 / NA	53	I316R (ATA→AGA)	yehA	putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/27);  total (26/27)

GGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCC  >  NC_000913/2187321‑2187614                                                                                                                                                    |                                                                                                                                                   gggAAATGCAGACGCCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACTTTGATTAATGCGGTTGCCTCCAGTGCGACAGATTTAACTTTGTTTGTATCTTAGACGTAGTAACTGGCTGTTCTCg                                                                                                                                                 <  1:32742/151‑1 (MQ=255)  ggAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCgg                                                                                                                                                >  1:875599/1‑151 (MQ=255)   gAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGa                                                                                                                                               <  2:20465/151‑1 (MQ=255)   gAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGa                                                                                                                                               <  2:452812/151‑1 (MQ=255)      aTGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAAtt                                                                                                                                            <  2:430665/151‑1 (MQ=255)        gCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTg                                                                                                                                           <  1:65632/150‑1 (MQ=255)               ccTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGtt                                                                                                                                   >  1:720970/1‑151 (MQ=255)                 ttGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGtt                                                                                                                                   >  2:67210/1‑149 (MQ=255)                   gTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGCTTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTAt                                                                                                                               >  1:463864/1‑151 (MQ=255)                      aTTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTg                                                                                                                            >  1:884698/1‑151 (MQ=255)                       ttttttATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGt                                                                                                                           >  1:689958/1‑151 (MQ=255)                          tttATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTa                                                                                                                          >  1:33828/1‑149 (MQ=255)                           ttattattaGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTaa                                                                                                                         <  1:586089/149‑1 (MQ=255)                             attattaGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGtt                                                                                                                      <  2:120907/150‑1 (MQ=255)                             attattaGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGtt                                                                                                                      <  2:200760/150‑1 (MQ=255)                               tattaGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTg                                                                                                                   >  2:331186/1‑151 (MQ=255)                                  taGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGAc                                                                                                                 <  2:513461/150‑1 (MQ=255)                                          tGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTcgccgc                                                                                                         >  1:347602/1‑150 (MQ=255)                                              cGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACt                                                                                                     <  1:122612/150‑1 (MQ=255)                                               gTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTc                                                                                                   >  2:697296/1‑151 (MQ=255)                                                  ttgcttgTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCtt                                                                                                >  2:195585/1‑151 (MQ=255)                                                  ttgcttgTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCtt                                                                                                >  1:799420/1‑151 (MQ=255)                                                    gcttgTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCttt                                                                                               <  1:379464/150‑1 (MQ=255)                                                      ttgTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTaa                                                                                             >  2:297992/1‑150 (MQ=255)                                                      ttgTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTaa                                                                                             >  2:320360/1‑150 (MQ=255)                                                              aGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGtat                                                                                     <  1:214196/150‑1 (MQ=255)                                                                  gTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCt                                                                                >  1:604572/1‑151 (MQ=255)                                                                   tATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCtt                                                                               >  1:369581/1‑151 (MQ=255)                                                                           gTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTActgctg                                                                       <  2:812233/151‑1 (MQ=255)                                                                              gATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGtt                                                                     >  1:93503/1‑150 (MQ=255)                                                                                   aTGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGat                                                               <  2:560662/151‑1 (MQ=255)                                                                                       ggTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATAttt                                                           <  1:733159/151‑1 (MQ=255)                                                                                           gCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTaaaa                                                       <  2:575146/151‑1 (MQ=255)                                                                                               ccAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAgtgt                                                   <  2:640785/151‑1 (MQ=255)                                                                                               ccAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAgtg                                                    <  1:153272/150‑1 (MQ=255)                                                                                                    gcgcCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTg                                              >  1:434284/1‑151 (MQ=255)                                                                                                     cgcCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTgg                                             <  1:378256/151‑1 (MQ=255)                                                                                                       ccAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCg                                           >  1:128395/1‑151 (MQ=255)                                                                                                        cAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTTCTGTTAGTGATATTTAAAAGTGTGACTGGCg                                           >  2:331450/1‑150 (MQ=255)                                                                                                         aGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGaa                                         <  1:318695/151‑1 (MQ=255)                                                                                                          gATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAAc                                        >  2:65294/1‑151 (MQ=255)                                                                                                               aaCTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCAtt                                    <  1:139244/150‑1 (MQ=255)                                                                                                                 cTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACt                                 >  1:501629/1‑151 (MQ=255)                                                                                                                     gtttgtATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTc                              <  1:820663/150‑1 (MQ=255)                                                                                                                        tgtATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCa                             >  2:443063/1‑148 (MQ=255)                                                                                                                        tgtATCGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAg                           <  1:120604/150‑1 (MQ=255)                                                                                                                             cGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATaa                      <  2:303361/150‑1 (MQ=255)                                                                                                                             cGTAGACGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATaa                      <  1:464356/150‑1 (MQ=255)                                                                                                                                  aCGTAGTAACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCt                 >  1:550075/1‑150 (MQ=255)                                                                                                                                    gtagtaACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACg              <  1:615558/151‑1 (MQ=255)                                                                                                                                       gtaACTGGCTGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGcc            >  1:678666/1‑150 (MQ=255)                                                                                                                                                tGTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTcc  >  1:125591/1‑151 (MQ=255)                                                                                                                                                 gTTCTCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTc   <  2:593871/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCC  >  NC_000913/2187321‑2187614

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