Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,485,906	C→A	A16S (GCC→TCC)	pepA ←	aminopeptidase A/I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,485,906	0	C	A	100.0%	86.7 / NA	22	A16S (GCC→TCC)	pepA	aminopeptidase A/I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/9);  total (13/9)

TCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485861‑4485951                                              |                                              tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                            <  1:1434602/49‑1 (MQ=255)   ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                          >  1:1700930/1‑49 (MQ=255)       aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                      >  1:350972/1‑49 (MQ=255)       aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                      >  1:1548773/1‑49 (MQ=255)       aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                      <  1:1662994/49‑1 (MQ=255)          ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                   <  1:381191/49‑1 (MQ=255)             gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                <  1:1986623/49‑1 (MQ=255)              aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                               >  1:1905940/1‑49 (MQ=255)                gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                             >  1:1750724/1‑49 (MQ=255)                   gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                           <  1:1137492/48‑1 (MQ=255)                      cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                       >  1:1748903/1‑49 (MQ=255)                       gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                      >  1:2100917/1‑49 (MQ=255)                        aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                     <  1:2074204/49‑1 (MQ=255)                        aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                     >  1:1193346/1‑49 (MQ=255)                         aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                    >  1:1727792/1‑49 (MQ=255)                          gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                   <  1:112823/49‑1 (MQ=255)                             gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                <  1:919641/49‑1 (MQ=255)                              cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa               <  1:906880/49‑1 (MQ=255)                                 acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac            >  1:1575504/1‑49 (MQ=255)                                     cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa        >  1:1261092/1‑49 (MQ=255)                                         gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc     >  1:294396/1‑48 (MQ=255)                                           aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt  >  1:535754/1‑49 (MQ=255)                                              |                                              TCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485861‑4485951

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