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

Only 100 of 124 total aligned reads displayed.

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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

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

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

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

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

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

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

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

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

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

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:5056667/1‑50 (MQ=255)

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          >  1:792922/1‑50 (MQ=255)

    ggAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCg                                          <  1:4286704/50‑1 (MQ=255)

     gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         >  1:7260327/1‑50 (MQ=255)

     gagaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGc                                         <  1:3232707/50‑1 (MQ=255)

      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:1495973/1‑50 (MQ=255)

      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:5264899/1‑50 (MQ=255)

      agaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCt                                        >  1:8519162/1‑50 (MQ=255)

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

       gaAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTg                                       >  1:6128117/1‑50 (MQ=255)

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

           ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   <  1:7804361/50‑1 (MQ=255)

           ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   <  1:33543/50‑1 (MQ=255)

            gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  <  1:7861144/50‑1 (MQ=255)

            gCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtct                                  <  1:7654597/50‑1 (MQ=255)

             cGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtctc                                 >  1:1392259/1‑50 (MQ=255)

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

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

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

               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:1007255/50‑1 (MQ=255)

               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:6494068/50‑1 (MQ=255)

               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               <  1:6466444/50‑1 (MQ=255)

                cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              >  1:5618233/1‑50 (MQ=255)

                cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              <  1:803708/50‑1 (MQ=255)

                cGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCgg                              >  1:4449008/1‑50 (MQ=255)

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

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

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

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

                  tGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGc                            <  1:8284372/50‑1 (MQ=255)

                    gTTTGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCCCCGGGCTa                          >  1:1382135/1‑50 (MQ=255)

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

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

                    gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTa                          >  1:6428747/1‑50 (MQ=255)

                    gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTa                          >  1:7256995/1‑50 (MQ=255)

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

                        gAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCg                       <  1:3993640/49‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                            aCGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  >  1:1801602/1‑50 (MQ=255)

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

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

                             cGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtttt                  >  1:8096099/1‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                  acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            <  1:925254/50‑1 (MQ=255)

                                  acgacgATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt            <  1:5017397/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         <  1:4481731/50‑1 (MQ=255)

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

                                     acgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa         <  1:7314631/50‑1 (MQ=255)

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

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

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

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

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

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

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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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