BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	313,795	A→G	14.0%	P98P (CCT→CCC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	313,795	0	A	G	14.0%	255.3 / 29.1	107	P98P (CCT→CCC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base A (37/55); new base G (6/9); total (43/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 108 total aligned reads displayed.

ATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAA  >  Exported/313727‑313865

atatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                                      <  1:3929537/71‑1 (MQ=255)

atatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                                      <  1:4099394/71‑1 (MQ=255)

 tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                                      >  1:3301504/1‑70 (MQ=255)

 tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                                      <  1:4106053/70‑1 (MQ=255)

 tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                                     >  1:3319957/1‑71 (MQ=255)

 tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                                     >  1:488604/1‑71 (MQ=255)

 tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                                     >  1:5692641/1‑71 (MQ=255)

 tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                                     >  1:2168727/1‑71 (MQ=255)

     ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                  <  1:2673492/70‑1 (MQ=255)

     ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                                  <  1:2647068/70‑1 (MQ=255)

            gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAg                                                                       <  1:726102/58‑1 (MQ=255)

               tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                                          >  1:2990771/1‑68 (MQ=255)

                 ggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                       <  1:1591286/69‑1 (MQ=255)

                   ttttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATcg                                                                  <  1:4541670/55‑2 (MQ=255)

                   cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                                 <  1:1300814/57‑1 (MQ=255)

                   cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCaa                                                                 <  1:2323984/57‑1 (MQ=255)

                    tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                         <  1:3241212/64‑1 (MQ=255)

                    tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                         <  1:4010275/64‑1 (MQ=255)

                    tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                         <  1:2361173/64‑1 (MQ=255)

                    tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                         <  1:5227791/64‑1 (MQ=255)

                     ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                                                >  1:4953910/1‑56 (MQ=255)

                     ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                                                >  1:3474709/1‑56 (MQ=255)

                     ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg                                                 >  1:1668823/1‑68 (MQ=255)

                     ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggact                                                  <  1:4523712/70‑3 (MQ=255)

                     ttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                                         <  1:121238/63‑1 (MQ=255)

                      tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTgg                                                            <  1:1020961/59‑2 (MQ=255)

                      tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg                                                 <  1:762063/70‑4 (MQ=255)

                      tCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggactg                                                 <  1:1499681/70‑4 (MQ=255)

                      tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGaca                                                >  1:3964632/1‑71 (MQ=255)

                       cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:3554642/62‑1 (MQ=255)

                       cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGacac                                               >  1:1295142/1‑71 (MQ=255)

                        gCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTa                                                          >  1:23704/1‑59 (MQ=255)

                           aGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGc                                           >  1:4880651/1‑71 (MQ=255)

                             aTCACCTTTCGCACTGTCGTATGTCGGCTGATAGGTCGCAGGATCAATTTGATAGAGAGAGACACCTGCt                                          <  1:2381469/70‑1 (MQ=255)

                             aTCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCtt                                         >  1:499912/1‑71 (MQ=255)

                              tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        <  1:3853621/71‑1 (MQ=255)

                              tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        <  1:5018371/71‑1 (MQ=255)

                               cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        >  1:3171989/1‑70 (MQ=255)

                               cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        >  1:2200629/1‑70 (MQ=255)

                               cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        >  1:67833/1‑70 (MQ=255)

                               cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        >  1:3260999/1‑70 (MQ=255)

                                 ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTTATagagaga                                                    <  1:2828308/56‑1 (MQ=255)

                                 ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                                    <  1:4006488/56‑1 (MQ=255)

                                  cTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                                    <  1:304061/55‑1 (MQ=255)

                                  cTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                  >  1:3375539/1‑57 (MQ=255)

                                   tttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacaggga                                                    <  1:367153/54‑1 (MQ=255)

                                   tttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacaggga                                                    <  1:2690401/54‑1 (MQ=255)

                                   tttCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacaggga                                                    <  1:5173001/54‑1 (MQ=255)

                                   tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:3441785/50‑1 (MQ=255)

                                   tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:4040099/50‑1 (MQ=255)

                                    ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGc                                           >  1:5028477/1‑62 (MQ=255)

                                     tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTgg                                                            >  1:4448477/1‑43 (MQ=255)

                                     tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTgg                                                            >  1:3643475/1‑43 (MQ=255)

                                     tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGTATCAATCTGATAGAGAGAGacac                                               >  1:512899/1‑57 (MQ=255)

                                     tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                                              >  1:4390208/1‑42 (MQ=255)

                                     tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                                       >  1:817059/1‑49 (MQ=255)

                                     tCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGacac                                               >  1:2912821/1‑57 (MQ=255)

                                       gCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                                              <  1:2331018/40‑1 (MQ=255)

                                           tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:2967019/42‑1 (MQ=255)

                                           tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:40323/42‑1 (MQ=255)

                                           tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:2973631/42‑1 (MQ=255)

                                           tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                                  <  1:2271121/64‑1 (MQ=255)

                                           tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTc                                  <  1:5569640/64‑1 (MQ=255)

                                             tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagaga                                                      >  1:5534865/1‑42 (MQ=255)

                                             tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCt                               <  1:5526433/65‑1 (MQ=255)

                                             tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttc                         <  1:3383050/71‑1 (MQ=255)

                                             tCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACct                           <  1:3160888/69‑1 (MQ=255)

                                               gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttctt                       >  1:5369406/1‑71 (MQ=255)

                                               gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACcttct                        <  1:4555267/70‑1 (MQ=255)

                                                 aTGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt                      >  1:2100569/1‑70 (MQ=255)

                                                  tGTCGCCTGATAGGTCGCAGGGTCAATCGGATaga                                                        <  1:741158/35‑1 (MQ=38)

                                                  tGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:2251736/35‑1 (MQ=255)

                                                  tGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                                        <  1:3258420/35‑1 (MQ=255)

                                                  tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGATAGACACCTGCTTc                                        <  1:680537/51‑1 (MQ=255)

                                                  tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        <  1:4658673/51‑1 (MQ=255)

                                                  tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                    <  1:1209712/71‑1 (MQ=255)

                                                  tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                    <  1:1669876/71‑1 (MQ=255)

                                                   gTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                    >  1:4420989/1‑70 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATagaga                                                      <  1:529635/35‑1 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        <  1:1296612/49‑1 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                                        <  1:3872741/49‑1 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg                     >  1:4657609/1‑68 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg                     >  1:4620502/1‑68 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                  >  1:4796847/1‑71 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                  >  1:884865/1‑71 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                  >  1:2910723/1‑71 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                  >  1:4658335/1‑71 (MQ=255)

                                                    tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa                  >  1:4803369/1‑71 (MQ=255)

                                                      gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                 >  1:4009006/1‑70 (MQ=255)

                                                      gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                 <  1:5370912/70‑1 (MQ=255)

                                                      gCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAAt                 <  1:1600183/70‑1 (MQ=255)

                                                        cTGATAGGTCGCGGGATCGATCTGGTACAGGGACTGGCCTGCTTGCACATCGCTGCCTTCAGTGAAAttgc              <  1:1309565/71‑3 (MQ=255)

                                                        cTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc              <  1:2847420/71‑1 (MQ=255)

                                                        cTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc              <  1:2489528/71‑1 (MQ=255)

                                                          gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGa                    <  1:1383000/63‑1 (MQ=255)

                                                          gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc            >  1:1979692/1‑71 (MQ=255)

                                                          gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc            >  1:5075756/1‑71 (MQ=255)

                                                                  gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt    <  1:261144/71‑1 (MQ=255)

                                                                  gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt    <  1:772659/71‑1 (MQ=255)

                                                                    aGGATCAATCTGAAAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGTTTCAGGATaa  <  1:2256437/71‑1 (MQ=255)

ATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAA  >  Exported/313727‑313865

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