BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	313,861	G→T	17.2%	I76I (ATC→ATA) ‡	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	313,861	0	G	T	17.2%	246.2 / 27.6	93	I76I (ATC→ATA) ‡	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base G (42/35); new base T (10/6); total (52/41)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

CGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTG  >  Exported/313792‑313931

cGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGa                                                                       >  1:895848/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)

       tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATccc                                                                >  1:3630095/1‑71 (MQ=255)

       tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATcc                                                                 >  1:88970/1‑70 (MQ=255)

             tGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAAc                                                          >  1:5186021/1‑71 (MQ=255)

               aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTa                                                            <  1:472724/67‑1 (MQ=255)

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

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

                  gagagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAg                                                     <  1:3355624/71‑1 (MQ=255)

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

                       agaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAAt                                                                   >  1:4562670/1‑52 (MQ=255)

                       agaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                              >  1:3336319/1‑57 (MQ=255)

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

                            ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt                                                                      <  1:688721/44‑1 (MQ=255)

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

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

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

                              tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                    >  1:2894055/1‑44 (MQ=255)

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

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

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

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

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

                                cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                    >  1:403855/1‑42 (MQ=255)

                                cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                    >  1:5225794/1‑42 (MQ=255)

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

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

                                cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                                         >  1:4422835/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)

                                       gTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                                    >  1:1187704/1‑35 (MQ=255)

                                       gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                              >  1:4344865/1‑41 (MQ=255)

                                       gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                              >  1:4953699/1‑41 (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)

                                              ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCt                                                             <  1:991978/35‑1 (MQ=255)

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

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

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

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

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

                                                ttcttTGAAATTACGCTTCAGGATAATCCCGCTaa                                                           >  1:3772504/1‑35 (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: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: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:3156673/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

CGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTG  >  Exported/313792‑313931

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