BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,531	C→T	7.2%	K186K (AAG→AAA)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,531	0	C	T	7.2%	286.5 / 8.0	97	K186K (AAG→AAA)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base C (56/34); new base T (2/5); total (58/39)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 101 total aligned reads displayed.

AAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAAC  >  minE/313462‑313599

aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGAc                                                                      >  1:2868680/1‑70 (MQ=255)

aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGAc                                                                      >  1:4574650/1‑70 (MQ=255)

aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                                     >  1:3475514/1‑71 (MQ=255)

aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                                     >  1:1856073/1‑71 (MQ=255)

aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                                     >  1:1450698/1‑71 (MQ=255)

aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                                     >  1:2749826/1‑71 (MQ=255)

aaGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                                     >  1:4640144/1‑71 (MQ=255)

aaGTTGCTGCACGGGTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACt                                                                     <  1:1732072/71‑1 (MQ=255)

     gctgcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTAcc                                                                 <  1:540531/70‑1 (MQ=255)

     gctgcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTAcc                                                                 <  1:3693696/70‑1 (MQ=255)

        gcACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAAt                                                              >  1:4312313/1‑70 (MQ=255)

         cACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                            >  1:5399426/1‑71 (MQ=255)

          aCGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                            >  1:1981536/1‑70 (MQ=255)

          aCGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCg                                                           <  1:1841616/71‑1 (MQ=255)

            ggTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGAc                                                         >  1:1172146/1‑71 (MQ=255)

             gTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                            <  1:1699255/67‑1 (MQ=255)

             gTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                            <  1:3544861/67‑1 (MQ=255)

               tGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGa                                                          >  1:2772511/1‑67 (MQ=255)

               tGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGa                                                          >  1:2442808/1‑67 (MQ=255)

               tGCCAGCGCAGTCGCCTGACCGCTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGa                                                          >  1:3024802/1‑67 (MQ=255)

                 ccAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCaa                                                               >  1:3156769/1‑60 (MQ=255)

                 ccAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTa                                                    <  1:5544201/71‑1 (MQ=255)

                   aGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAAt                                                  >  1:3808787/1‑71 (MQ=255)

                    gcgcAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATc                                                 >  1:5354033/1‑71 (MQ=255)

                      gcAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATc                                                 >  1:3150097/1‑69 (MQ=255)

                      gcAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCgg                                               <  1:2841279/71‑1 (MQ=255)

                       cAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                              <  1:2631577/71‑1 (MQ=255)

                       cAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGA‑CGCTAATCGga                                              <  1:4390026/70‑1 (MQ=255)

                         gTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                              >  1:1120198/1‑69 (MQ=255)

                            gCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAgg                                          <  1:4347586/70‑1 (MQ=255)

                                 aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCt                                                     <  1:5126496/54‑1 (MQ=255)

                                 aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGa                                       >  1:1387643/1‑68 (MQ=255)

                                 aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACtt                                    >  1:1000774/1‑71 (MQ=255)

                                 aCCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACt                                     <  1:4078213/70‑1 (MQ=255)

                                  ccGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCgg                                               <  1:2058600/59‑1 (MQ=255)

                                  ccGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCgg                                               <  1:2773031/59‑1 (MQ=255)

                                   cGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGc                                                            >  1:2714395/1‑45 (MQ=255)

                                   cGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTaa                                                   >  1:1363686/1‑54 (MQ=255)

                                    gTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGgag                                             >  1:2922155/1‑59 (MQ=255)

                                     ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACttt                                   <  1:1507394/68‑1 (MQ=255)

                                     ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACttt                                   <  1:3908390/68‑1 (MQ=255)

                                     ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgt                                >  1:1044453/1‑71 (MQ=255)

                                     ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgt                                >  1:1079468/1‑71 (MQ=255)

                                     ttCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgt                                >  1:4168834/1‑71 (MQ=255)

                                      tCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                               >  1:284897/1‑71 (MQ=255)

                                      tCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                               >  1:278007/1‑71 (MQ=255)

                                       cTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                               >  1:3060659/1‑70 (MQ=255)

                                          tACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtgt                              <  1:3958640/68‑1 (MQ=255)

                                          tACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTa                             <  1:1234323/69‑1 (MQ=255)

                                          tACCAATGCGCCTTCCGTCACGTGCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtgt                              <  1:2560667/68‑1 (MQ=255)

                                           aCCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGc                          >  1:4053291/1‑71 (MQ=255)

                                              aaTGCGGCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAg                       <  1:4957902/71‑1 (MQ=255)

                                                tGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                               >  1:1280405/1‑61 (MQ=255)

                                                tGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtg                               >  1:2424949/1‑61 (MQ=255)

                                                 gcgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                              >  1:877363/1‑45 (MQ=255)

                                                  cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                              >  1:609626/1‑44 (MQ=255)

                                                  cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGga                                              >  1:885729/1‑44 (MQ=255)

                                                  cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt                    >  1:5604312/1‑70 (MQ=255)

                                                  cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt                    >  1:1815974/1‑70 (MQ=255)

                                                  cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt                    <  1:555686/70‑1 (MQ=255)

                                                  cgcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTg                   <  1:3291585/71‑1 (MQ=255)

                                                   gcCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTg                   <  1:345774/70‑1 (MQ=255)

                                                     cTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATc                >  1:5429045/1‑71 (MQ=255)

                                                      ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGc                          <  1:4492323/60‑1 (MQ=255)

                                                      ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGc                          <  1:5609237/60‑1 (MQ=255)

                                                      ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGc                          <  1:410310/60‑1 (MQ=255)

                                                      ttCGGTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGAt                 <  1:2698399/69‑1 (MQ=255)

                                                      ttCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAg                       <  1:333139/63‑1 (MQ=255)

                                                      ttCCGTCACGTTCGACTTACCAATACGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATc                >  1:1768711/1‑70 (MQ=255)

                                                       tcggTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGa                  >  1:1674694/4‑67 (MQ=255)

                                                       tcggTCACAGTCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGa                  >  1:2488414/4‑67 (MQ=255)

                                                        ccGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgc             >  1:747901/1‑71 (MQ=255)

                                                        ccGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgc             >  1:5309105/1‑71 (MQ=255)

                                                         cGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg            >  1:3528495/1‑71 (MQ=255)

                                                         cGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg            <  1:2819622/71‑1 (MQ=255)

                                                         cGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg            <  1:2231851/71‑1 (MQ=255)

                                                          gTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTg                                  <  1:4627500/48‑1 (MQ=255)

                                                          gTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTg                                  <  1:1119593/48‑1 (MQ=255)

                                                           tCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGgtgt                              <  1:5474630/51‑1 (MQ=255)

                                                           tCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCa          <  1:2370381/71‑1 (MQ=255)

                                                            cACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGa                      >  1:105230/1‑58 (MQ=255)

                                                            cACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAg         >  1:1840239/1‑71 (MQ=255)

                                                             aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTg                   >  1:3201007/1‑60 (MQ=255)

                                                             aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGAt                 >  1:5252971/1‑62 (MQ=255)

                                                             aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGt        >  1:3249714/1‑71 (MQ=255)

                                                             aCGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACATTGGTGTAAGcc                         >  1:378175/1‑54 (MQ=255)

                                                              cGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACtt                                    <  1:5043541/42‑1 (MQ=255)

                                                                ttCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGc                          <  1:4100273/50‑1 (MQ=255)

                                                                ttCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTc     <  1:5258923/71‑1 (MQ=255)

                                                                 tCGATTTGCCGATACGTCCGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGAtt                    <  1:905778/55‑1 (MQ=255)

                                                                 tCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt                    >  1:3766842/1‑55 (MQ=255)

                                                                 tCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt                    >  1:3422879/1‑55 (MQ=255)

                                                                 tCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt                    >  1:3141002/1‑55 (MQ=255)

                                                                 tCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGAtt                    >  1:2593473/1‑55 (MQ=255)

                                                                 tCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGa                  >  1:3952600/1‑57 (MQ=255)

                                                                  cGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTg                                  >  1:5520484/1‑40 (MQ=255)

                                                                  cGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTg                                  >  1:28298/1‑40 (MQ=255)

                                                                  cGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCaa   >  1:717827/1‑71 (MQ=255)

                                                                     cTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCa                        <  1:3053031/47‑1 (MQ=255)

                                                                     cTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCTCGCAGTTTCAAc  <  1:3527718/69‑1 (MQ=255)

AAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACCGTTCTGTACCAATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAAC  >  minE/313462‑313599

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