Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	806,505	G→A	S259S (AGC→AGT)	ybhC ←	acyl‑CoA thioesterase, lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	806,505	0	G	A	100.0%	174.9 / NA	51	S259S (AGC→AGT)	ybhC	acyl‑CoA thioesterase, lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/25);  total (26/25)

TCACCACGCGGAATTCGGTGTTATCGAACACCACTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCGCTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCATCTACGCTATCGCCCAGCGTGTTTTC  >  NC_000913/806372‑806636                                                                                                                                       |                                                                                                                                    tCACCACGCGGAATTCGGTGTTATCGAACACCACTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGtt                                                                                                                                <  2:1336248/139‑1 (MQ=255)    ccacGCGGAATTCGGTGTTATCGAACACCACTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GTTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGt                                                                                                                             >  2:1386268/1‑139 (MQ=255)    ccacGCGGAATTCGGTGTAATCGAACACCACTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGt                                                                                                                             <  2:313758/139‑1 (MQ=255)                            acacCACTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GTTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGAcc                                                                                                     <  1:1386268/139‑1 (MQ=255)                             caccacTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGa                                                                                                                           >  1:662325/1‑116 (MQ=255)                             caccacTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGa                                                                                                                           <  2:662325/116‑1 (MQ=255)                              accacTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCga                                                                                                   <  1:1046508/139‑1 (MQ=255)                                  cTGCGCCGCGACCAGAACCGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAAt                                                                                               >  1:943888/1‑139 (MQ=255)                                  cTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTATGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAAt                                                                                               >  2:516260/1‑139 (MQ=255)                                  cTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAAt                                                                                               >  1:586449/1‑139 (MQ=255)                                       ccgcGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTaa                                                                                          <  1:230502/139‑1 (MQ=255)                                        cgcgACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAAc                                                                                         >  1:1453039/1‑139 (MQ=255)                                         gcgACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACg                                                                                        <  2:943888/139‑1 (MQ=255)                                          cgACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTTGGGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACg                                                                                        >  2:236502/1‑139 (MQ=255)                                          cgACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACgt                                                                                       >  2:350132/1‑139 (MQ=255)                                          cgACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACgt                                                                                       >  2:828820/1‑139 (MQ=255)                                             ccAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTATGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACgttgt                                                                                    <  1:516260/139‑1 (MQ=255)                                             ccAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACgttgt                                                                                    <  1:350132/139‑1 (MQ=255)                                             ccAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACgttgt                                                                                    <  2:586449/139‑1 (MQ=255)                                             ccAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACgttgt                                                                                    <  1:710613/139‑1 (MQ=255)                                             ccAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACgttgt                                                                                    <  2:1021710/139‑1 (MQ=255)                                                       tatCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTAc                                                                          >  1:508211/1‑139 (MQ=255)                                                       tatCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTAc                                                                          >  1:1110892/1‑139 (MQ=255)                                                        atCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTAcc                                                                         >  2:791870/1‑139 (MQ=255)                                                             caTCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTc                                                                    <  1:828820/139‑1 (MQ=255)                                                               tCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCAc                                                                  <  1:1005986/139‑1 (MQ=255)                                                                   cTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATc                                                              >  2:968397/1‑139 (MQ=255)                                                                   cTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATc                                                              >  1:482992/1‑139 (MQ=255)                                                                    ttCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCa                                                             >  2:185829/1‑139 (MQ=255)                                                                      cAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCAc                                                                  <  2:1033770/132‑1 (MQ=255)                                                                      cAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCAc                                                                  >  1:1033770/1‑132 (MQ=255)                                                                                    ggTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCAcc                                              >  1:1318071/1‑139 (MQ=255)                                                                                       cACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGt                                                                     <  1:4872/113‑1 (MQ=255)                                                                                       cACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGt                                                                     >  2:4872/1‑113 (MQ=255)                                                                                         ccAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCgg                                         >  1:1510785/1‑139 (MQ=255)                                                                                          cAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTaa                                                                                          <  1:450604/89‑1 (MQ=255)                                                                                          cAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTaa                                                                                          >  2:450604/1‑89 (MQ=255)                                                                                               tACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTAcc                                                                         >  1:375062/1‑101 (MQ=255)                                                                                               tACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTAcc                                                                         <  2:375062/101‑1 (MQ=255)                                                                                                 cgcgGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGt                                                                                                                             <  2:1256557/47‑1 (MQ=255)                                                                                                 cgcgGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGt                                                                                                                             >  1:1256557/1‑47 (MQ=255)                                                                                                  gcgGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTAcc                                <  2:1510785/139‑1 (MQ=255)                                                                                                  gcgGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTAcc                                <  2:1507688/139‑1 (MQ=255)                                                                                                   cgGCTGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCt                               <  1:791870/139‑1 (MQ=255)                                                                                                       tGACGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCAt                           >  2:301804/1‑139 (MQ=255)                                                                                                          cGATTCGTTTCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCAGATGGTTACCTGCATCTa                        <  2:1453039/139‑1 (MQ=255)                                                                                                                   tCCAGACGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCATCTACGCTATCGc               >  2:20035/1‑139 (MQ=255)                                                                                                                        aCGGTTCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCATCTACGCTATCGCCCAGc          >  1:196574/1‑139 (MQ=255)                                                                                                                            ttCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCATCTACGCTATCGCCCAGCgtgt      >  1:831912/1‑139 (MQ=255)                                                                                                                            ttCTGCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCATCTACGCTATCGCCCAGCgtgt      <  2:831912/139‑1 (MQ=255)                                                                                                                                gCACACCACTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCATCTACGCTATCGCCCAGCGTGTTTTc  <  2:482992/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    TCACCACGCGGAATTCGGTGTTATCGAACACCACTGCGCCGCGACCAGAAACGATATCCACATCCCCTTCAATGTAGCTGTT‑GGTCACCAGCGTACGCGGCTGACGATTCGTTTCCAGACGGTTCTGCACACCGCTGTTGGTGACAAAGAAGGTGTTCTGACGACCGAGAATGTTAACGTTGTTAATCTGTACCTGGTCACCATCAGTACGCAGTGCCACCGCCGGATGGTTACCTGCATCTACGCTATCGCCCAGCGTGTTTTC  >  NC_000913/806372‑806636

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: ‑