BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	313,876	T→C	19.2%	Q71Q (CAA→CAG)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	313,876	0	T	C	19.2%	201.9 / 41.1	99	Q71Q (CAA→CAG)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base T (39/41); new base C (13/6); total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 102 total aligned reads displayed.

ATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAG  >  Exported/313807‑313943

aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACtt                                                                    <  1:4302614/71‑1 (MQ=255)

aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACt                                                                     <  1:4365912/70‑1 (MQ=255)

    agagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAgg                                                                <  1:2898845/71‑1 (MQ=255)

          acacCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                                          <  1:5298/71‑1 (MQ=255)

             ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                                                   <  1:755252/59‑1 (MQ=255)

             ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                                                   <  1:320484/59‑1 (MQ=255)

             ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACg                                                             <  1:3472635/65‑1 (MQ=255)

               tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGa                                                               >  1:1226219/1‑61 (MQ=255)

               tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTc                                                       <  1:268844/69‑1 (MQ=255)

               tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTc                                                       <  1:390887/69‑1 (MQ=255)

                gCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGaa                                                           >  1:3313676/1‑64 (MQ=255)

                gCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGaa                                                           >  1:5434545/1‑64 (MQ=255)

                 cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                                                     >  1:3209556/1‑69 (MQ=255)

                 cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                                                     >  1:4422835/1‑69 (MQ=255)

                 cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                                                     >  1:1906437/1‑69 (MQ=255)

                    cGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATc                                                >  1:4141687/1‑71 (MQ=255)

                      acatCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACt                                                         <  1:4917131/57‑1 (MQ=255)

                      acatCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACt                                                         <  1:5261606/57‑1 (MQ=255)

                      aTGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAg                                                                 <  1:3123355/52‑1 (MQ=255)

                       tGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                                             <  1:4271291/71‑1 (MQ=255)

                        gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                                             >  1:4521836/1‑70 (MQ=255)

                        atCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGc                                                    <  1:3393494/62‑1 (MQ=255)

                         tCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACgat                                            >  1:2845283/1‑68 (MQ=255)

                           gCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                                        <  1:2650266/56‑1 (MQ=255)

                           gCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAgg                                         >  1:4051084/1‑71 (MQ=255)

                            cTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGc                                        <  1:254615/71‑1 (MQ=255)

                              aCCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGt                                            <  1:5230191/65‑1 (MQ=255)

                              aCCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGt                                            <  1:1228423/65‑1 (MQ=255)

                               ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                                          <  1:5564794/50‑1 (MQ=255)

                               ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                                          <  1:2626268/50‑1 (MQ=255)

                               ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt                                     >  1:1869642/1‑71 (MQ=255)

                               ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt                                     >  1:1876190/1‑71 (MQ=255)

                               ccTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCatt                                     <  1:2516331/71‑3 (MQ=255)

                                 ttcttTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                                                   <  1:5387738/39‑1 (MQ=255)

                                  tcttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACt                                                         >  1:1237771/1‑48 (MQ=255)

                                   cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                                        >  1:3856700/1‑48 (MQ=255)

                                   cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                                        >  1:4605907/1‑48 (MQ=255)

                                   cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTg                                 >  1:3984667/1‑71 (MQ=255)

                                    tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                                >  1:2760457/1‑71 (MQ=255)

                                    tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                                >  1:513687/1‑71 (MQ=255)

                                     ttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATcc                                               <  1:2144178/55‑1 (MQ=255)

                                     gtGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCTGCTATACGATAAGCATTGGTGCg                               <  1:2343094/70‑1 (MQ=255)

                                      tGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCgg                              >  1:1375392/1‑70 (MQ=255)

                                       gAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGa                              <  1:5551224/70‑1 (MQ=255)

                                        aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGt                                  >  1:1981234/1‑65 (MQ=255)

                                        aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                                 >  1:536100/1‑66 (MQ=255)

                                        aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                                 >  1:3001702/1‑66 (MQ=255)

                                        aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                                 >  1:1572776/1‑66 (MQ=255)

                                        aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCggc                             >  1:4205515/1‑68 (MQ=255)

                                        aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc                            <  1:4455268/71‑1 (MQ=255)

                                        aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc                            >  1:2215143/1‑71 (MQ=255)

                                        aaaTTACTCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGAc                             <  1:5693475/70‑1 (MQ=255)

                                        aaaTTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGAcc                            >  1:2856905/1‑71 (MQ=255)

                                          aTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                    >  1:5135467/1‑61 (MQ=255)

                                          aTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                    >  1:5356662/1‑61 (MQ=255)

                                          aTTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCg                          <  1:2536126/71‑1 (MQ=255)

                                           ttACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCgg                         <  1:4605918/71‑1 (MQ=255)

                                           ttACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCgg                         <  1:2171877/71‑1 (MQ=255)

                                             aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                                              >  1:4708280/1‑48 (MQ=255)

                                             aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                                              >  1:4701667/1‑48 (MQ=255)

                                             aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                                              >  1:1544449/1‑48 (MQ=255)

                                             aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGa                        <  1:2308429/70‑1 (MQ=255)

                                               gCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGc                     >  1:139169/1‑71 (MQ=255)

                                                cTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCt                    <  1:1138572/71‑1 (MQ=255)

                                                  tCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                                             <  1:2520358/44‑1 (MQ=255)

                                                  tCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTc                   >  1:4285273/1‑70 (MQ=255)

                                                   cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCgg                 >  1:723581/1‑71 (MQ=255)

                                                   cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCgg                 >  1:3156673/1‑71 (MQ=255)

                                                   cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCg                  <  1:5147570/70‑1 (MQ=255)

                                                    aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTa                                           >  1:1272766/1‑44 (MQ=255)

                                                    aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                                    >  1:887874/1‑51 (MQ=255)

                                                    aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                                    >  1:1081018/1‑51 (MQ=255)

                                                    aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                                    >  1:2580969/1‑51 (MQ=255)

                                                     ggATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGtt               >  1:1133632/1‑71 (MQ=255)

                                                      gATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtg              <  1:663073/71‑1 (MQ=255)

                                                        tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                                       >  1:5552838/1‑44 (MQ=255)

                                                        tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                                       >  1:480781/1‑44 (MQ=255)

                                                        tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:157456/71‑1 (MQ=255)

                                                        tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:4601654/71‑1 (MQ=255)

                                                        tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:4781529/71‑1 (MQ=255)

                                                        tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:2803890/71‑1 (MQ=255)

                                                          aTCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGa           <  1:28518/70‑1 (MQ=255)

                                                          aTCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                    >  1:1602933/1‑45 (MQ=255)

                                                          aTCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                    >  1:3751927/1‑45 (MQ=255)

                                                          aTCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                    >  1:180262/1‑45 (MQ=255)

                                                           tCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTa                                           >  1:4622318/1‑37 (MQ=255)

                                                             ccGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTg       >  1:4630114/1‑71 (MQ=255)

                                                               gCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCa     <  1:5687358/71‑1 (MQ=255)

                                                                cTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag    <  1:2953582/71‑1 (MQ=255)

                                                                cTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag    <  1:2961360/71‑1 (MQ=255)

                                                                cTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag    <  1:4744225/71‑1 (MQ=255)

                                                                 tAACTTGAGGAGGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:4120954/62‑1 (MQ=255)

                                                                 tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                                       >  1:2529336/1‑35 (MQ=255)

                                                                 tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGtt               >  1:464566/1‑59 (MQ=255)

                                                                 tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:1502636/62‑1 (MQ=255)

                                                                 tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:2541091/62‑1 (MQ=255)

                                                                  aaCTTGAGGACGAACTTCTGCGGTCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCagag  >  1:5040791/1‑71 (MQ=255)

                                                                  aaCTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:2107542/61‑1 (MQ=255)

                                                                    cTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGAt          >  1:4652778/1‑61 (MQ=255)

                                                                     ttGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg            <  1:1048885/58‑1 (MQ=255)

ATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAG  >  Exported/313807‑313943

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