Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,720,213	A→C	E87D (GAA→GAC)	patZ →	peptidyl‑lysine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,720,213	0	A	C	100.0%	128.3 / NA	37	E87D (GAA→GAC)	patZ	peptidyl‑lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/16);  total (21/16)

GATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTG  >  NC_000913/2720136‑2720290                                                                              |                                                                              gATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAcgag                                                                            <  2:486236/81‑1 (MQ=255) gATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAcgag                                                                            <  2:486235/81‑1 (MQ=255)   tatTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCt                                                                          >  1:486301/1‑81 (MQ=255)     ttGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTc                                                                         <  2:486237/80‑1 (MQ=255)      tGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCgg                                                                       >  1:486278/1‑81 (MQ=255)        ccAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGc                                                                      <  2:486238/80‑1 (MQ=255)         cAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGc                                                                      >  1:486275/1‑79 (MQ=255)          aGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCgag                                                                   <  1:486239/81‑1 (MQ=255)            cTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCgaga                                                                  >  1:486265/1‑80 (MQ=255)               gCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGaaa                                                                <  1:486240/79‑1 (MQ=255)               gCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAgg                                                              >  2:486287/1‑81 (MQ=255)                    ttACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGt                                                          >  1:486284/1‑80 (MQ=255)                     tACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGt                                                          <  2:486242/79‑1 (MQ=255)                     tACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGt                                                          <  2:486241/79‑1 (MQ=255)                       cacCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTaaaa                                                      <  2:486243/81‑1 (MQ=255)                        acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAAc                                                     >  2:486280/1‑81 (MQ=255)                            gACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTg                                                  >  2:486302/1‑80 (MQ=255)                               cTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCat                                               >  2:486309/1‑80 (MQ=255)                                      ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                       <  1:486245/81‑1 (MQ=255)                                      ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                       <  1:486244/81‑1 (MQ=255)                                            gTACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGccc                                 >  2:486288/1‑81 (MQ=255)                                             tACCAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGccc                                 <  2:486246/80‑1 (MQ=255)                                                cAATGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGc                             <  2:486247/81‑1 (MQ=255)                                                 aaTGCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCa                            >  2:486286/1‑81 (MQ=255)                                                    gCCAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCg                         >  1:486285/1‑81 (MQ=255)                                                     ccAGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCg                         <  2:486248/80‑1 (MQ=255)                                                       aGCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCaa                      >  2:486290/1‑81 (MQ=255)                                                        gCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCaa                      >  1:486303/1‑80 (MQ=255)                                                        gCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCaa                      >  1:486282/1‑80 (MQ=255)                                                        gCCGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAac                     <  2:486249/81‑1 (MQ=255)                                                          cGTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAacac                   <  2:486250/81‑1 (MQ=255)                                                           gTAATCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAacac                   >  2:486292/1‑80 (MQ=255)                                                             aaTCTTGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACg                  >  1:486293/1‑79 (MQ=255)                                                                 ttGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtc            >  1:486305/1‑81 (MQ=255)                                                                  tGCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtct           >  2:486289/1‑81 (MQ=255)                                                                   gCTCTTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtc            <  1:486251/79‑1 (MQ=255)                                                                    ctctTCTGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtctc          <  1:486252/80‑1 (MQ=255)                                                                      cttctGGACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCcgc       >  1:486291/1‑81 (MQ=255)                                                                           ggACGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctg  >  2:486304/1‑81 (MQ=255)                                                                              |                                                                              GATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTG  >  NC_000913/2720136‑2720290

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