BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	313,907	C→T	19.4%	S61N (AGT→AAT)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	313,907	0	C	T	19.4%	149.7 / 27.9	72	S61N (AGT→AAT)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base C (25/33); new base T (12/2); total (37/35)
Fisher's exact test for biased strand distribution p-value = 6.18e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAACGGCGG  >  minE/313838‑313976

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

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

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

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

     tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                                                                 >  1:513687/1‑71 (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:3001702/1‑66 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                  tAACTTGAGGAGGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                             <  1:4120954/62‑1 (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)

                                        gAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGAt                                           >  1:4345115/1‑58 (MQ=255)

                                           gACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag                                     <  1:2579096/61‑1 (MQ=255)

                                                    cTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGt                   <  1:3090748/70‑1 (MQ=255)

                                                      gCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGtttt                         <  1:1925149/62‑1 (MQ=255)

                                                       cGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTtac                <  1:4439050/70‑1 (MQ=255)

                                                            cGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGcc          <  1:4939065/71‑1 (MQ=255)

                                                            cGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGcc          <  1:3288105/71‑1 (MQ=255)

                                                                 ggCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGtt                           <  1:108045/49‑1 (MQ=255)

                                                                  gCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTc                               <  1:3546895/44‑1 (MQ=255)

                                                                  gCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAAcggc    <  1:1012473/71‑1 (MQ=255)

                                                                  gCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAAcggc    <  1:1390951/71‑1 (MQ=255)

                                                                    aCTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCaa        >  1:3745896/1‑65 (MQ=255)

                                                                    aCTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAAcggcgg  >  1:3806753/1‑71 (MQ=255)

                                                                    aCTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAAcggcgg  >  1:3998447/1‑71 (MQ=255)

                                                                    aCTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAAcggcgg  >  1:4151366/1‑71 (MQ=255)

                                                                    aCTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAAcggcgg  >  1:4838113/1‑71 (MQ=255)

                                                                    aCTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAAcggcgg  >  1:4494780/1‑71 (MQ=255)

CCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCCAACGGCGG  >  minE/313838‑313976

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