BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,421,748	T→C	16.8%	K196K (AAA→AAG)	wzxE ←	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,421,748	0	T	C	16.8%	265.8 / 46.1	113	K196K (AAA→AAG)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base T (58/36); new base C (10/9); total (68/45)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 113 total aligned reads displayed.

CTCAACTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCC  >  minE/2421683‑2421816

cTCAACTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTg                                                                   <  1:1518618/69‑1 (MQ=255)

cTCAACTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGa                                                                  <  1:2122622/70‑1 (MQ=255)

cTCAACTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGa                                                                  <  1:2167917/70‑1 (MQ=255)

cTCAACTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGa                                                                  <  1:3546677/70‑1 (MQ=255)

cTCAACTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGCTTg                                                                   <  1:3384988/69‑1 (MQ=255)

   aaCTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGa                                                                  >  1:4137338/1‑67 (MQ=255)

   aaCTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGa                                                                  >  1:3997079/1‑67 (MQ=255)

   aaCTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCaa                                                              <  1:3732553/71‑1 (MQ=255)

   aaCTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCaa                                                              >  1:3407669/1‑71 (MQ=255)

   aaCTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCaa                                                              >  1:3797161/1‑71 (MQ=255)

   aaCTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCaa                                                              >  1:630300/1‑71 (MQ=255)

   aaCTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCa                                                               <  1:4073720/70‑1 (MQ=255)

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

       gccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAAc                                                             <  1:2610666/68‑1 (MQ=255)

       gccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACat                                                           >  1:771666/1‑70 (MQ=255)

        ccctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatga                                                         >  1:2386564/1‑71 (MQ=255)

         cctgccAGACCGTTATCCCAGCTGTGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgat                                                        <  1:1993553/71‑1 (MQ=255)

          ctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGa                                                                  <  1:1115747/60‑1 (MQ=255)

          ctgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGCTTGATCAACatgatg                                                       >  1:377883/1‑71 (MQ=255)

           tgccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATgg                                                      >  1:3195995/1‑71 (MQ=255)

            gccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatgatg                                                       >  1:3790407/1‑69 (MQ=255)

            gccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatga                                                         >  1:525877/1‑67 (MQ=255)

            gccAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatga                                                         >  1:1865314/1‑67 (MQ=255)

             ccAGACCGTTATCCAAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACatg                                                          <  1:1941437/65‑1 (MQ=255)

              cAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATggcgg                                                   >  1:942598/1‑71 (MQ=255)

                 aCCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACa                                                            >  1:2822683/1‑59 (MQ=255)

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

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

                  ccGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAg                                                >  1:3385171/1‑70 (MQ=255)

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

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

                    gTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCa                                                               <  1:1381681/53‑1 (MQ=255)

                    gTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCa                                                               <  1:2624617/53‑1 (MQ=255)

                     ttATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAAc                                                             <  1:334856/54‑1 (MQ=255)

                     ttATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAAt                                            >  1:3699144/1‑71 (MQ=255)

                      tATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAAc                                                             >  1:1885445/1‑53 (MQ=255)

                      tATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAAc                                                             >  1:3597131/1‑53 (MQ=255)

                      tATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGAAACCACGTTTGATCAAc                                                             >  1:1196785/1‑53 (MQ=255)

                       aTCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAAt                                            >  1:3097171/1‑69 (MQ=255)

                       aTCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTa                                          >  1:3289952/1‑71 (MQ=255)

                       aTCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTa                                          >  1:3046933/1‑71 (MQ=255)

                        tCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTac                                         <  1:2636505/71‑1 (MQ=255)

                        tCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTac                                         <  1:1907525/71‑1 (MQ=255)

                          ccAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTacca                                       <  1:1375321/71‑1 (MQ=255)

                          ccAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTacc                                        <  1:1147011/70‑1 (MQ=255)

                            aGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGCTTGATCAACATGATGGCGGCAGGAATTac                                         <  1:2400960/67‑1 (MQ=255)

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

                                 ggTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAgcgc                                >  1:3884771/1‑71 (MQ=255)

                                 ggTTTCAGATAGCTTAACGGGATGACACCACGCTTGATCAACa                                                            <  1:297312/43‑1 (MQ=255)

                                   tttCAGATAGCTTAACGGGATGACACCACGCTTGATCAACat                                                           <  1:2975091/42‑1 (MQ=255)

                                    ttCAGATAGCTTAACGGGATGACACCACGCTTGATCAACat                                                           <  1:727747/41‑1 (MQ=255)

                                     tCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTacca                                       >  1:1108696/1‑60 (MQ=255)

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

                                     tCAGATAGCTTAACGGGATGACACCACGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCggg                             >  1:3852449/1‑70 (MQ=255)

                                     tCAGATAGCTTAACGGGATGACACCACGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCggg                             >  1:3178670/1‑70 (MQ=255)

                                        gATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATc                         <  1:4077281/71‑1 (MQ=255)

                                         aTAGCTTAACGGGATGACACCACGTTTGATCAACa                                                            >  1:824669/1‑35 (MQ=255)

                                         aTAGCTTAACGGGATGACACCACGTTTGATCAACATGATggcggc                                                  <  1:3233762/45‑1 (MQ=255)

                                         aTAGCTTAACGGGATGACACCACGTTTGATCAACATGATggcggc                                                  <  1:226062/45‑1 (MQ=255)

                                          tAGCTTAACGGGATGACACGACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGaa                           >  1:1333143/1‑67 (MQ=255)

                                          tAGCTTAACGGGATGACACCACGTTTGATCACCatga                                                         >  1:411399/1‑37 (MQ=255)

                                          tAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGaa                           >  1:2616958/1‑67 (MQ=255)

                                          tAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGaa                           >  1:3594550/1‑67 (MQ=255)

                                          tAGCTTAACGGGATGACACCACGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAg                       >  1:2357676/1‑71 (MQ=255)

                                           aggttAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgc                      >  1:3751917/4‑71 (MQ=255)

                                           aGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTaccac                                      >  1:884164/1‑55 (MQ=255)

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

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

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

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

                                           aGCTTAACGGGATGACACCACGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgc                      >  1:147720/1‑71 (MQ=255)

                                            gCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgcg                     >  1:296509/1‑71 (MQ=255)

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

                                               tAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGAGCGGGAATCAgcgc                    >  1:1884368/1‑69 (MQ=255)

                                                aaCGGGATGACACCACGTTTGATCAACATGATGGCGGCAg                                                >  1:1742801/1‑40 (MQ=255)

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

                                                  cGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGAc               <  1:3211440/71‑1 (MQ=255)

                                                   gggATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGAcc              <  1:2802565/71‑1 (MQ=255)

                                                     gATGACACCACGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGAcc              >  1:4046973/1‑69 (MQ=255)

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

                                                        gACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAAt                          <  1:2505743/54‑1 (MQ=255)

                                                        gACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgc                      <  1:835838/58‑1 (MQ=255)

                                                          caccacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATc                         <  1:3842570/53‑1 (MQ=255)

                                                           accacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATc                         >  1:3157899/1‑52 (MQ=255)

                                                           accacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACg      >  1:1046451/1‑71 (MQ=255)

                                                           accacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACg      >  1:826087/1‑71 (MQ=255)

                                                            ccacGTTTGATCAACATGATGGCGGCAGGAATTac                                         <  1:865134/35‑1 (MQ=255)

                                                            ccacGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAgc                      >  1:1160684/1‑54 (MQ=255)

                                                            ccacGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGa                            <  1:3588672/48‑1 (MQ=255)

                                                            ccacGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGc     <  1:3283634/71‑1 (MQ=255)

                                                              acGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATc                         >  1:3119267/1‑49 (MQ=255)

                                                              acGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc   <  1:3048922/71‑1 (MQ=255)

                                                              acGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc   >  1:387364/1‑71 (MQ=255)

                                                              acGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCaa        >  1:525066/1‑66 (MQ=255)

                                                              acGCTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCaa        >  1:3972756/1‑66 (MQ=255)

                                                              acGCTTGATCAACATGATGGCGGCAGGAATTACAACCAGCGCGGGAATCAGCGcc                   >  1:3365570/1‑55 (MQ=255)

                                                               cGTTTGATCAACATGATGGCGGCAGGAATTACCTCCAGCGCGGGAATCAGCGCCAGACCCAGCaa        >  1:2290450/1‑65 (MQ=255)

                                                               cGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCa                  >  1:494414/1‑55 (MQ=255)

                                                               cGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCaa        >  1:706756/1‑65 (MQ=255)

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

CTCAACTGCCCTGCCAGACCGTTATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCC  >  minE/2421683‑2421816

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