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	87.0%	362.7 / 45.3	131	D212D (GAT→GAC)	trpS	tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base A (9/8); new base G (66/48); total (75/56)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

Only 100 of 131 total aligned reads displayed.

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTGGA  >  NC_000913/3512956‑3513050

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

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

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

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

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

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

  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                              <  2:1063881/49‑1 (MQ=255)

  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                              <  2:806070/49‑1 (MQ=255)

  ggAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCt                                              <  2:4102579/49‑1 (MQ=255)

    aGTCAGTTACCGCACGTTTTATTTTCTTCACTACCGATTTCGGGTCTTcc                                           <  2:4491215/50‑1 (MQ=255)

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

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

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

    aGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTc                                            >  1:4433119/1‑49 (MQ=255)

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

        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagc                                        <  1:1994189/49‑1 (MQ=255)

        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                         >  1:1548235/1‑48 (MQ=255)

        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcag                                         >  1:563693/1‑48 (MQ=255)

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

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

         gTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCcagcag                                      >  2:2904531/1‑50 (MQ=255)

          tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                     <  2:3935705/50‑1 (MQ=255)

          tGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAgg                                     <  1:579860/50‑1 (MQ=255)

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

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

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

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

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

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

             ccGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGcc                                   >  1:914172/1‑49 (MQ=255)

             ccGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGcc                                   >  2:1777573/1‑49 (MQ=255)

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

               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                                 >  1:3713976/1‑49 (MQ=255)

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

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

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

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

                 aCGTTTGATTTTCTTCAGTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                              <  1:3986848/50‑1 (MQ=38)

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

                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                              >  1:4555343/1‑50 (MQ=255)

                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                              >  1:3940372/1‑50 (MQ=255)

                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                             >  1:4334462/1‑50 (MQ=255)

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

                   gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                            <  2:340445/50‑1 (MQ=255)

                   gTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACg                            <  1:3030009/50‑1 (MQ=255)

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

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

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

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

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                       >  1:1850218/1‑50 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                       >  1:742743/1‑50 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                       <  2:1410883/50‑1 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTAtt                       >  2:4164147/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

                            tcttcACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCg                    <  2:2985541/49‑1 (MQ=255)

                            tcttcACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCg                    <  1:1814753/49‑1 (MQ=255)

                             cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGa                   <  2:825234/49‑1 (MQ=255)

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

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

                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt               >  1:3701833/1‑50 (MQ=255)

                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt               >  2:2316111/1‑50 (MQ=255)

                                cACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGAtt                 <  1:3325039/48‑1 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTAt               <  1:1325774/49‑1 (MQ=255)

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

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

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

                                 aCTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATc              <  1:2463088/50‑1 (MQ=255)

                                  cTACCTATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg             <  2:351410/50‑1 (MQ=255)

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

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

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

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

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

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

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

                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt            <  1:3419143/48‑1 (MQ=255)

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

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

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

                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGCTTATCGTCAGa        >  2:4792762/1‑50 (MQ=255)

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

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

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

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

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

                                          ttCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACtt     >  2:2494449/1‑50 (MQ=255)

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

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

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

                                           tCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg    >  2:4913575/1‑50 (MQ=255)

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

                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg    <  2:4469843/49‑1 (MQ=255)

                                             ggATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTGGa  <  2:3947031/50‑1 (MQ=255)

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTGGA  >  NC_000913/3512956‑3513050

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