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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952                                                |                                               aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:2697433/1‑50 (MQ=255) aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:5418379/1‑50 (MQ=255)  aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             >  1:5833596/1‑50 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:3411279/49‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          <  1:5613295/50‑1 (MQ=255)      gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         >  1:1181214/1‑50 (MQ=255)      gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         >  1:3137549/1‑50 (MQ=255)      gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         >  1:2375906/1‑50 (MQ=255)      gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         >  1:1955603/1‑50 (MQ=255)        gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  1:1495127/50‑1 (MQ=255)        gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  1:3409540/50‑1 (MQ=255)        gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       >  1:3001638/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      <  1:5058270/50‑1 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  1:3804637/1‑50 (MQ=255)          aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                     >  1:1836332/1‑50 (MQ=255)          aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                     >  1:482074/1‑50 (MQ=255)          aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                     >  1:3035051/1‑50 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    <  1:4415937/50‑1 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    <  1:5700067/50‑1 (MQ=255)              cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                 >  1:4628569/1‑50 (MQ=255)              cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                 >  1:1415717/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:3388680/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:5798996/50‑1 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:3452716/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:1270529/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:1398548/1‑50 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:1061137/50‑1 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:2410670/1‑50 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:2741560/50‑1 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:2362172/50‑1 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:1732284/1‑50 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:249473/1‑50 (MQ=255)                   tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            >  1:1305813/1‑50 (MQ=255)                   tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            <  1:3045534/50‑1 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:4836295/49‑1 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTa                          <  1:2111405/50‑1 (MQ=255)                       tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:2147485/1‑50 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:2148681/49‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:2708517/49‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:5328136/49‑1 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       <  1:4169588/49‑1 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:4640433/1‑50 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:4712213/1‑50 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     <  1:3364197/50‑1 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     >  1:4939349/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  1:4111128/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    <  1:1947188/50‑1 (MQ=255)                            gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   <  1:1307947/50‑1 (MQ=255)                            gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   >  1:5683543/1‑50 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:3842562/50‑1 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:3359494/50‑1 (MQ=255)                               gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:2933876/1‑50 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              <  1:4001053/50‑1 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              >  1:5378370/1‑50 (MQ=255)                                  cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             <  1:1798553/50‑1 (MQ=255)                                  cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:4675264/1‑50 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:3654227/1‑49 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:6006863/1‑49 (MQ=255)                                    cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           <  1:3748590/50‑1 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:4783006/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:4816585/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:842932/1‑49 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:302484/1‑50 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:4824517/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:5239395/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:246791/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:4006519/1‑50 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  1:5060835/1‑50 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTAACGCTTTTTACACTAAACt       >  1:4336585/1‑50 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:1989794/1‑50 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:5413430/1‑50 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc     <  1:3028503/50‑1 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:3622212/1‑50 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:997550/1‑50 (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 < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑