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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485859‑4485951                                                |                                              aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                             >  1:1223857/1‑50 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                       <  1:1199971/50‑1 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                     <  1:282507/50‑1 (MQ=255)          aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                    >  1:2126336/1‑50 (MQ=255)            ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                  <  1:589846/50‑1 (MQ=255)              cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                <  1:2258294/50‑1 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                               >  1:2388720/1‑50 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                              >  1:2687193/1‑50 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                             >  1:1666992/1‑50 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                            >  1:535999/1‑50 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                          <  1:44154/49‑1 (MQ=255)                      ttCGAAGACGCCCACGACGATGGAGGAACTCCGCTGTTTCTCCGGGCTAc                        >  1:2297056/1‑50 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                        >  1:153234/1‑50 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                        <  1:2084020/50‑1 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                        >  1:454262/1‑50 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                        <  1:745000/50‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                       <  1:1721002/49‑1 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                    >  1:569568/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                   >  1:1117979/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                   >  1:896079/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                   >  1:1129484/1‑50 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                 >  1:2159670/1‑50 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                 <  1:1950256/50‑1 (MQ=255)                                cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac              <  1:2590353/50‑1 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca             <  1:875834/50‑1 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca             <  1:1403397/50‑1 (MQ=255)                                    cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa          <  1:1692522/50‑1 (MQ=255)                                    cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa          <  1:347719/50‑1 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa         <  1:1628828/50‑1 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa          >  1:2328507/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa          >  1:484872/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa          >  1:757545/1‑49 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa        <  1:1154011/50‑1 (MQ=255)                                       cgaTGCAGGAACTTCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc       <  1:112756/50‑1 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc     >  1:357504/1‑50 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc     >  1:1430179/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt  <  1:904601/50‑1 (MQ=255)                                                |                                              AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485859‑4485951

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