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%	112.3 / NA	31	A16S (GCC→TCC)	pepA	aminopeptidase A/I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/15);  total (16/15)

ATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485860‑4485951                                               |                                              aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                            >  1:1355062/1‑50 (MQ=255)  tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                           <  1:2326033/50‑1 (MQ=255)    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                         >  1:1006029/1‑50 (MQ=255)     gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                        >  1:619725/1‑50 (MQ=255)      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                       <  1:1932253/50‑1 (MQ=255)      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                       <  1:1287515/50‑1 (MQ=255)        aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                     >  1:1268018/1‑50 (MQ=255)         aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                    >  1:1006406/1‑50 (MQ=255)         aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                    >  1:2042311/1‑50 (MQ=255)          aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                   <  1:1730505/50‑1 (MQ=255)            gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                 <  1:1799033/50‑1 (MQ=255)             cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                >  1:1123321/1‑50 (MQ=255)              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                               >  1:561893/1‑50 (MQ=255)               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                              >  1:2176177/1‑50 (MQ=255)                  tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                           <  1:923259/50‑1 (MQ=255)                  tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                           <  1:1752872/50‑1 (MQ=255)                   ggTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                           >  1:358694/1‑49 (MQ=255)                   ggTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                           >  1:45063/1‑49 (MQ=255)                         aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                    >  1:1034383/1‑50 (MQ=255)                           gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                  <  1:441196/50‑1 (MQ=255)                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                 >  1:1608622/1‑50 (MQ=255)                             cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                <  1:1462804/50‑1 (MQ=255)                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa               <  1:530516/50‑1 (MQ=255)                                 cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac            <  1:1317683/50‑1 (MQ=255)                                   cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt           >  1:1116411/1‑49 (MQ=255)                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa        <  1:1294421/50‑1 (MQ=255)                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa        <  1:849412/50‑1 (MQ=255)                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt      >  1:1817616/1‑50 (MQ=255)                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt      >  1:1920291/1‑50 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa   <  1:1169200/49‑1 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt  <  1:2088450/50‑1 (MQ=255)                                               |                                              ATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485860‑4485951

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