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	98.5%	248.9 / ‑6.3	67	P163S (CCA→TCA)	asnS	asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (31/35); minor base T (0/1); total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAG  >  NC_000913/988454‑988544

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

ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGtt                                           >  1:6036422/1‑50 (MQ=255)

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

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

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

  aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTg                                          >  1:2003591/1‑49 (MQ=255)

  aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                         <  1:230539/50‑1 (MQ=255)

  aaCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                         <  1:3755748/50‑1 (MQ=255)

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

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

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

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

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

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

     aTTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAc                                      <  1:3810654/50‑1 (MQ=255)

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

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

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

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

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

          gCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                  >  1:3477462/1‑49 (MQ=255)

          gCCTGCACCTTCGGTATAAGATGCGGTAATCAGTGACGTTGAAGCCCag                                  <  1:3797321/49‑1 (MQ=255)

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

           ccTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCag                                  >  1:4640155/1‑48 (MQ=255)

            cTGCACCTTCGGTATCAGATGCGGTAATCAGTGTCGTTGAAACCCagaag                               <  1:1935460/50‑1 (MQ=255)

             tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                               >  1:5342113/1‑49 (MQ=255)

             tGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaag                               >  1:1672862/1‑49 (MQ=255)

              gCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAgaagaa                             >  1:3552835/1‑50 (MQ=255)

               cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCagaaga                              <  1:5514013/48‑1 (MQ=255)

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

               cACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAAt                            <  1:5368296/50‑1 (MQ=255)

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

                aCCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATc                           <  1:6275086/50‑1 (MQ=255)

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

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

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

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

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

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

                   ttCGGTATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCt                        >  1:6419594/1‑50 (MQ=255)

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

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

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

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

                        tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                   >  1:1496720/1‑50 (MQ=255)

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

                          tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                 <  1:3355867/50‑1 (MQ=255)

                          tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                 <  1:2781668/50‑1 (MQ=255)

                          tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGtt                 <  1:6471755/50‑1 (MQ=255)

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

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

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

                               tGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaag             >  1:547144/1‑49 (MQ=255)

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

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

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

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

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

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

                                      aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg     <  1:2635056/50‑1 (MQ=255)

                                      aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg     <  1:2652096/50‑1 (MQ=255)

                                      aTCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATg     <  1:3613995/50‑1 (MQ=255)

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

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

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

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

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

GGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAG  >  NC_000913/988454‑988544

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