BRESEQ :: Evidence

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

ATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485860‑4485951

aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                             >  1:621566/1‑49 (MQ=255)

   ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                          >  1:598304/1‑49 (MQ=255)

   ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                          >  1:1142623/1‑49 (MQ=255)

   ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                           <  1:2097993/48‑1 (MQ=255)

    gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                         >  1:2076192/1‑49 (MQ=255)

       aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                      >  1:1033079/1‑49 (MQ=255)

       aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                      <  1:1880218/49‑1 (MQ=255)

        aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                     >  1:1061769/1‑49 (MQ=255)

          ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                   >  1:1792293/1‑49 (MQ=255)

            cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                 >  1:998730/1‑49 (MQ=255)

            cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                 <  1:685303/49‑1 (MQ=255)

             gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                <  1:1295032/49‑1 (MQ=255)

              aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                               >  1:1588741/1‑49 (MQ=255)

              aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                               <  1:1144781/49‑1 (MQ=255)

                gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                             >  1:1017022/1‑49 (MQ=255)

                  ggTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                           >  1:441189/1‑49 (MQ=255)

                   gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                          <  1:1677726/49‑1 (MQ=255)

                   gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCt                          <  1:1227025/49‑1 (MQ=255)

                       gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                       <  1:1280685/48‑1 (MQ=255)

                        aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                     >  1:1742213/1‑49 (MQ=255)

                             gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                 >  1:1337100/1‑48 (MQ=255)

                                  cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac            >  1:1004338/1‑48 (MQ=255)

                                   gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt           >  1:756734/1‑48 (MQ=255)

                                     cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa        <  1:205326/49‑1 (MQ=255)

                                     cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa        <  1:1104600/49‑1 (MQ=255)

                                      gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc       >  1:1183239/1‑49 (MQ=255)

                                       aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt      <  1:638529/49‑1 (MQ=255)

                                       aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt      >  1:1223230/1‑49 (MQ=255)

                                        tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc     <  1:48770/49‑1 (MQ=255)

                                         gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc     >  1:1255051/1‑48 (MQ=255)

                                          cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc    <  1:710863/48‑1 (MQ=255)

                                          cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa   <  1:1894983/49‑1 (MQ=255)

                                          cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa   <  1:1063821/49‑1 (MQ=255)

                                          cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa   <  1:776592/49‑1 (MQ=255)

                                           aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt  <  1:1963735/49‑1 (MQ=255)

                                           aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt  >  1:215254/1‑49 (MQ=255)

ATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485860‑4485951

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