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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952                                                |                                               aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              <  1:1203942/50‑1 (MQ=255)  aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             >  1:1701007/1‑50 (MQ=255)  aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             >  1:2459303/1‑50 (MQ=255)   tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            >  1:2341188/1‑50 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:433797/49‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:1314502/49‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:4981399/1‑50 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:3796920/1‑50 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:2569047/1‑50 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:5953731/1‑50 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:5907131/1‑50 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:2995067/1‑50 (MQ=255)        gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  1:1905856/50‑1 (MQ=255)        gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  1:97469/50‑1 (MQ=255)            ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   >  1:4851980/1‑50 (MQ=255)             gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  >  1:2251668/1‑50 (MQ=255)             gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  <  1:4167044/50‑1 (MQ=255)             gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  >  1:1806020/1‑50 (MQ=255)              cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                 >  1:4949379/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:2241274/50‑1 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:5927535/50‑1 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:5367181/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:5771830/1‑50 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              >  1:4489404/1‑50 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:4291429/50‑1 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:1687408/1‑50 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:836481/1‑50 (MQ=255)                   tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            <  1:3049959/50‑1 (MQ=255)                   tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            <  1:5196045/50‑1 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:3036186/49‑1 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTa                          >  1:2282604/1‑50 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         >  1:2198594/1‑50 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         <  1:4261476/50‑1 (MQ=255)                       tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:5107119/1‑50 (MQ=255)                       tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:2865921/1‑50 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:5802006/49‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:5215359/1‑50 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:4286290/1‑50 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:596459/1‑50 (MQ=255)                           aGACGCCCACGGCAATGCAGGAACGCCTCTGTTTCTCCGGGCTACCGCtt                    <  1:1047225/50‑1 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  1:2756787/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    <  1:4858715/50‑1 (MQ=255)                            gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   <  1:1064185/50‑1 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:4075978/50‑1 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:6230534/50‑1 (MQ=255)                              cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:3549310/50‑1 (MQ=255)                              cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:576647/50‑1 (MQ=255)                                cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac               <  1:6030930/50‑1 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              >  1:3973695/1‑50 (MQ=255)                                  cACGACGATGCAGGAACTCCGCTGTTTCTACGGGCTACTGCTTTTTacac             >  1:3379073/1‑50 (MQ=38)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:4375855/1‑49 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:4064174/1‑49 (MQ=255)                                    cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            >  1:2460293/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:4182496/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:4066437/1‑49 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         <  1:5678706/50‑1 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:4675485/1‑50 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  1:1897067/50‑1 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  1:1360505/1‑50 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      <  1:3404592/50‑1 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:1286770/1‑50 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:4120460/1‑50 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:6341348/1‑50 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc     <  1:876368/50‑1 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:987634/1‑49 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:2016429/1‑49 (MQ=255)                                           gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    >  1:4307481/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:3435511/50‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:3687432/50‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:2891189/50‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:5264453/50‑1 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:2413285/1‑50 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:2513169/1‑50 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:5308883/1‑50 (MQ=255)                                                |                                               AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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