BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	988,499	G→A	P163S (CCA→TCA)	asnS ←	asparagine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	988,499	0	G	A	100.0%	536.7 / NA	133	P163S (CCA→TCA)	asnS	asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (61/72); total (61/72)

Only 100 of 133 total aligned reads displayed.

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545

gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              >  1:4016056/1‑50 (MQ=255)

gcgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACg                                              <  2:1952561/50‑1 (MQ=255)

  ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:4754691/50‑1 (MQ=255)

  ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  1:5039290/50‑1 (MQ=255)

  ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            <  2:1146782/50‑1 (MQ=255)

   gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                            >  2:3023313/1‑49 (MQ=255)

   gAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                           <  2:3704032/50‑1 (MQ=255)

     aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaa                                         <  2:4265069/50‑1 (MQ=255)

     aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaa                                         <  2:3498732/50‑1 (MQ=255)

      cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        >  1:733968/1‑50 (MQ=255)

      cATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGaaa                                        <  1:1917980/50‑1 (MQ=255)

       aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                       >  2:120061/1‑50 (MQ=255)

        tttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      <  1:4300945/50‑1 (MQ=255)

         ttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                      >  1:885183/1‑49 (MQ=255)

          tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  1:1771006/49‑1 (MQ=255)

          tCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAccc                                     <  2:3991450/49‑1 (MQ=255)

           cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCa                                    <  2:4161797/49‑1 (MQ=255)

            gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:4088338/50‑1 (MQ=255)

            gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCaga                                  <  1:4961264/50‑1 (MQ=255)

            gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  2:2063814/49‑1 (MQ=255)

            gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   <  1:2521933/49‑1 (MQ=255)

             ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaa                                 <  2:2343576/50‑1 (MQ=255)

              cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                <  1:5628400/50‑1 (MQ=255)

              cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                                <  2:2500913/50‑1 (MQ=255)

              cTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                   >  1:1111602/1‑47 (MQ=255)

                 cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                              <  2:1894657/49‑1 (MQ=255)

                 cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                             >  1:4386173/1‑50 (MQ=255)

                  aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            >  1:4857955/1‑50 (MQ=255)

                  aCCTACGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                            >  1:1352087/1‑50 (MQ=255)

                   ccTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATcc                           <  2:846965/50‑1 (MQ=255)

                    cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          <  1:2143547/50‑1 (MQ=255)

                    cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  2:4212637/1‑50 (MQ=255)

                    cTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATccc                          >  1:5340455/1‑50 (MQ=255)

                      tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  2:168087/1‑50 (MQ=255)

                      tCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        >  1:3505367/1‑50 (MQ=255)

                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  1:1258146/48‑1 (MQ=255)

                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  2:5224343/48‑1 (MQ=255)

                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  2:3634031/48‑1 (MQ=255)

                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                         <  2:4455129/48‑1 (MQ=255)

                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  1:293448/49‑1 (MQ=255)

                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                        <  2:1588382/49‑1 (MQ=255)

                       cGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGc                       >  1:3923242/1‑50 (MQ=255)

                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:2198071/1‑50 (MQ=255)

                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:856164/1‑50 (MQ=255)

                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:5651239/1‑50 (MQ=255)

                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:5251383/1‑50 (MQ=255)

                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  2:4564736/50‑1 (MQ=255)

                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      >  1:4538584/1‑50 (MQ=255)

                        ggTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCt                      <  2:1766386/50‑1 (MQ=255)

                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:1240568/50‑1 (MQ=255)

                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     >  1:351676/1‑50 (MQ=255)

                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     >  1:336696/1‑50 (MQ=255)

                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:320240/50‑1 (MQ=255)

                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:4481131/50‑1 (MQ=255)

                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  1:4215563/50‑1 (MQ=255)

                         gTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTc                     <  2:1136790/50‑1 (MQ=255)

                           aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  2:2354795/1‑50 (MQ=255)

                           aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  2:467301/1‑50 (MQ=255)

                           aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  1:4610852/1‑50 (MQ=255)

                           aTCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   >  2:5625408/1‑50 (MQ=255)

                            tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                   <  2:1236321/49‑1 (MQ=255)

                             cAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                  >  1:2825983/1‑49 (MQ=255)

                              aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:5441349/1‑49 (MQ=255)

                              aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:2917773/1‑49 (MQ=255)

                              aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 <  2:3090433/49‑1 (MQ=255)

                              aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 <  1:4113388/49‑1 (MQ=255)

                               gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaaa               <  2:3264414/50‑1 (MQ=255)

                               gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaaa               <  2:5685876/50‑1 (MQ=255)

                               gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:1009416/1‑48 (MQ=255)

                               gATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:4831126/1‑48 (MQ=255)

                                aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTa                 >  1:2526723/1‑47 (MQ=255)

                                aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  2:38011/50‑1 (MQ=255)

                                aTGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              <  1:3073197/50‑1 (MQ=255)

                                 tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaaga             <  2:2437029/50‑1 (MQ=255)

                                 tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaaga             <  1:1935632/50‑1 (MQ=255)

                                  gCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag              >  1:1409926/1‑48 (MQ=255)

                                   cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           <  2:605514/50‑1 (MQ=255)

                                   cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaag           >  1:4912026/1‑50 (MQ=255)

                                     gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc          >  2:697124/1‑49 (MQ=255)

                                     gTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCg         <  2:4888232/50‑1 (MQ=255)

                                      tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  2:3226133/1‑50 (MQ=255)

                                      tAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGa        >  2:3280389/1‑50 (MQ=255)

                                       aaTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt       >  1:493360/1‑50 (MQ=255)

                                        aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:4978709/1‑50 (MQ=255)

                                         tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  2:4473234/1‑49 (MQ=255)

                                         tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGc     <  2:1189117/50‑1 (MQ=255)

                                         tCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGc     <  1:4522764/50‑1 (MQ=255)

                                          cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  2:4836611/1‑48 (MQ=255)

                                          cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:4380214/1‑48 (MQ=255)

                                          cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  2:2388281/1‑48 (MQ=255)

                                          cAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg      >  1:2554594/1‑48 (MQ=255)

                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:610615/49‑1 (MQ=255)

                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:4016014/49‑1 (MQ=255)

                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  2:1940436/49‑1 (MQ=255)

                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa    <  1:246250/49‑1 (MQ=255)

                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   >  2:4002545/1‑50 (MQ=255)

                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  1:1877275/50‑1 (MQ=255)

                                           aGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAg   <  2:1173861/50‑1 (MQ=255)

                                            gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  >  1:3639862/1‑50 (MQ=255)

                                            gTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgc  <  2:3091119/50‑1 (MQ=255)

GCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGC  >  NC_000913/988452‑988545

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