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%	101.4 / NA	28	Q56R (CAA→CGA)	ybhC	acyl‑CoA thioesterase, lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/16);  total (12/16)

CGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCA  >  NC_000913/806973‑807256                                                                                                                                               |                                                                                                                                              cGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCg                                                                                                                                       <  1:543220/151‑1 (MQ=255)   cATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGaa                                                                                                                                     <  1:608289/151‑1 (MQ=255)          cgccgcCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTgc                                                                                                                              <  2:442978/151‑1 (MQ=255)             cgccTGAATGGTGGTATGCGTTACTCCTTGAGTGGCCGCCGTGCCGACAAGAAAGTCAGGTTGCTCAGGCAGGTTAATCGGGTAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTgcgc                                                                                                                            <  2:439084/150‑1 (MQ=255)                              gCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGaa                                                                                                           <  1:465782/150‑1 (MQ=255)                               cGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATc                                                                                                         <  2:58246/151‑1 (MQ=255)                                         tGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCgaga                                                                                                 <  1:611413/149‑1 (MQ=255)                                          gAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTgcgc                                                                                                                            <  1:423287/121‑1 (MQ=255)                                          gAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTgcgc                                                                                                                            >  2:423280/1‑121 (MQ=255)                                               cccGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTAc                                                                                         >  1:387860/1‑151 (MQ=255)                                                     gggCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGgcg                                                                                   <  1:608128/151‑1 (MQ=255)                                                         cGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCgg                                                                                <  1:125500/150‑1 (MQ=255)                                                         cGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCgg                                                                                >  2:125497/1‑150 (MQ=255)                                                             aacaaAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTg                                                                            >  2:233398/1‑150 (MQ=255)                                                                  aGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATc                                                                         >  2:317251/1‑148 (MQ=255)                                                                                          tCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTg                                               >  2:471268/1‑150 (MQ=255)                                                                                             gggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCt                                           <  1:58673/151‑1 (MQ=255)                                                                                              ggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTa                                          <  2:54465/151‑1 (MQ=255)                                                                                              ggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTa                                          <  2:76024/151‑1 (MQ=255)                                                                                                   gATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAgg                                      >  2:452087/1‑150 (MQ=255)                                                                                                       ccACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTc                                 <  2:63235/151‑1 (MQ=255)                                                                                                        cACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCa                                <  2:298634/151‑1 (MQ=255)                                                                                                               cagcaCCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTc                           >  2:189871/1‑149 (MQ=255)                                                                                                                 gcaCCTGGTGTCAGGGATGCAAAATAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGc                       >  1:594565/1‑151 (MQ=255)                                                                                                                                     aaaaTAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGcc   >  1:445036/1‑151 (MQ=255)                                                                                                                                      aaaTAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCa  <  1:420590/151‑1 (MQ=255)                                                                                                                                        aTAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCCAAAGCCAATGCCAGCGcc   >  1:15409/1‑148 (MQ=255)                                                                                                                                        aTAGTGTCGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGcc   >  1:520087/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                              CGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCA  >  NC_000913/806973‑807256

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