BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,421,757	C→T	17.8%	M193I (ATG→ATA)	wzxE ←	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,421,757	0	C	T	17.8%	294.4 / 50.3	129	M193I (ATG→ATA)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base C (55/50); new base T (9/14); total (65/64)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 130 total aligned reads displayed.

CTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAC  >  minE/2421688‑2421824

ctgccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACa                                                                    >  1:2509731/1‑71 (MQ=255)

 tgccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACat                                                                   >  1:2049725/1‑71 (MQ=255)

   ccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATatg                                                                  <  1:3924542/70‑1 (MQ=255)

   ccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatg                                                                  <  1:2365830/70‑1 (MQ=255)

   ccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACat                                                                   <  1:1104021/69‑1 (MQ=255)

   ccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACat                                                                   <  1:714189/69‑1 (MQ=255)

    cctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTTATCAACatgat                                                                <  1:2391947/71‑1 (MQ=255)

    cctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATatgat                                                                >  1:2081010/1‑71 (MQ=255)

     ctgccAGGCCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgatg                                                               >  1:3200476/1‑71 (MQ=255)

     ctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgatg                                                               <  1:2666403/71‑1 (MQ=255)

     ctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgatg                                                               >  1:1176283/1‑71 (MQ=255)

     ctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgatg                                                               <  1:3133496/71‑1 (MQ=255)

      tgccAGCCCGTTATCCCAGCTGGGTTTCAGATAGCTTACCGGGATGACACCACGTTTGATCAACatgatg                                                               >  1:913931/1‑70 (MQ=255)

      tgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGCCACCACGTTTGATCAATATGATgg                                                              >  1:2260375/1‑71 (MQ=255)

      tgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATatgat                                                                >  1:1242547/1‑69 (MQ=255)

         cAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATATGATggcgg                                                           >  1:796655/1‑71 (MQ=255)

         cAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATATGATggcgg                                                           >  1:2247932/1‑71 (MQ=255)

          aGACCGTTATCCCAGCTGGGTTACAGATAGCTTAACGGGATGACACCACGTTTGATCAATATGATggcgg                                                           >  1:174144/1‑70 (MQ=255)

           gACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCa                                                         <  1:2520549/71‑1 (MQ=255)

             ccGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATATGATGGCGGCAgg                                                       <  1:2184520/71‑1 (MQ=255)

             ccGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAgg                                                       <  1:1157941/71‑1 (MQ=255)

              cGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAgg                                                       >  1:1129579/1‑70 (MQ=255)

               gTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGaa                                                     <  1:413013/71‑1 (MQ=255)

               gTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGaa                                                     <  1:1626905/71‑1 (MQ=255)

                ttATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATATGATGGCGGCAGGAAt                                                    <  1:291802/71‑1 (MQ=255)

                ttATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgat                                                                >  1:3296436/1‑59 (MQ=255)

                ttATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAAt                                                    <  1:3157047/71‑1 (MQ=255)

                ttATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAAt                                                    <  1:1882454/71‑1 (MQ=255)

                ttATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAAt                                                    <  1:650781/71‑1 (MQ=255)

                 tATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATgg                                                              >  1:2067671/1‑60 (MQ=255)

                  aTCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACat                                                                   >  1:1290978/1‑54 (MQ=255)

                  aTCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATggcg                                                            >  1:1160744/1‑61 (MQ=255)

                  aTCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATggc                                                             >  1:2107904/1‑60 (MQ=255)

                  aTCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAAtt                                                   <  1:1443528/70‑1 (MQ=255)

                   tCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAATATGATGGCGGCAGGAATTac                                                 >  1:704674/1‑71 (MQ=255)

                   tCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgat                                                                <  1:2189816/56‑1 (MQ=255)

                    cccAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTac                                                 >  1:3266251/1‑70 (MQ=255)

                    accAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTac                                                 >  1:2258183/2‑70 (MQ=255)

                     ccAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAgg                                                       <  1:3809871/63‑1 (MQ=255)

                     ccAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTacca                                               >  1:956548/1‑71 (MQ=255)

                      cAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTaccac                                              >  1:2957477/1‑71 (MQ=255)

                       aGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTaccacc                                             <  1:2776151/71‑1 (MQ=255)

                         cTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAg                                           >  1:3511740/1‑71 (MQ=255)

                           gggTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCa                                                         >  1:3696721/1‑55 (MQ=255)

                           gggTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAgg                                                       <  1:3861077/57‑1 (MQ=255)

                           gggTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCcgcg                                         >  1:3394890/1‑71 (MQ=255)

                           gggTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAg                                           <  1:3220810/69‑1 (MQ=255)

                            ggTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAgcg                                         <  1:3542778/70‑1 (MQ=255)

                             gTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGcgcg                                       <  1:1477631/71‑1 (MQ=255)

                              tttCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTaccac                                              <  1:1120805/63‑1 (MQ=255)

                                tCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATggcggc                                                          >  1:1887727/1‑49 (MQ=255)

                                tCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCggg                                     >  1:215686/1‑70 (MQ=255)

                                 cAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGa                                    <  1:3285998/70‑1 (MQ=255)

                                 cAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGa                                    <  1:1114327/70‑1 (MQ=255)

                                  aGATAGCTTAACGGGATGACACCACGTTTGATCAATATGATGGCGGCAGGaa                                                     <  1:2819490/52‑1 (MQ=255)

                                  aGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAAt                                  >  1:3253820/1‑71 (MQ=255)

                                  aGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAAt                                  <  1:3165903/71‑1 (MQ=255)

                                    aTAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTacca                                               >  1:2678015/1‑56 (MQ=255)

                                    aTAGCTTAAAGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAAt                                  >  1:742496/1‑69 (MQ=255)

                                     tAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCa                                <  1:3356590/70‑1 (MQ=255)

                                     tAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCa                                <  1:82728/70‑1 (MQ=255)

                                      aGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgc                              >  1:1258453/1‑71 (MQ=255)

                                      aGCTTAACGGGATGACACCACGTTTGATCAA‑A‑‑ATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgcgc                            >  1:564909/1‑70 (MQ=255)

                                       gCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAgg                                                       >  1:196624/1‑45 (MQ=255)

                                       gCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAg                                                        <  1:1677612/44‑1 (MQ=255)

                                       gCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgc                              <  1:3579524/70‑1 (MQ=255)

                                         ttAACGGGATGACACCACGTTTGATCAACATGATGGCGGCa                                                         <  1:3084701/41‑1 (MQ=255)

                                         ttAACGGGATGACACCACGTTTGATCAACATGATGGCGGCa                                                         <  1:3096779/41‑1 (MQ=255)

                                         ttAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgcgc                            >  1:3168145/1‑70 (MQ=255)

                                          tAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgcgc                            >  1:1859958/1‑69 (MQ=255)

                                           aaCGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAgcg                                         >  1:3665563/1‑55 (MQ=255)

                                           aaCGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgcgc                            >  1:543618/1‑68 (MQ=255)

                                           aaCGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgcgc                            >  1:42679/1‑68 (MQ=255)

                                           aaCGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCa                          <  1:230995/70‑1 (MQ=255)

                                             cGGGATGACACCACGTTTGATCAATATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGAc                       <  1:428002/71‑1 (MQ=255)

                                              gggATGACACCACGTTTGATCAATATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGAcc                      <  1:2422058/71‑1 (MQ=255)

                                              gggATGACACCACGTTTGATCAATATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGAcc                      <  1:2765753/71‑1 (MQ=255)

                                              gggATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGGcc                      <  1:174291/71‑1 (MQ=255)

                                               ggATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTAccacca                                            >  1:2808306/1‑48 (MQ=255)

                                                 aTGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAg                   >  1:3919747/1‑71 (MQ=255)

                                                 aTGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAg                   >  1:2167611/1‑71 (MQ=255)

                                                  tGACACAACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAg                   >  1:505597/1‑70 (MQ=255)

                                                    acacCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCa                          <  1:1517653/61‑1 (MQ=255)

                                                    acacCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCaa                >  1:2009477/1‑71 (MQ=255)

                                                     caccacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCC‑‑ACCCAGCAACGc             <  1:2413810/71‑1 (MQ=255)

                                                       ccacGTTTGATCAATATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCAAGCAACGc             >  1:2038618/1‑71 (MQ=255)

                                                       ccacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGc             <  1:9424/71‑1 (MQ=255)

                                                        cacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACg              >  1:1630504/1‑69 (MQ=255)

                                                        cacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGcc            <  1:3755684/71‑1 (MQ=255)

                                                         acGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGc             >  1:1556169/1‑69 (MQ=255)

                                                          cGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGcccc          >  1:3010671/1‑71 (MQ=255)

                                                          cGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGcccc          >  1:3931241/1‑71 (MQ=255)

                                                                aTCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCAt     <  1:1728716/70‑1 (MQ=255)

                                                                aTCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCAt     <  1:2665489/70‑1 (MQ=255)

                                                                 tCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGa                                    >  1:3830099/1‑38 (MQ=255)

                                                                 tCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATaa   >  1:902320/1‑71 (MQ=255)

                                                                 tCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATaa   <  1:2028589/71‑1 (MQ=255)

                                                                   aaCATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAg                   >  1:3688910/1‑53 (MQ=255)

                                                                   aaCATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAc  <  1:2104517/70‑1 (MQ=255)

                                                                   aaCATGATGGCGGCAGGAAATACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAc  >  1:2043185/1‑70 (MQ=255)

CTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAC  >  minE/2421688‑2421824

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