BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,093,223	A→G	V72V (GTT→GTC)	rhaM ←	L‑rhamnose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,223	0	A	G	90.5%	414.7 / 25.5	127	V72V (GTT→GTC)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand): ref base A (4/8); new base G (64/51); total (68/59)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 128 total aligned reads displayed.

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCTC  >  NC_000913/4093177‑4093270

ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                              <  1:4155276/50‑1 (MQ=255)

ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                              <  2:4435074/50‑1 (MQ=255)

 gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                             <  2:3242458/50‑1 (MQ=255)

 gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                             <  1:4269920/50‑1 (MQ=255)

 gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                             >  1:1406474/1‑50 (MQ=255)

  ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                             >  1:3832001/1‑49 (MQ=255)

  ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACAt                                             >  1:2344306/1‑49 (MQ=255)

   tAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                             >  1:5381578/1‑48 (MQ=255)

    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                           <  2:1699867/49‑1 (MQ=255)

    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                           <  1:1073861/49‑1 (MQ=255)

    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                           <  2:3951200/49‑1 (MQ=255)

     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                         <  2:3651703/50‑1 (MQ=255)

     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                         <  2:2689509/50‑1 (MQ=255)

      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                        <  2:5125631/50‑1 (MQ=255)

      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                        <  1:4313494/50‑1 (MQ=255)

      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                        <  1:2511225/50‑1 (MQ=255)

      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTg                                        <  1:604479/50‑1 (MQ=255)

       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                       >  2:2075118/1‑50 (MQ=255)

        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                      >  1:4244162/1‑50 (MQ=255)

        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                      <  2:3126111/50‑1 (MQ=255)

        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                      <  2:3698201/50‑1 (MQ=255)

         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                     <  1:2121302/50‑1 (MQ=255)

         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                     >  1:4581308/1‑50 (MQ=255)

         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                     >  1:721547/1‑50 (MQ=255)

          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                    <  2:1526315/50‑1 (MQ=255)

           aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                   <  2:1986592/50‑1 (MQ=255)

            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                   >  1:3835089/1‑49 (MQ=255)

            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                  <  1:3267453/50‑1 (MQ=255)

             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                 >  1:5361934/1‑50 (MQ=255)

             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                 >  1:5280390/1‑50 (MQ=255)

             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                 >  2:4326647/1‑50 (MQ=255)

              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                                >  1:4967244/1‑50 (MQ=255)

              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                                >  1:1616010/1‑50 (MQ=255)

              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                                >  1:4482098/1‑50 (MQ=255)

               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                               >  1:1705725/1‑50 (MQ=255)

               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                               >  1:3043986/1‑50 (MQ=255)

               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                               <  2:2496503/50‑1 (MQ=255)

                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCGACCGc                             <  2:853889/50‑1 (MQ=255)

                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                              >  1:2279511/1‑49 (MQ=255)

                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGc                             <  1:4143839/50‑1 (MQ=255)

                  cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                            >  1:2578207/1‑50 (MQ=255)

                  cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                            >  2:3876668/1‑50 (MQ=255)

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCACCCGCAt                           <  2:3588230/50‑1 (MQ=255)

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                            <  2:2713241/49‑1 (MQ=255)

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                           <  2:562877/50‑1 (MQ=255)

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                           >  1:3745427/1‑50 (MQ=255)

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                           >  2:3954950/1‑50 (MQ=255)

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                           <  2:2614665/50‑1 (MQ=255)

                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                          >  2:1735676/1‑50 (MQ=255)

                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  2:3393252/1‑50 (MQ=255)

                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  2:1427244/1‑50 (MQ=255)

                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  2:4024782/1‑50 (MQ=255)

                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  2:3665874/1‑49 (MQ=255)

                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  2:2505022/1‑49 (MQ=255)

                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  1:1441208/1‑49 (MQ=255)

                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  1:1223961/1‑49 (MQ=255)

                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                        <  1:2255484/49‑1 (MQ=255)

                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                        <  1:3059160/49‑1 (MQ=255)

                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                         >  1:4978361/1‑48 (MQ=255)

                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                      >  2:1070006/1‑50 (MQ=255)

                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                      <  1:1103098/50‑1 (MQ=255)

                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                     <  2:1413269/49‑1 (MQ=255)

                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                    >  1:4435309/1‑50 (MQ=255)

                           atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                   >  2:5317026/1‑50 (MQ=255)

                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                  <  2:1652140/50‑1 (MQ=255)

                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                  <  2:5043193/50‑1 (MQ=255)

                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                  <  1:156420/50‑1 (MQ=255)

                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                  <  1:1951093/50‑1 (MQ=255)

                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                  <  2:2583968/50‑1 (MQ=255)

                             ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                 >  1:1168095/1‑50 (MQ=255)

                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct                >  1:1989548/1‑50 (MQ=255)

                               ccaccaACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGttct                <  1:3031434/49‑1 (MQ=255)

                                caccaACGTTGGCAGACATCAGTGCTGGCCTCCGCATTCCAGCGttctt               <  2:292053/49‑1 (MQ=255)

                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt               <  2:1288058/49‑1 (MQ=255)

                                caccaACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGttcttc              <  2:4559742/50‑1 (MQ=255)

                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt               >  2:2940207/1‑48 (MQ=255)

                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt               >  2:4992647/1‑48 (MQ=255)

                                 accaACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGttctt               >  1:5346146/1‑48 (MQ=255)

                                  ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt               >  2:215400/1‑47 (MQ=255)

                                  ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg            <  2:5297864/50‑1 (MQ=255)

                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa           <  2:4892660/50‑1 (MQ=255)

                                    aaCGTTGGCAGACCTCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt          >  2:3589995/1‑50 (MQ=39)

                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt          >  1:1867691/1‑50 (MQ=255)

                                    aaCGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt          <  2:4486153/50‑1 (MQ=255)

                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt         >  1:3896291/1‑50 (MQ=255)

                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt         >  1:3967742/1‑50 (MQ=255)

                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt         >  2:2651144/1‑50 (MQ=255)

                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc        >  2:4052247/1‑49 (MQ=255)

                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc        >  2:415784/1‑49 (MQ=255)

                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa       <  1:4250264/50‑1 (MQ=255)

                                         tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt     >  1:3265224/1‑50 (MQ=255)

                                         tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt     >  2:3066366/1‑50 (MQ=255)

                                          ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc    >  2:5295990/1‑50 (MQ=255)

                                          ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc    >  2:2929284/1‑50 (MQ=255)

                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct   >  1:3733188/1‑50 (MQ=255)

                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct   >  1:119824/1‑50 (MQ=255)

                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct   >  1:4719424/1‑50 (MQ=255)

                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc    <  1:1997380/49‑1 (MQ=255)

                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc    <  2:771847/49‑1 (MQ=255)

                                             aaaCATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtctc  >  2:1009742/1‑49 (MQ=255)

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCTC  >  NC_000913/4093177‑4093270

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