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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952                                                |                                               aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:2172794/1‑50 (MQ=255) aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:3137428/1‑50 (MQ=255)  aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             >  1:4671801/1‑50 (MQ=255)   tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            >  1:1176015/1‑50 (MQ=255)    cGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           >  1:319387/1‑50 (MQ=255)    cGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           >  1:465011/1‑50 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:1074076/49‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:4069496/49‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:5339370/49‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:4722706/49‑1 (MQ=255)      gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         <  1:1067050/50‑1 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        <  1:5586113/50‑1 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:2415946/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  1:3374737/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  1:333960/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  1:2221056/1‑50 (MQ=255)          aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                     <  1:5419037/50‑1 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    >  1:5458064/1‑50 (MQ=255)            ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   >  1:2288677/1‑50 (MQ=255)              cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                 <  1:4825962/50‑1 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:4539108/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:3748951/50‑1 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:3616015/1‑50 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:3797406/1‑50 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:5503907/50‑1 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:1938122/1‑50 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:17991/49‑1 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:2949728/49‑1 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:4833665/49‑1 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         <  1:1485361/50‑1 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         >  1:4880957/1‑50 (MQ=255)                       tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:4112460/1‑50 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:4822845/1‑50 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       <  1:4928309/49‑1 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       <  1:2321244/49‑1 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:3381736/1‑50 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     >  1:2994091/1‑50 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     >  1:5443430/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  1:3472619/1‑50 (MQ=255)                           aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    <  1:4240220/50‑1 (MQ=255)                            gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   >  1:1606860/1‑50 (MQ=255)                            gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   >  1:4283997/1‑50 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:4556519/50‑1 (MQ=255)                             aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  >  1:1109623/1‑50 (MQ=255)                              cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:813209/50‑1 (MQ=255)                              cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:769081/50‑1 (MQ=255)                               gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:1066393/1‑50 (MQ=255)                               gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:2916796/1‑50 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              >  1:4239803/1‑50 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              <  1:1876445/50‑1 (MQ=255)                                  cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:5316576/1‑50 (MQ=255)                                  cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:4825588/1‑50 (MQ=255)                                  cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:5425078/1‑50 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:5200631/1‑49 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:2447283/1‑49 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:1403799/1‑49 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            <  1:4757526/50‑1 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            <  1:2009094/50‑1 (MQ=255)                                    cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           <  1:1861747/50‑1 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:3403738/50‑1 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:5385906/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:4256859/1‑49 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:1155049/1‑50 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:2750931/1‑50 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         <  1:2483645/50‑1 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         <  1:4282973/50‑1 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  1:551365/1‑50 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  1:1316099/50‑1 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      <  1:2171981/50‑1 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      <  1:2932229/50‑1 (MQ=255)                                         aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:2136541/1‑50 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc     <  1:3620289/50‑1 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:3242951/1‑49 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:1504257/1‑49 (MQ=255)                                          tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:3789466/1‑49 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    <  1:4601445/49‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    <  1:3415985/49‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:1478329/50‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:4658450/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:673337/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:758653/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:3883989/50‑1 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:2165610/50‑1 (MQ=255)                                                |                                               AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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