Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,167,586	T→C	T45A (ACG→GCG)	ygiS ←	putative deoxycholate binding periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,167,586	0	T	C	100.0%	158.1 / NA	45	T45A (ACG→GCG)	ygiS	putative deoxycholate binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/23);  total (22/23)

GCCGTCCAGTATCTCCCAGCGTTCAGCCTGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGTACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGTATACA  >  NC_000913/3167452‑3167717                                                                                                                                       |                                                                                                                                    gccgTCCAGTATCTCCCAGCGTTCAGCCTGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTg                                                                                                                                <  2:55149/140‑1 (MQ=255)               cccAGCGTTCAGCCTGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGAttg                                                                                                                  >  1:165949/1‑140 (MQ=255)               cccAGCGTTCAGCCTGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGAttg                                                                                                                  >  2:40418/1‑140 (MQ=255)                      ttCAGCCTGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAAc                                                                                                           >  1:716768/1‑140 (MQ=255)                        cAGCCTGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGa                                                                                                         <  2:165949/140‑1 (MQ=255)                             tGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGAc                                                                                                    <  1:176999/140‑1 (MQ=255)                              gAGCGGGCTGCACCTGG‑CTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACtt                                                                                                  >  1:825363/1‑140 (MQ=255)                                        cACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCgggg                                                                                         <  2:716768/140‑1 (MQ=255)                                             ggCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGGCTTGTTGCGGGGCGAGc                                                                                    >  2:769224/1‑140 (MQ=255)                                                   cACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGgtgtg                                                                              >  2:599235/1‑140 (MQ=255)                                                         ccatccatACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCg                                                                        <  1:5282/140‑1 (MQ=255)                                                          catccatACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCgg                                                                       >  2:656066/1‑140 (MQ=255)                                                                tACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGt                                                                 >  2:355589/1‑140 (MQ=255)                                                                   cAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTcag                                                              <  2:207473/140‑1 (MQ=255)                                                                   cAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTcag                                                              <  2:684538/140‑1 (MQ=255)                                                                   cAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTcag                                                              >  1:684538/1‑140 (MQ=255)                                                                           cAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTag                                                      >  2:341564/1‑140 (MQ=255)                                                                             aaCAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTagag                                                    <  2:455213/140‑1 (MQ=255)                                                                                aGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAgagg                                                 >  2:838776/1‑140 (MQ=255)                                                                                     cAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCg                                            <  1:100213/140‑1 (MQ=255)                                                                                        cacaATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAgaggag                                               <  1:704521/134‑1 (MQ=255)                                                                                        cacaATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAgaggag                                               >  2:704521/1‑134 (MQ=255)                                                                                            aTCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTAc                                     >  1:731003/1‑140 (MQ=255)                                                                                            aTCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTAc                                     >  1:354308/1‑140 (MQ=255)                                                                                                gcgcGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAgg                                 <  1:769224/140‑1 (MQ=255)                                                                                                 cgcgGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGc                                >  2:853601/1‑140 (MQ=255)                                                                                                 cgcgGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGc                                >  1:623298/1‑140 (MQ=255)                                                                                                  gcgGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTa                                      <  2:223912/133‑1 (MQ=255)                                                                                                  gcgGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTa                                      >  1:223912/1‑133 (MQ=255)                                                                                                              tcctccACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGAc                            <  2:744502/131‑1 (MQ=255)                                                                                                              tcctccACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGAc                            >  1:744502/1‑131 (MQ=255)                                                                                                                ctccACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCataata                 >  2:291199/1‑140 (MQ=255)                                                                                                                ctccACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCataata                 >  1:605986/1‑140 (MQ=255)                                                                                                                   cACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGGGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAAt              <  2:623298/140‑1 (MQ=255)                                                                                                                   cACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAAt              <  2:354308/140‑1 (MQ=255)                                                                                                                   cACCTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACAAGGCTGACCAGCCATAATAAAt              >  2:647487/1‑140 (MQ=255)                                                                                                                     ccTTTTGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGGACGTCAGCAGCGTAAGGGGGAGCCGCACTTAACAGGCAGACCAGCCATAATAAAttt            >  2:199139/1‑140 (MQ=255)                                                                                                                       ttttGCGGGTCGAGCGCCCCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTcg          <  1:199139/140‑1 (MQ=255)                                                                                                                        tttGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTcgc         <  1:656066/140‑1 (MQ=255)                                                                                                                        tttGCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTcgc         <  1:291199/140‑1 (MQ=255)                                                                                                                           gCGGGTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGta      <  2:825363/140‑1 (MQ=255)                                                                                                                             gggTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGtata    <  1:599235/140‑1 (MQ=255)                                                                                                                             gggTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGtata    <  1:341564/140‑1 (MQ=255)                                                                                                                               gTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGTATACa  <  1:838776/140‑1 (MQ=255)                                                                                                                               gTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGTATACa  <  1:853601/140‑1 (MQ=255)                                                                                                                               gTCGAGCGCACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGTATACa  <  1:355589/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    GCCGTCCAGTATCTCCCAGCGTTCAGCCTGAGCGGGCTGCACCTGGCCTTCACCGTCCATCCATACCAGACCTTCAAACAGATCCAGCACAATCTGCGCGGCAGTATTCTCCTCCACCTTTTGCGGGTCGAGCGTACCTGGGTCGCTATGATTGTTGTAACGAAAGACTTGTTGCGGGGCGAGCGGTGTGTTGGCGGGAACGTCAGCAGCGTAGAGAGGAGCCGCACTTACCAGGCTGACCAGCCATAATAAATTTCGCGTATACA  >  NC_000913/3167452‑3167717

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