BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,576	C→A	5.9%	L171L (CTG→CTT)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,576	0	C	A	5.9%	331.0 / 4.1	101	L171L (CTG→CTT)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base C (53/42); new base A (2/4); total (55/46)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 105 total aligned reads displayed.

AATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCA  >  minE/313508‑313645

aaTGCGGCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAg                                                                     <  1:4957902/71‑1 (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)

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

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

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

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

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

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

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

           cGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg                                                          <  1:2819622/71‑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)

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

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

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

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

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

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

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

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

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

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

                          aCCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCaa                                                 >  1:1506094/1‑65 (MQ=255)

                          aCCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCaa                                                 >  1:2968328/1‑65 (MQ=255)

                          aCCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAg                                           >  1:1059598/1‑71 (MQ=255)

                            cAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCt                                         <  1:3454192/71‑1 (MQ=255)

                             aaTGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCtt                                        <  1:5120344/71‑1 (MQ=255)

                              aTGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTc                                                   <  1:4571945/59‑1 (MQ=255)

                              aTGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTc                                                   <  1:1477636/59‑1 (MQ=255)

                              aTGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTc                                                   <  1:5623541/59‑1 (MQ=255)

                               tGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCcgcg                                                          >  1:1829267/1‑51 (MQ=255)

                               tGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCaa                                                 >  1:3440844/1‑60 (MQ=255)

                               tGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTc                                      <  1:4584586/71‑1 (MQ=255)

                                gCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCaa                                                 >  1:2261863/1‑59 (MQ=255)

                                gCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCaa                                                 >  1:712816/1‑59 (MQ=255)

                                gCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcg                                     >  1:5637393/1‑71 (MQ=255)

                                 cGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGa                                                                    >  1:3357281/1‑39 (MQ=255)

                                 cGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGa                                                                    >  1:5145249/1‑39 (MQ=255)

                                 cGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCa                                            >  1:3694373/1‑63 (MQ=255)

                                   aCCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCa                                                        <  1:2506763/49‑1 (MQ=255)

                                   aCCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCa                                                        <  1:823128/49‑1 (MQ=255)

                                   aCCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCa                                                        <  1:4850801/49‑1 (MQ=255)

                                   aCCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccg                                  <  1:3844898/71‑1 (MQ=255)

                                    ccGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccgc                                 <  1:4261138/71‑1 (MQ=255)

                                    ccGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccg                                  >  1:159968/1‑70 (MQ=255)

                                     cGCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGATGCGAGCGCTTTCGACTGTGGCTTTtgcggc                                 <  1:3684844/70‑1 (MQ=255)

                                     cGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCa                                <  1:137215/71‑1 (MQ=255)

                                      gCTAATTGGCGCAGTGACTTTGGTATAAGCAAGATTGATGCGAGCGCTTTCGACTGTGGCTTTTGCGGCaa                               <  1:2618231/71‑2 (MQ=255)

                                      gCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccgc                                 <  1:5635521/69‑1 (MQ=255)

                                      gCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccgc                                 <  1:5641863/69‑1 (MQ=255)

                                      gCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgcc                                   <  1:4520171/67‑1 (MQ=255)

                                        tAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccgc                                 >  1:4573751/1‑67 (MQ=255)

                                        tAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccgc                                 >  1:3958319/1‑67 (MQ=255)

                                        tAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGt                              <  1:4904374/70‑1 (MQ=255)

                                          aTCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGtt                             <  1:700084/69‑1 (MQ=255)

                                          aTCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGtt                             <  1:2694304/69‑1 (MQ=255)

                                           tCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAg                                           <  1:418016/54‑1 (MQ=255)

                                           tCGGAGAGGTGACTTAGGTGTAAGCCAGATTGATCCGCGCAGTTTCaa                                                 >  1:2007864/1‑48 (MQ=255)

                                            cGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcg                                     >  1:4558197/1‑59 (MQ=255)

                                             ggAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTa                            >  1:523444/1‑67 (MQ=255)

                                             ggAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTa                            >  1:2674627/1‑67 (MQ=255)

                                               atagGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCCCCGCAGTTACCgcag                      <  1:4828774/69‑1 (MQ=255)

                                               agagGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCgcag                      <  1:4566644/71‑1 (MQ=255)

                                               agagGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCgcag                      <  1:5159677/71‑1 (MQ=255)

                                               agagGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCgcag                      <  1:2974496/71‑1 (MQ=255)

                                               agagGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCgcag                      <  1:4123146/71‑1 (MQ=255)

                                                gagGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCgcag                      >  1:2062819/1‑70 (MQ=255)

                                                   gTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGttt                                                    <  1:3608245/37‑1 (MQ=255)

                                                   gTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGttt                                                    <  1:1053314/37‑1 (MQ=255)

                                                   gTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTc                                                   <  1:1313283/38‑1 (MQ=255)

                                                     gACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCg                >  1:5040767/1‑71 (MQ=255)

                                                     gACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCg                >  1:3447020/1‑71 (MQ=255)

                                                     gACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCg                >  1:4287956/1‑71 (MQ=255)

                                                      aCTTTGGTGTAAGCCAGATTGATCCTCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCg                <  1:3390142/70‑1 (MQ=255)

                                                          tgctgtAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGt           >  1:367355/4‑71 (MQ=255)

                                                          tGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCttt                                       >  1:4177671/1‑43 (MQ=255)

                                                          tGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCttt                                       >  1:4552531/1‑43 (MQ=255)

                                                          tGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGt           >  1:1085495/1‑71 (MQ=255)

                                                           ggTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGtt          >  1:2886451/1‑71 (MQ=255)

                                                             tgtAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGt                              <  1:3618412/49‑1 (MQ=255)

                                                               tAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccgc                                 >  1:4470734/1‑44 (MQ=255)

                                                               tAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgccgc                                 >  1:1534975/1‑44 (MQ=255)

                                                               tAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTgcgc      <  1:4792377/71‑1 (MQ=255)

                                                                 aGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgc                                    >  1:1604172/1‑39 (MQ=255)

                                                                 aGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgc                                    >  1:232424/1‑39 (MQ=255)

                                                                 aGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTcgc                                    >  1:3572066/1‑39 (MQ=255)

                                                                 aGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCAt    >  1:140671/1‑71 (MQ=255)

                                                                   ccAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCa  <  1:4609226/71‑1 (MQ=255)

                                                                   ccAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCa  <  1:2526721/71‑1 (MQ=255)

                                                                    cAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCa  <  1:1269102/70‑1 (MQ=255)

AATGCGCCTTCCGTCACGTTCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCA  >  minE/313508‑313645

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