BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,825	T→C	8.6%	L244L (TTG→CTG)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,825	0	T	C	8.6%	287.6 / 13.4	104	L244L (TTG→CTG)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (40/55); new base C (6/3); total (46/58)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

Only 100 of 105 total aligned reads displayed.

TCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAA  >  minE/2570757‑2570893

tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCtt                                                                     <  1:124460/70‑1 (MQ=255)

  aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                                                                  >  1:2726988/1‑71 (MQ=255)

   ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:96965/71‑1 (MQ=255)

   ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:3065483/71‑1 (MQ=255)

   ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:3303062/71‑1 (MQ=255)

   ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:3667828/71‑1 (MQ=255)

   ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:4652072/71‑1 (MQ=255)

   ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:5155473/71‑1 (MQ=255)

   ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:1277911/71‑1 (MQ=255)

   ttCGCGCGCCGAGCAAAGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 <  1:1593534/71‑1 (MQ=255)

    tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                                                >  1:4405786/1‑71 (MQ=255)

     cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                                                 >  1:4819656/1‑69 (MQ=255)

         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:2791965/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

                    ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                                        >  1:4733665/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)

                          gTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                                                                  <  1:4329799/47‑1 (MQ=255)

                           ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                                                <  1:3557122/48‑1 (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)

                                 aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                                                >  1:4093893/1‑42 (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:554302/59‑1 (MQ=255)

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

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

                                   gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                                             <  1:261161/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:1625402/70‑1 (MQ=255)

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

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

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

                                        gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                       >  1:1758223/1‑60 (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)

                                                   gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                              <  1:2733484/58‑1 (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:802483/1‑41 (MQ=255)

                                                      cATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                            <  1:600178/41‑1 (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:1697807/52‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                            cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt        <  1:2527354/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:4446137/71‑1 (MQ=255)

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

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

TCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAA  >  minE/2570757‑2570893

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