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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952                                                |                                               aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:4601741/1‑50 (MQ=255) aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:412246/1‑50 (MQ=255)  aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             <  1:2094694/50‑1 (MQ=255)  aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             >  1:84450/1‑50 (MQ=255)   tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            <  1:1151327/50‑1 (MQ=255)   tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            >  1:4691992/1‑50 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:1188401/49‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          <  1:525552/50‑1 (MQ=255)     ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:1875428/1‑50 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:4777776/1‑50 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:1537254/1‑50 (MQ=255)       agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        <  1:3618138/50‑1 (MQ=255)        gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  1:1248656/50‑1 (MQ=255)        gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       >  1:4327958/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  1:97988/1‑50 (MQ=255)         aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      <  1:3040287/50‑1 (MQ=255)          aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                     >  1:3277064/1‑50 (MQ=255)           aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    <  1:3878394/50‑1 (MQ=255)            ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   <  1:3039180/50‑1 (MQ=255)             gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  >  1:3930935/1‑50 (MQ=255)             gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  <  1:2084152/50‑1 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:4252605/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:1230761/1‑50 (MQ=255)               gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:1159698/1‑50 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:3255102/50‑1 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:1099518/1‑50 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:2495532/50‑1 (MQ=255)                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:1400770/50‑1 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              >  1:1817461/1‑50 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:1752853/50‑1 (MQ=255)                 cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              >  1:3247945/1‑50 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             <  1:4489727/50‑1 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             <  1:3190028/50‑1 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             <  1:1418857/50‑1 (MQ=255)                  gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:2054253/1‑50 (MQ=255)                   tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            <  1:3513924/50‑1 (MQ=255)                     gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:2797510/49‑1 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         <  1:4489417/50‑1 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         >  1:3070606/1‑50 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         >  1:3573084/1‑50 (MQ=255)                      ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         <  1:4669051/50‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:1965468/49‑1 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:31483/1‑50 (MQ=255)                        cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:419930/1‑50 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       <  1:1517282/49‑1 (MQ=255)                         gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:1019740/1‑50 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     <  1:1225573/50‑1 (MQ=255)                          aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     <  1:2202957/50‑1 (MQ=255)                              cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:2536223/50‑1 (MQ=255)                              cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:2457413/50‑1 (MQ=255)                               gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                <  1:2104879/50‑1 (MQ=255)                                cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac               <  1:343231/50‑1 (MQ=255)                                cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:4359234/1‑49 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              <  1:1705088/50‑1 (MQ=255)                                 ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              >  1:4558744/1‑50 (MQ=255)                                  cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:5096621/1‑50 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            <  1:1560548/50‑1 (MQ=255)                                   acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            <  1:4712852/50‑1 (MQ=255)                                    cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            >  1:3269127/1‑49 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:3821941/50‑1 (MQ=255)                                     gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:4643661/50‑1 (MQ=255)                                      acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:5103389/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:4382234/1‑50 (MQ=255)                                       cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        <  1:5036545/50‑1 (MQ=255)                                        gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  1:3960292/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    <  1:1770841/49‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    <  1:874137/49‑1 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:1418207/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:3580936/1‑50 (MQ=255)                                            cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:2524733/1‑50 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:2575373/50‑1 (MQ=255)                                             aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:1206299/50‑1 (MQ=255)                                                |                                               AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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: ‑