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

Only 100 of 130 total aligned reads displayed.

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:1110927/1‑50 (MQ=255)

aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              >  1:3056729/1‑50 (MQ=255)

aaTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAAc                                              <  1:2653361/50‑1 (MQ=255)

 aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             >  2:2736879/1‑50 (MQ=255)

 aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             <  1:1526700/50‑1 (MQ=255)

 aTCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACt                                             <  1:2949062/50‑1 (MQ=255)

  tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            >  1:105597/1‑50 (MQ=255)

  tCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTc                                            <  2:497891/50‑1 (MQ=255)

   cGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           >  2:1678933/1‑50 (MQ=255)

   cGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           >  2:558782/1‑50 (MQ=255)

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcct                                          >  1:2534586/1‑49 (MQ=255)

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:3578266/49‑1 (MQ=255)

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:4094104/49‑1 (MQ=255)

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:1864856/49‑1 (MQ=255)

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTcc                                           <  1:1757104/49‑1 (MQ=255)

      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        <  2:3566115/50‑1 (MQ=255)

       gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  1:112473/50‑1 (MQ=255)

       gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       <  2:3190056/50‑1 (MQ=255)

        aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      >  2:4261021/1‑50 (MQ=255)

         aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGTTGtc                                     >  2:667765/1‑49 (MQ=255)

         aaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGtt                                     >  1:3831318/1‑50 (MQ=255)

          aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    >  2:3041304/1‑50 (MQ=255)

          aGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGttt                                    <  2:1041553/50‑1 (MQ=255)

            gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  >  2:138763/1‑50 (MQ=255)

              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:3443197/50‑1 (MQ=255)

              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  1:2848994/50‑1 (MQ=255)

              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                <  2:1547805/50‑1 (MQ=255)

                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  1:783094/1‑50 (MQ=255)

                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  2:1068685/1‑50 (MQ=255)

                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             <  1:270133/50‑1 (MQ=255)

                 gTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCggg                             >  2:2932863/1‑50 (MQ=255)

                 gTGGTTCGAAGACGCACACGACGATGTAGGAACTCCGCTGTTTCCCCggg                             >  2:3614911/1‑50 (MQ=255)

                  tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            >  1:1430838/1‑50 (MQ=255)

                     ttCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAc                         <  1:405258/50‑1 (MQ=255)

                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        <  2:2535420/49‑1 (MQ=255)

                       cGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       >  1:4271385/1‑50 (MQ=255)

                        gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGc                      >  1:3832770/1‑50 (MQ=255)

                         aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     >  1:4354168/1‑50 (MQ=255)

                         aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     <  1:947807/50‑1 (MQ=255)

                         aaGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCt                     >  2:525945/1‑50 (MQ=255)

                          aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    >  2:143567/1‑50 (MQ=255)

                           gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   <  2:2160594/50‑1 (MQ=255)

                           gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   <  1:1361969/50‑1 (MQ=255)

                           gACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttt                   <  1:2908430/50‑1 (MQ=255)

                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  2:3462344/50‑1 (MQ=255)

                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  <  2:3977567/50‑1 (MQ=255)

                             cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  1:2754184/50‑1 (MQ=255)

                             cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCttttt                 <  2:4430811/50‑1 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                <  2:1414096/50‑1 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:390968/1‑50 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:2133653/1‑50 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                <  2:2656127/50‑1 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                <  1:4152087/50‑1 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  2:3413715/1‑50 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                <  1:1861990/50‑1 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:3961058/1‑50 (MQ=255)

                              gCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTa                >  1:2821224/1‑50 (MQ=255)

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

                                ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              <  1:2620557/50‑1 (MQ=255)

                                ccACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTaca              >  1:2341467/1‑50 (MQ=255)

                                 cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:2715606/1‑50 (MQ=255)

                                 cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             >  1:2427292/1‑50 (MQ=255)

                                 cACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac             <  1:4546185/50‑1 (MQ=255)

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

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

                                   cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            >  1:4400475/1‑49 (MQ=255)

                                   cgacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           <  2:2097593/50‑1 (MQ=255)

                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  2:729496/1‑49 (MQ=255)

                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  2:3304607/1‑49 (MQ=255)

                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa           >  1:320762/1‑49 (MQ=255)

                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         >  2:4198403/1‑50 (MQ=255)

                                      cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        >  1:3858778/1‑50 (MQ=255)

                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  1:2077886/50‑1 (MQ=255)

                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  2:2129504/1‑50 (MQ=255)

                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  2:2987362/50‑1 (MQ=255)

                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  2:1469436/1‑50 (MQ=255)

                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       >  2:3139832/1‑50 (MQ=255)

                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  1:4107958/50‑1 (MQ=255)

                                       gaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt       <  1:602057/50‑1 (MQ=255)

                                        aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACta      >  1:796295/1‑49 (MQ=255)

                                        aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      <  2:2513656/50‑1 (MQ=255)

                                         tGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc     <  2:3929249/50‑1 (MQ=255)

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

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

                                          gCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCa    >  2:708805/1‑50 (MQ=255)

                                           cAGGAACTCCGTTGTTCCCCCGGGCTCCCGTTTTTTACACAAAACCCCAt   <  2:1957524/50‑1 (MQ=255)

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

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

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

                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  2:3388714/50‑1 (MQ=255)

                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  2:3497014/1‑50 (MQ=255)

                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  2:1920397/50‑1 (MQ=255)

                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   >  1:3363892/1‑50 (MQ=255)

                                           cAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt   <  1:2322112/50‑1 (MQ=255)

                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:1008707/1‑50 (MQ=255)

                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:1972345/50‑1 (MQ=255)

                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  2:2239852/1‑50 (MQ=255)

                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  1:3429655/1‑50 (MQ=255)

                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  <  1:529572/50‑1 (MQ=255)

                                            aGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg  >  2:3395336/1‑50 (MQ=255)

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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