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

ATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485860‑4485952                                               |                                               aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             <  1:926679/50‑1 (MQ=255)   cGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           >  1:1368310/1‑50 (MQ=255)    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:4586801/49‑1 (MQ=255)    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:3624229/49‑1 (MQ=255)    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:6203635/1‑50 (MQ=255)     gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         <  1:1978474/50‑1 (MQ=255)      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:2738027/1‑50 (MQ=255)      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        <  1:6861943/50‑1 (MQ=255)      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:6864024/1‑50 (MQ=255)       gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  1:2199266/50‑1 (MQ=255)        aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      <  1:1354683/50‑1 (MQ=255)          aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    >  1:1667474/1‑50 (MQ=255)          aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    >  1:803189/1‑50 (MQ=255)          aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    <  1:3572706/50‑1 (MQ=255)           ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   <  1:5726601/50‑1 (MQ=255)           ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   <  1:395720/50‑1 (MQ=255)            gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  <  1:3981540/50‑1 (MQ=255)            gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  <  1:6043028/50‑1 (MQ=255)             cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                 >  1:77537/1‑50 (MQ=255)              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:221212/50‑1 (MQ=255)              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:1341573/1‑50 (MQ=255)               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:4976017/50‑1 (MQ=255)               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:1891984/1‑50 (MQ=255)                cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:2503509/50‑1 (MQ=255)                 gTGGTTCGTAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:1569694/1‑50 (MQ=255)                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             <  1:464729/50‑1 (MQ=255)                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:2979223/1‑50 (MQ=255)                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:5446825/1‑50 (MQ=255)                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:2108469/1‑50 (MQ=255)                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:6858822/1‑50 (MQ=255)                  tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            >  1:3552945/1‑50 (MQ=255)                  tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            >  1:325593/1‑50 (MQ=255)                   ggTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            >  1:981279/1‑49 (MQ=255)                     ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         >  1:6068890/1‑50 (MQ=255)                      tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:5622661/1‑50 (MQ=255)                      tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:6869079/1‑50 (MQ=255)                      tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:99371/1‑50 (MQ=255)                      tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:1719472/1‑50 (MQ=255)                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:262591/49‑1 (MQ=255)                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:283349/49‑1 (MQ=255)                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:3976571/1‑50 (MQ=255)                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:512328/1‑50 (MQ=255)                        gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:5557173/1‑50 (MQ=255)                        gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:1859531/1‑50 (MQ=255)                          aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  1:441952/1‑50 (MQ=255)                           gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   <  1:3332297/50‑1 (MQ=255)                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:5627684/50‑1 (MQ=255)                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:4285952/50‑1 (MQ=255)                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:2655263/50‑1 (MQ=255)                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  1:6115560/50‑1 (MQ=255)                             cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:6261438/50‑1 (MQ=255)                             cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  >  1:2192906/1‑49 (MQ=255)                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:1682659/1‑50 (MQ=255)                               cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac               <  1:57699/50‑1 (MQ=255)                               cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac               <  1:2514912/50‑1 (MQ=255)                               cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac               <  1:4164766/50‑1 (MQ=255)                               cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:2496221/1‑49 (MQ=255)                                ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              >  1:4999724/1‑50 (MQ=255)                                ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              <  1:1671557/50‑1 (MQ=255)                                   cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            >  1:2756813/1‑49 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:6219778/50‑1 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:2844389/50‑1 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:882250/50‑1 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:5450825/50‑1 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:3209182/50‑1 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:1052378/1‑49 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:5283398/1‑49 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:3676037/1‑49 (MQ=255)                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:5157372/1‑49 (MQ=255)                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:2669779/1‑50 (MQ=255)                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  1:3763546/1‑50 (MQ=255)                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         <  1:247510/50‑1 (MQ=255)                                      cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        <  1:2193185/50‑1 (MQ=255)                                      cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        <  1:4910180/50‑1 (MQ=255)                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  1:1950040/1‑50 (MQ=255)                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  1:2845886/50‑1 (MQ=255)                                        aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:259666/1‑50 (MQ=255)                                        aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:3400798/1‑50 (MQ=255)                                         tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc     <  1:6127494/50‑1 (MQ=255)                                         tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:5740404/1‑49 (MQ=255)                                         tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:3384707/1‑49 (MQ=255)                                         tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:1621933/1‑49 (MQ=255)                                          gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    >  1:2021741/1‑50 (MQ=255)                                          gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    >  1:6716259/1‑50 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    <  1:6113390/49‑1 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:4005993/50‑1 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:5587323/1‑50 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:1648629/1‑50 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:817237/1‑50 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:1535742/50‑1 (MQ=255)                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:2964397/1‑50 (MQ=255)                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:5747257/50‑1 (MQ=255)                                               |                                               ATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485860‑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: ‑