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	92.2%	434.1 / 23.0	129	P163S (CCA→TCA)	asnS	asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base G (7/3); new base A (60/59); total (67/62)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 131 total aligned reads displayed.

ACGCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGCGC  >  NC_000913/988450‑988547

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

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

   cgGAACATTTTGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                               <  2:2198370/50‑1 (MQ=255)

   cgGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGt                                               >  1:4885319/1‑50 (MQ=255)

    ggAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGtt                                              >  2:330575/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

       aCATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGa                                            >  1:3372347/1‑49 (MQ=255)

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

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

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

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

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

          tttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                        >  1:2449603/1‑50 (MQ=255)

          tttCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGACGTTGAAAcc                                        >  1:2225679/1‑50 (MQ=255)

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

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

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

             cGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAAcc                                        >  2:2135861/1‑47 (MQ=255)

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

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

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

                cTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCag                                     >  2:2328669/1‑47 (MQ=255)

                cTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCag                                     >  2:1470204/1‑47 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                        tCGGGATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTg                          >  2:2388906/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                            tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                      <  1:2435105/50‑1 (MQ=255)

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

                            tATCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCg                      <  1:3567904/50‑1 (MQ=255)

                              tCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGtt                    <  2:3564550/50‑1 (MQ=255)

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

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

                              tCAGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGt                     >  2:1669138/1‑49 (MQ=255)

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

                                aGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTa                   >  1:1588726/1‑49 (MQ=255)

                                aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaa                  <  2:4489244/50‑1 (MQ=255)

                                aGATGCGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTaa                  <  1:4736317/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                     cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa              <  1:1716919/49‑1 (MQ=255)

                                     cGGTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAaagaa              <  2:1917507/49‑1 (MQ=255)

                                      ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc            >  2:5293365/1‑50 (MQ=255)

                                      ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc            >  1:2061770/1‑50 (MQ=255)

                                      ggTAATCAGTGACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGc            <  1:3887818/50‑1 (MQ=255)

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

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

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

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

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

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

                                          aTCAGTAACGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGAt         <  1:5302039/49‑1 (MQ=255)

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

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

                                             aGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCa      <  2:2262602/49‑1 (MQ=255)

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

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

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

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

                                                ggCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAgcgc  >  1:3815784/1‑50 (MQ=255)

ACGCGGAACATTTCGCCTGCACCTTCGGTATCAGATGCGGTAATCAGTGGCGTTGAAACCCAGAAGAATCCCTGCTCGTTAAAGAAGCGATGCAGCGC  >  NC_000913/988450‑988547

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