BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,513,003	A→G	D212D (GAT→GAC)	trpS ←	tryptophan‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,513,003	0	A	G	99.5%	752.2 / NA	190	D212D (GAT→GAC)	trpS	tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (101/88); minor base T (0/1); total (101/89)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 190 total aligned reads displayed.

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             >  2:3137970/1‑50 (MQ=255)

 cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  2:3403681/1‑50 (MQ=255)

 cgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                            >  1:2958922/1‑50 (MQ=255)

  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  2:3689145/1‑50 (MQ=255)

  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  2:1901785/1‑50 (MQ=255)

  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           >  1:2338902/1‑50 (MQ=255)

  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCtt                                           <  2:1122813/50‑1 (MQ=255)

    aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:148668/1‑49 (MQ=255)

    aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:58056/1‑49 (MQ=255)

     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  2:3020273/50‑1 (MQ=255)

     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  2:5254350/50‑1 (MQ=255)

     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  2:5598334/50‑1 (MQ=255)

     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGACTTc                                          >  1:4732520/1‑48 (MQ=38)

      tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  2:4664477/1‑50 (MQ=255)

      tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:4666852/1‑50 (MQ=255)

      tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  2:3746952/1‑50 (MQ=255)

      tCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       <  2:317713/50‑1 (MQ=255)

       cAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                       >  1:3649126/1‑49 (MQ=255)

        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGTTCTTCcagc                                      <  2:1080440/49‑1 (MQ=255)

         gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    >  1:898257/1‑50 (MQ=255)

         gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    <  1:5686586/50‑1 (MQ=255)

         gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagcag                                    >  2:1269638/1‑50 (MQ=255)

          tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   >  1:789358/1‑50 (MQ=255)

          tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   >  1:218792/1‑50 (MQ=255)

          tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                   >  2:5261085/1‑50 (MQ=255)

           gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  2:1475431/1‑50 (MQ=255)

            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  2:3657245/1‑50 (MQ=255)

            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  2:392769/50‑1 (MQ=255)

            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  <  1:5175327/49‑1 (MQ=255)

             ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 >  1:4164462/1‑49 (MQ=255)

              cGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:4822356/1‑50 (MQ=255)

               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  2:3085669/50‑1 (MQ=255)

               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:5358969/50‑1 (MQ=255)

                cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  2:2908608/1‑49 (MQ=255)

                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  1:4291750/50‑1 (MQ=255)

                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  2:5709603/50‑1 (MQ=255)

                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  2:182563/50‑1 (MQ=255)

                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           <  1:3888426/50‑1 (MQ=255)

                   gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:4434441/1‑50 (MQ=255)

                    tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  2:672996/50‑1 (MQ=255)

                    tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  2:2204895/1‑50 (MQ=255)

                    tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         <  2:1462907/50‑1 (MQ=255)

                    tttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:973804/1‑50 (MQ=255)

                     ttGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGt                         >  1:5731829/1‑49 (MQ=255)

                      tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  2:2398169/49‑1 (MQ=255)

                      tGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGtt                        <  1:4036511/49‑1 (MQ=255)

                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTa                       <  2:2648244/49‑1 (MQ=255)

                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:3241492/1‑50 (MQ=255)

                       gATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      >  1:5534436/1‑50 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:2624908/49‑1 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  2:2000719/49‑1 (MQ=255)

                        aTTTTCTTCACTACCGATTACGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:5303991/50‑1 (MQ=255)

                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  2:5187110/50‑1 (MQ=255)

                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    <  2:2823884/50‑1 (MQ=255)

                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:1854728/1‑50 (MQ=255)

                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  2:1538765/1‑50 (MQ=255)

                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:2472432/1‑50 (MQ=255)

                         ttttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTg                    >  1:5615095/1‑50 (MQ=255)

                          tttCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:298514/1‑50 (MQ=255)

                           ttcttcTCTACCGAATTCGGGTCTTCCAGCAGAACGATAACGATAttac                   <  2:1926452/49‑3 (MQ=255)

                           ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  2:1697531/49‑1 (MQ=255)

                           ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   >  1:396295/1‑49 (MQ=255)

                           ttcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGc                   <  1:1209250/49‑1 (MQ=255)

                            tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  2:4485281/49‑1 (MQ=255)

                            tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  2:618223/49‑1 (MQ=255)

                            tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  2:3749083/49‑1 (MQ=255)

                            tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:4372826/49‑1 (MQ=255)

                            tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGa                 >  2:5364822/1‑50 (MQ=255)

                             cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:248239/1‑50 (MQ=255)

                              ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:3354129/1‑49 (MQ=255)

                               tcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              >  1:1163476/1‑50 (MQ=255)

                               tcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              >  1:3442877/1‑50 (MQ=255)

                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  2:4623556/48‑1 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  2:2723543/1‑50 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  2:2283675/1‑50 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  2:5384133/1‑50 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:4348346/50‑1 (MQ=255)

                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  2:2921841/50‑1 (MQ=255)

                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:1093830/1‑50 (MQ=255)

                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:2060687/50‑1 (MQ=255)

                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:286160/1‑50 (MQ=255)

                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           <  1:3849629/50‑1 (MQ=255)

                                   tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:157051/1‑50 (MQ=255)

                                   tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  2:4899181/1‑50 (MQ=255)

                                    aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:5438193/1‑49 (MQ=255)

                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:5240532/49‑1 (MQ=255)

                                      cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  2:777978/1‑47 (MQ=255)

                                      cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:4091140/1‑47 (MQ=255)

                                      cGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAg       <  2:2084200/50‑1 (MQ=255)

                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  2:3005293/50‑1 (MQ=255)

                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  2:4054834/1‑50 (MQ=255)

                                       aaTTTCGGGTCTTCCAGCAAGCCGATAACGTTATTGCGATTATCGTCAGa      >  2:4186987/2‑50 (MQ=38)

                                        aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     <  1:3632916/50‑1 (MQ=255)

                                        aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     <  1:307730/50‑1 (MQ=255)

                                        aTTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGAc     >  1:2459103/1‑50 (MQ=255)

                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  1:5117994/1‑50 (MQ=255)

                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    >  1:1739254/1‑50 (MQ=255)

                                         tttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACt    <  2:594592/50‑1 (MQ=255)

                                           tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  2:3266230/1‑50 (MQ=255)

                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  2:2133686/1‑49 (MQ=255)

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

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