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

Only 100 of 118 total aligned reads displayed.

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

        aaaGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGt                                      <  1:3587855/50‑1 (MQ=255)

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

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

           ggCGACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTtc                                   >  1:717942/1‑50 (MQ=255)

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

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

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

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

              gACGTGGTTCGCAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:4706939/1‑50 (MQ=255)

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

              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:6496768/1‑50 (MQ=255)

              gACGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTcc                                >  1:2957890/1‑50 (MQ=255)

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

               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:1599636/1‑50 (MQ=255)

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

               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:7999305/1‑50 (MQ=255)

               aCGTGGTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCg                               >  1:2069290/1‑50 (MQ=255)

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

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

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

                    gTTCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTT‑‑TCCGGGCTAc                         <  1:3301440/49‑1 (MQ=255)

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

                     ttCGAAGACGCCCACGACGATGCAGGAACTCCGCCGTTTCTCCGGGCTAc                         <  1:3308253/50‑1 (MQ=255)

                      tCGAAGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTAcc                        >  1:8314716/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

                          aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    <  1:1244473/50‑1 (MQ=255)

                          aGACGCCCACGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCtt                    <  1:5625661/50‑1 (MQ=255)

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

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

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

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

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

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

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

                               cccAAGACGATGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac               <  1:6575960/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:6338281/50‑1 (MQ=255)

                                    gacgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa          <  1:2447626/50‑1 (MQ=255)

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

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

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

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

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

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

                                      cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        <  1:2623348/50‑1 (MQ=255)

                                      cgaTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc        <  1:3764008/50‑1 (MQ=255)

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

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

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

                                        aTGCAGGAACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc      >  1:5192622/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

AATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATG  >  NC_000913/4485859‑4485952

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