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_0009134,485,9060CA100.0% 102.3 / NA 26A16S (GCC→TCC) pepAaminopeptidase A/I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/14);  total (12/14)

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952
                                               |                                              
aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:568241/1‑50 (MQ=255)
  tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            >  1:1714549/1‑50 (MQ=255)
   cGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           >  1:5595/1‑50 (MQ=255)
       gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       >  1:1196242/1‑50 (MQ=255)
         aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                     <  1:2005734/50‑1 (MQ=255)
           ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   >  1:1981184/1‑50 (MQ=255)
            gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  <  1:1810085/50‑1 (MQ=255)
              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:769978/1‑50 (MQ=255)
              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:714632/50‑1 (MQ=255)
               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:1604017/50‑1 (MQ=255)
                cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:1573312/50‑1 (MQ=255)
                    gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                           <  1:657607/49‑1 (MQ=255)
                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  1:1635700/49‑1 (MQ=255)
                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:1006933/1‑50 (MQ=255)
                         aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     >  1:1209715/1‑50 (MQ=255)
                           gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   >  1:1214226/1‑50 (MQ=255)
                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:1849363/1‑50 (MQ=255)
                               cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac               <  1:950917/50‑1 (MQ=255)
                                 cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             <  1:2122966/50‑1 (MQ=255)
                                  acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:458885/1‑49 (MQ=255)
                                   cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           <  1:24108/50‑1 (MQ=255)
                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:2240193/50‑1 (MQ=255)
                                      cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        <  1:562510/50‑1 (MQ=255)
                                      cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:1805032/1‑50 (MQ=255)
                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACGCCAt   <  1:1859062/50‑1 (MQ=255)
                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:59618/50‑1 (MQ=255)
                                               |                                              
AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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