Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	807,115	T→C	Q56R (CAA→CGA)	ybhC ←	acyl‑CoA thioesterase, lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	807,115	0	T	C	100.0%	92.6 / NA	26	Q56R (CAA→CGA)	ybhC	acyl‑CoA thioesterase, lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/15);  total (11/15)

GCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCG  >  NC_000913/806971‑807253                                                                                                                                                 |                                                                                                                                           gCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCAt                                                                                                                                      <  1:265181/151‑1 (MQ=255)   cGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCg                                                                                                                                    >  1:527597/1‑151 (MQ=255)          aCCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTc                                                                                                                              <  1:192400/150‑1 (MQ=255)                gccTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGctt                                                                                                                       >  2:239316/1‑151 (MQ=255)                                     acacCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGc                                                                                                  <  2:245977/151‑1 (MQ=255)                                                gCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGt                                                                                        <  2:125752/150‑1 (MQ=255)                                                gCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGt                                                                                        <  1:586933/150‑1 (MQ=255)                                                gCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTa                                                                                       <  2:529723/151‑1 (MQ=255)                                                 cccGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGTTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTAc                                                                                      <  1:440067/151‑1 (MQ=255)                                                                acaaAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGAt                                                                       <  2:563809/151‑1 (MQ=255)                                                                       cAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAgg                                                                >  2:642073/1‑151 (MQ=255)                                                                       cAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATGCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTTCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAgg                                                                <  1:171211/151‑1 (MQ=255)                                                                          gTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACg                                                             <  1:290980/151‑1 (MQ=255)                                                                          gTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACg                                                             >  2:449340/1‑151 (MQ=255)                                                                                   gcaggGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATcggg                                                    <  1:44376/151‑1 (MQ=255)                                                                                   gcaggGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATcggg                                                    >  1:290515/1‑151 (MQ=255)                                                                                    caggGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATcgggc                                                   <  2:290516/151‑1 (MQ=255)                                                                                                  aaGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACa                                     <  2:563565/151‑1 (MQ=255)                                                                                                    ggATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAgg                                   >  2:425085/1‑151 (MQ=255)                                                                                                     gATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGTCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAgg                                   >  1:336317/1‑150 (MQ=255)                                                                                                      aTTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCg                                 >  2:285544/1‑151 (MQ=255)                                                                                                             gCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCg                          >  1:311594/1‑151 (MQ=255)                                                                                                             gCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGCCa                             >  2:296286/1‑148 (MQ=255)                                                                                                                         ggTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCaaaa               <  1:119771/150‑1 (MQ=255)                                                                                                                                     gCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAgc   <  1:642063/150‑1 (MQ=255)                                                                                                                                      cAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAgcg  >  1:192218/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                           GCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCG  >  NC_000913/806971‑807253

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