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

GGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485863‑4485951

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

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

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

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

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

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

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

      aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                    >  1:305935/1‑49 (MQ=255)

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

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

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

          gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                >  1:2273424/1‑49 (MQ=255)

          gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                >  1:885604/1‑49 (MQ=255)

          gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                >  1:1767214/1‑49 (MQ=255)

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

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

            cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                              >  1:1166040/1‑49 (MQ=255)

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

             gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                             <  1:654931/49‑1 (MQ=255)

                gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                           >  1:1762986/1‑48 (MQ=255)

                 ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTa                         >  1:1135741/1‑49 (MQ=255)

                 ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTa                         >  1:1578233/1‑49 (MQ=255)

                   cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                        <  1:1988468/48‑1 (MQ=255)

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

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

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

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

                    gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                      >  1:2498732/1‑49 (MQ=255)

                    gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                      >  1:1319733/1‑49 (MQ=255)

                       gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                   >  1:1855343/1‑49 (MQ=255)

                        aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                  <  1:2040885/49‑1 (MQ=255)

                         cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                 <  1:875436/49‑1 (MQ=255)

                          gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                <  1:1513884/49‑1 (MQ=255)

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

                           cccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa               >  1:2350805/1‑49 (MQ=255)

                             cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca             >  1:2209101/1‑49 (MQ=255)

                             cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca             <  1:99716/49‑1 (MQ=255)

                              acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac            <  1:1261402/49‑1 (MQ=255)

                              acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac            >  1:2054289/1‑49 (MQ=255)

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

                                 acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa         <  1:2274165/49‑1 (MQ=255)

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

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

                                     tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc     >  1:2446152/1‑49 (MQ=255)

                                     tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc     >  1:1965898/1‑49 (MQ=255)

                                       cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa   >  1:2576304/1‑49 (MQ=255)

                                       cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa   >  1:2200900/1‑49 (MQ=255)

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

GGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAT  >  NC_000913/4485863‑4485951

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