BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,833	G→A	14.5%	L246L (CTG→CTA)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,833	0	G	A	14.5%	252.6 / 28.3	103	L246L (CTG→CTA)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (41/47); new base A (9/6); total (50/53)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 105 total aligned reads displayed.

CGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCT  >  minE/2570766‑2570901

cgcCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGt                                                                   >  1:2187208/1‑71 (MQ=255)

cgcCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGt                                                                   >  1:4854762/1‑71 (MQ=255)

 gcagAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                                                  >  1:5095792/4‑71 (MQ=255)

 gcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                                                  >  1:5027913/1‑71 (MQ=255)

    gAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                                                               <  1:4523482/71‑1 (MQ=255)

    gAGCAATGGCTACGTTACTCAGGTACTGAACCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTcc                                                                 >  1:2462052/1‑69 (MQ=255)

      gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                                                               <  1:3771478/69‑1 (MQ=255)

      gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                                                               <  1:2791965/69‑1 (MQ=255)

         aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg                                                          >  1:544441/1‑71 (MQ=255)

         aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg                                                          <  1:4950043/71‑1 (MQ=255)

         aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                                                  <  1:5461260/63‑1 (MQ=255)

         aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                                                  <  1:2249279/63‑1 (MQ=255)

          tGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATgg                                                         >  1:622311/1‑71 (MQ=255)

          tGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa                                                            <  1:198856/68‑1 (MQ=255)

           ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                                                >  1:4733665/1‑63 (MQ=255)

           ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                                                >  1:951436/1‑63 (MQ=255)

           ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATgg                                                         >  1:2486290/1‑70 (MQ=255)

           ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                                                        >  1:4296935/1‑71 (MQ=255)

           ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                                                        >  1:5129572/1‑71 (MQ=255)

            gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                                                       <  1:2328269/71‑1 (MQ=255)

            gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                                                       >  1:3392954/1‑71 (MQ=255)

             cTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCt                                                      <  1:4172809/71‑1 (MQ=255)

                cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                                                     <  1:1319461/69‑1 (MQ=255)

                cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    <  1:4657486/70‑1 (MQ=255)

                cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCa                                                   >  1:2203160/1‑71 (MQ=255)

                  ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                                                  <  1:2937265/54‑1 (MQ=255)

                  ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATc                                                 >  1:508147/1‑71 (MQ=255)

                    aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                                                                      >  1:2121406/1‑48 (MQ=255)

                    aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                                                                      >  1:2601556/1‑48 (MQ=255)

                    aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATcc                                                >  1:415778/1‑70 (MQ=255)

                        aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCg                                                                >  1:126606/1‑50 (MQ=255)

                        aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                                                       >  1:4223055/1‑59 (MQ=255)

                        aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                           >  1:2878984/1‑71 (MQ=255)

                         ggTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                           >  1:3890918/1‑70 (MQ=255)

                          gTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCa                                         <  1:2235105/71‑1 (MQ=255)

                          gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTg                                                             <  1:632709/51‑1 (MQ=255)

                          gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                                                     <  1:261161/59‑1 (MQ=255)

                          gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                                                     <  1:2622058/59‑1 (MQ=255)

                          gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                                                     <  1:4955373/59‑1 (MQ=255)

                          gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                                                     <  1:554302/59‑1 (MQ=255)

                            aCTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCag                                        >  1:5457882/1‑69 (MQ=255)

                             cTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCaga                                       <  1:1211790/70‑3 (MQ=255)

                             cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                      <  1:3691402/71‑1 (MQ=255)

                             cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                       <  1:3172040/70‑1 (MQ=255)

                             cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                       <  1:1549848/70‑1 (MQ=255)

                             cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                       <  1:2417283/70‑1 (MQ=255)

                             cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                       <  1:1625402/70‑1 (MQ=255)

                               gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                         >  1:1108039/1‑66 (MQ=255)

                                aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa                                     <  1:5462102/69‑1 (MQ=255)

                                aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                    <  1:4349864/70‑1 (MQ=255)

                                 tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCAt                                                  >  1:2717923/1‑55 (MQ=255)

                                 tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACgg                                  >  1:4744885/1‑71 (MQ=255)

                                   cGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                <  1:2159368/71‑1 (MQ=255)

                                   cGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                <  1:2631997/71‑1 (MQ=255)

                                     cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                <  1:1483664/69‑1 (MQ=255)

                                       cAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCt                                                      <  1:1420553/45‑1 (MQ=255)

                                         ggTACATACGCAGCTGCCTTGCCGCTGCGTACGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGt                          >  1:1046851/1‑71 (MQ=255)

                                          gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTc                         <  1:1125899/71‑1 (MQ=255)

                                            aCATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    >  1:1756262/1‑42 (MQ=255)

                                            aCATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCa                                                   >  1:1705886/1‑43 (MQ=255)

                                            aCATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGc                       >  1:4283257/1‑71 (MQ=255)

                                             cATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    <  1:600178/41‑1 (MQ=255)

                                             cATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    >  1:802483/1‑41 (MQ=255)

                                               tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    <  1:1082178/39‑1 (MQ=255)

                                               tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                    <  1:396009/39‑1 (MQ=255)

                                               tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGc                                          <  1:2466572/49‑1 (MQ=255)

                                               tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                       <  1:5627644/52‑1 (MQ=255)

                                               tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                       <  1:1697807/52‑1 (MQ=255)

                                                aCGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCaa                   >  1:3306124/1‑71 (MQ=255)

                                                aCGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCaa                   >  1:5118750/1‑71 (MQ=255)

                                                 cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                    <  1:4641982/53‑1 (MQ=255)

                                                 cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                  >  1:2123252/1‑71 (MQ=255)

                                                 cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                  >  1:4691553/1‑71 (MQ=255)

                                                  gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                       <  1:4283010/49‑1 (MQ=255)

                                                   cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                <  1:2527354/71‑1 (MQ=255)

                                                   cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                >  1:1563921/1‑71 (MQ=255)

                                                   cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                <  1:2737867/71‑1 (MQ=255)

                                                    aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc               <  1:447287/71‑1 (MQ=255)

                                                    aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc               <  1:5590133/71‑1 (MQ=255)

                                                      cTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTctg             <  1:5453895/71‑3 (MQ=255)

                                                       tGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGa                                            >  1:3208350/1‑39 (MQ=255)

                                                       tGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGa                                            >  1:3548325/1‑39 (MQ=255)

                                                         ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa          <  1:4999601/71‑1 (MQ=255)

                                                         ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa          <  1:4446137/71‑1 (MQ=255)

                                                         ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCa           >  1:364113/1‑70 (MQ=255)

                                                          gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAAc                                    >  1:3872822/3‑44 (MQ=12)

                                                          gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCag          <  1:729192/68‑2 (MQ=255)

                                                            tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCTAAAACGGCAAATTGTCGCCCAAttt                >  1:2746410/1‑62 (MQ=255)

                                                            tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                >  1:1634008/1‑62 (MQ=255)

                                                            tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                >  1:3171260/1‑62 (MQ=255)

                                                            tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa        >  1:4206615/1‑70 (MQ=255)

                                                            tGCCGCTGCGTCCGGTGATGGTCTTCAT‑CCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                >  1:1560527/1‑61 (MQ=255)

                                                             gCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGt                          >  1:2903263/1‑51 (MQ=255)

                                                             gCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACt      >  1:1614448/1‑71 (MQ=255)

                                                              ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTc     <  1:1531971/71‑1 (MQ=255)

                                                              ccGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTc     >  1:352963/1‑71 (MQ=255)

                                                               cGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTcc    <  1:1330329/71‑2 (MQ=255)

                                                                gCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAtt                            <  1:1895200/46‑1 (MQ=255)

                                                                gCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGc   <  1:774436/71‑1 (MQ=255)

                                                                 cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt  >  1:3564061/1‑71 (MQ=255)

CGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCT  >  minE/2570766‑2570901

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