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

ACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTT  >  NC_000913/4485874‑4485937                                 |                                aCGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                               <  1:5196652/35‑1 (MQ=255)  cGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                              <  1:987210/35‑1 (MQ=255)   gTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                             <  1:5865814/35‑1 (MQ=255)   gTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                             >  1:2095813/1‑35 (MQ=255)   gTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                             >  1:2095734/1‑35 (MQ=255)   gTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                             >  1:7282586/1‑35 (MQ=255)    tGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                            >  1:4160418/1‑35 (MQ=255)    tGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                            >  1:6946281/1‑35 (MQ=255)      gTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                          >  1:3831566/1‑35 (MQ=255)      gTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                          >  1:6458267/1‑35 (MQ=39)         cGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                       >  1:7716470/1‑35 (MQ=255)         cGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                       >  1:1246081/1‑35 (MQ=255)          gAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                      <  1:4649681/35‑1 (MQ=39)          gAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                      <  1:7614586/35‑1 (MQ=255)           aaGACGCCTACGACGATGCAGGAACTCCGCTGttt                     >  1:5424265/1‑35 (MQ=37)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     >  1:5424132/1‑35 (MQ=255)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     >  1:4115199/1‑35 (MQ=255)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     <  1:4539287/35‑1 (MQ=255)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     <  1:4589362/35‑1 (MQ=39)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     >  1:952638/1‑35 (MQ=255)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     >  1:814371/1‑35 (MQ=255)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     >  1:804587/1‑35 (MQ=39)           aaGACGCCCACGACGATGCAGGAACTCCGCTGttt                     >  1:2143935/1‑35 (MQ=255)            aGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                    >  1:2631999/1‑35 (MQ=255)              aCGCCCACGACGATGCAGGAACTCCGCTGTTtctc                  >  1:5068006/1‑35 (MQ=255)              aCGCCCACGACGATGCAGGAACTCCGCTGTTtctc                  <  1:6714441/35‑1 (MQ=255)              aCGCCCACGACGATGCAGGAACTCCGCTGTTtctc                  <  1:1059316/35‑1 (MQ=39)              aCGCCCACGACGATGCAGGAACTCCGCTGTTtctc                  >  1:2707826/1‑35 (MQ=39)               cGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                 <  1:150463/35‑1 (MQ=39)                gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                <  1:3854774/35‑1 (MQ=39)                ccccACGACGATGCAGGAACTCCGCTGTTTCTCCg                <  1:365039/34‑1 (MQ=39)                 cccACGACGATGCAGGAACTCCGCTGTTTCTCCgg               >  1:3790545/1‑35 (MQ=255)                 cccACGACGATGCAGGAACTCCGCTGTTTCTCCgg               <  1:1315285/35‑1 (MQ=39)                 cccACGACGATGCAGGAACTCCGCTGTTTCTCCgg               >  1:5041195/1‑35 (MQ=255)                  ccACGACGATGCAGGAACTCCGCTGTTTCTCCggg              <  1:6383432/35‑1 (MQ=39)                   cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc             <  1:5915812/35‑1 (MQ=39)                   cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc             <  1:4667871/35‑1 (MQ=255)                    acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCt            >  1:5176680/1‑35 (MQ=255)                    acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCt            >  1:1600311/1‑35 (MQ=255)                      gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTAc          <  1:5789474/35‑1 (MQ=39)                       acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTAcc         <  1:3334335/35‑1 (MQ=255)                       acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTAcc         >  1:7539414/1‑35 (MQ=255)                       acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTAcc         >  1:2966205/1‑35 (MQ=255)                       acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTAcc         <  1:4486472/35‑1 (MQ=255)                        cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCg        >  1:3166498/1‑35 (MQ=255)                          aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt      >  1:425273/1‑35 (MQ=255)                          aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt      <  1:1632749/35‑1 (MQ=38)                          aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt      <  1:2411443/35‑1 (MQ=255)                          aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt      >  1:1564633/1‑35 (MQ=255)                          aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt      <  1:3412402/35‑1 (MQ=38)                          aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt      >  1:2546529/1‑35 (MQ=38)                           tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt     <  1:3728607/35‑1 (MQ=255)                           tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt     <  1:7482397/35‑1 (MQ=255)                           tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt     <  1:976312/35‑1 (MQ=255)                            gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt    <  1:6037893/35‑1 (MQ=255)                            gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt    <  1:5828001/35‑1 (MQ=38)                            gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt    <  1:1805060/35‑1 (MQ=255)                            gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt    >  1:4055655/1‑35 (MQ=38)                            gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt    <  1:4452585/35‑1 (MQ=255)                            gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt    <  1:160904/35‑1 (MQ=38)                             cAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt   >  1:1023569/1‑35 (MQ=255)                             cAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt   <  1:2507141/35‑1 (MQ=38)                             cAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt   <  1:582089/35‑1 (MQ=255)                             cAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt   >  1:3016159/1‑35 (MQ=38)                             cAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt   >  1:4230749/1‑35 (MQ=38)                              aGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt  <  1:5104204/35‑1 (MQ=255)                              aGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt  >  1:3859761/1‑35 (MQ=255)                                 |                                ACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTT  >  NC_000913/4485874‑4485937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑