Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	152,036	C→T	V66I (GTA→ATA)	yadL ←	putative fimbrial‑like adhesin protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	152,036	0	C	T	100.0%	210.2 / NA	65	V66I (GTA→ATA)	yadL	putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/36);  total (29/36)

AGCAGTTCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATACATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTGCCGCATAAGAGAGTACG  >  NC_000913/151905‑152163                                                                                                                                    |                                                                                                                                aGCAGTTCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCa                                                                                                                          <  1:979973/139‑1 (MQ=255)                tgtgGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGaa                                                                                                           <  2:1092173/139‑1 (MQ=255)                     cTCGCTTGGTTAATTATATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATcaacaa                                                                                                                <  2:248882/129‑1 (MQ=255)                     cTCGCTTGGTTAATTATATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATcaacaa                                                                                                                >  1:248882/1‑129 (MQ=255)                                    ttATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTc                                                                                                      >  1:1032363/1‑124 (MQ=255)                                    ttATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTc                                                                                                      <  2:1032363/124‑1 (MQ=255)                                      aTTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTa                                                                                     >  1:1095333/1‑139 (MQ=255)                                          gcgcTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTaa                                                                                 <  2:891180/139‑1 (MQ=255)                                            gcTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATg                                                                               >  1:92567/1‑139 (MQ=255)                                                    aTTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTa                                                                       <  1:556757/139‑1 (MQ=255)                                                          gATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTa                                                                 <  1:704173/139‑1 (MQ=255)                                                                  cacaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATcaacaa                                                                                                                >  2:991988/1‑84 (MQ=255)                                                                  cacaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATcaacaa                                                                                                                <  1:991988/84‑1 (MQ=255)                                                                  cacaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTc                                                         >  1:415976/1‑139 (MQ=255)                                                                   acaGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCt                                                        <  2:92567/139‑1 (MQ=255)                                                                     aGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                      <  2:535129/139‑1 (MQ=255)                                                                     aGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                      <  1:92369/139‑1 (MQ=255)                                                                      gTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCAtt                                                                                       >  2:899400/1‑105 (MQ=255)                                                                      gTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCAtt                                                                                       <  1:899400/105‑1 (MQ=255)                                                                      gTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTa                                                       <  1:132981/137‑1 (MQ=255)                                                                      gTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTa                                                       >  2:132981/1‑137 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGt                                                                                                       >  1:63788/1‑88 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGt                                                                                                       <  2:63788/88‑1 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGa                                                                                                  >  2:790518/1‑93 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGa                                                                                                  <  1:790518/93‑1 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAcc                                                               >  2:98156/1‑128 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAcc                                                               <  1:98156/128‑1 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACag                                                    >  1:1201202/1‑139 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACag                                                    >  1:121505/1‑139 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACag                                                    >  1:286176/1‑139 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACag                                                    >  1:955591/1‑139 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACag                                                    >  1:173548/1‑139 (MQ=255)                                                                       tCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACag                                                    >  2:642943/1‑139 (MQ=255)                                                                        cTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACagt                                                   >  1:982788/1‑139 (MQ=255)                                                                         tttGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGc                                                                           <  1:68852/114‑1 (MQ=255)                                                                         tttGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGc                                                                           >  2:68852/1‑114 (MQ=255)                                                                           tGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAcc                                                               >  1:1017045/1‑124 (MQ=255)                                                                           tGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAcc                                                               <  2:1017045/124‑1 (MQ=255)                                                                                 tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGaa                                                                                                           >  2:264383/1‑74 (MQ=255)                                                                                 tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGaa                                                                                                           <  1:264383/74‑1 (MQ=255)                                                                                 tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCa                                                                                         <  1:675157/92‑1 (MQ=255)                                                                                 tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCa                                                                                         >  2:675157/1‑92 (MQ=255)                                                                                 tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                      <  2:204865/127‑1 (MQ=255)                                                                                 tGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTAc                                                      >  1:204865/1‑127 (MQ=255)                                                                                    gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       <  2:1201202/139‑1 (MQ=255)                                                                                    gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       <  2:1102996/139‑1 (MQ=255)                                                                                    gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       <  2:286176/139‑1 (MQ=255)                                                                                    gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       <  2:415976/139‑1 (MQ=255)                                                                                    gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       <  2:1095333/139‑1 (MQ=255)                                                                                    gTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       <  1:642943/139‑1 (MQ=255)                                                                                           tGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCa                                <  2:396261/139‑1 (MQ=255)                                                                                                      tGGTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATActgtctgt                     <  2:982788/139‑1 (MQ=255)                                                                                                       ggTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACagtag                                                 <  1:605809/110‑1 (MQ=255)                                                                                                       ggTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACagtag                                                 >  2:605809/1‑110 (MQ=255)                                                                                                        gTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATcaacaa                                                                                                                <  1:231485/46‑1 (MQ=255)                                                                                                        gTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATcaacaa                                                                                                                >  2:231485/1‑46 (MQ=255)                                                                                                        gTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAccc                                                              >  2:1041657/1‑96 (MQ=255)                                                                                                        gTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTAccc                                                              <  1:1041657/96‑1 (MQ=255)                                                                                                        gTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTg                   <  2:955591/139‑1 (MQ=255)                                                                                                        gTCTTGGCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTg                   <  2:121505/139‑1 (MQ=255)                                                                                                             ggCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGaa                                                                                                           <  2:282013/46‑1 (MQ=255)                                                                                                             ggCATTGATTTCAGAAATATATATATCACCAAAATCAACAACAGaa                                                                                                           >  1:282013/1‑46 (MQ=255)                                                                                                                  tGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       <  2:1010250/109‑1 (MQ=255)                                                                                                                  tGATTTCAGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGt                                       >  1:1010250/1‑109 (MQ=255)                                                                                                                         aGAAATATATATATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTGCCGCATAAGAGAGTACg  >  2:640656/1‑139 (MQ=255)                                                                                                                                tatataTCACCAAAATCAACAACAGAAATGTCCTGAt                                                                                                 <  1:951232/37‑1 (MQ=255)                                                                                                                                tatataTCACCAAAATCAACAACAGAAATGTCCTGAt                                                                                                 >  2:951232/1‑37 (MQ=255)                                                                                                                                    |                                                                                                                                AGCAGTTCCCTCGCATGTGGCTCGCTTGGTTAATTTTATTTGCGCTTTTTTATTGGGGATACCCGCACAGTCTTTGAATTTGAGTTTGAATGTTTTTACTTTGGTCTTGGCATTGATTTCAGAAATATATACATCACCAAAATCAACAACAGAAATGTCCTGATCAGAGTCATTTACTAATGTAGCGGTACAAGTACCCATTTCTACAGTAGTAATAACGGTCAGTCCAATACTGTCTGTTGCCGCATAAGAGAGTACG  >  NC_000913/151905‑152163

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