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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952                                                |                                               aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:3001946/1‑50 (MQ=255) aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:3026507/1‑50 (MQ=255) aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              <  1:3428733/50‑1 (MQ=255)  aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             <  1:3898552/50‑1 (MQ=255)   tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            <  1:1346433/50‑1 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        <  1:3978644/50‑1 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        <  1:3079833/50‑1 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:3617507/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  1:2824687/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  1:375697/1‑50 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    >  1:1771326/1‑50 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    <  1:2603771/50‑1 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    >  1:1118856/1‑50 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    >  1:2026583/1‑50 (MQ=255)             gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  >  1:3949325/1‑50 (MQ=255)              cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                 >  1:1812930/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:3924016/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:2476195/1‑50 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:2562637/50‑1 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:2164217/50‑1 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              >  1:2384175/1‑50 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:1489132/1‑50 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             <  1:1785369/50‑1 (MQ=255)                    ggTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            >  1:4410255/1‑49 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:427111/49‑1 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         <  1:287272/50‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:3978696/49‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:4269129/49‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:1355503/49‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:1810910/49‑1 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       <  1:394063/49‑1 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       <  1:1638944/49‑1 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     >  1:1598454/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  1:666993/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  1:1915384/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  1:2935132/1‑50 (MQ=255)                            gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   <  1:3078230/50‑1 (MQ=255)                              cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:2819379/50‑1 (MQ=255)                               gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:1380201/1‑50 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              <  1:506856/50‑1 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:2896997/1‑49 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            <  1:1936688/50‑1 (MQ=255)                                    cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            >  1:66320/1‑49 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:2789860/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:1789656/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:1267656/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        <  1:3747258/50‑1 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:3356894/1‑50 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  1:95062/50‑1 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      <  1:2487919/50‑1 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      <  1:843116/50‑1 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:2745167/1‑49 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:1473845/1‑49 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:1047079/1‑49 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:4237612/1‑50 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:2465498/50‑1 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:1103746/50‑1 (MQ=255)                                                |                                               AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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