BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,148	C→T	13.9%	S351S (TCC→TCT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,148	0	C	T	13.9%	316.9 / 33.9	122	S351S (TCC→TCT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (59/46); new base T (10/7); total (69/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 126 total aligned reads displayed.

TTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTT  >  minE/2571078‑2571217

ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc                                                                       >  1:2782638/1‑71 (MQ=255)

     cTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCa                                                                  >  1:2231903/1‑71 (MQ=255)

        gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                                               >  1:1893495/1‑71 (MQ=255)

        gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAtt                                                                >  1:2102634/1‑70 (MQ=255)

        gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAt                                                                 <  1:181142/69‑1 (MQ=255)

        gACCCTAACGATGATATCGATGCCTTAACCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGaa                                                                    >  1:1334145/1‑65 (MQ=255)

          cccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                                               >  1:2810103/1‑69 (MQ=255)

           ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAg                                                            >  1:1721624/1‑71 (MQ=255)

           ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAg                                                            >  1:3908733/1‑71 (MQ=255)

             tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                                           >  1:2467528/1‑70 (MQ=255)

             tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                                           >  1:3405967/1‑70 (MQ=255)

              aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                         >  1:2523851/1‑71 (MQ=255)

              aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                         >  1:3696805/1‑71 (MQ=255)

                cGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCa                                                             >  1:828070/1‑65 (MQ=255)

                 gatgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGt                                                      >  1:1919951/1‑71 (MQ=255)

                  atgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:613427/1‑68 (MQ=255)

                   tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                    <  1:1883885/71‑1 (MQ=255)

                    gatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                    >  1:1149055/1‑70 (MQ=255)

                     atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                                                  >  1:1771985/1‑71 (MQ=255)

                       atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                         <  1:3848879/62‑1 (MQ=255)

                       atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                                                  <  1:1803250/69‑1 (MQ=255)

                        tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                         <  1:504025/61‑1 (MQ=255)

                        tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                <  1:1753831/70‑1 (MQ=255)

                         cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                                                     >  1:3213189/1‑48 (MQ=255)

                         cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATg                                                       >  1:3206653/1‑62 (MQ=255)

                          gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                              <  1:1608881/70‑1 (MQ=255)

                          gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                              <  1:445920/70‑1 (MQ=255)

                          gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                              <  1:2318931/70‑1 (MQ=255)

                           aTGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGcaa                                                         <  1:3331794/58‑3 (MQ=255)

                           aTGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                         <  1:3485451/58‑1 (MQ=255)

                            tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                                                     >  1:194779/1‑45 (MQ=255)

                            tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                                                     >  1:353665/1‑45 (MQ=255)

                               tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgc                                         >  1:2862143/2‑69 (MQ=255)

                               cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                <  1:2898723/63‑1 (MQ=255)

                               cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:707775/71‑1 (MQ=255)

                                ttACCCGACGGCATCTACGCCCACGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        >  1:3533788/1‑70 (MQ=255)

                                 tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        >  1:3348979/1‑69 (MQ=255)

                                  aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                       <  1:2904217/69‑1 (MQ=255)

                                  aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg                                     >  1:21810/1‑71 (MQ=255)

                                   cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGTTGAAAg                                     <  1:508836/70‑1 (MQ=255)

                                   cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg                                     <  1:3333110/70‑1 (MQ=255)

                                     agACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc                                  >  1:3886495/2‑71 (MQ=255)

                                      gACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc                                  >  1:1792418/1‑70 (MQ=255)

                                      gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct                                 >  1:3001520/1‑71 (MQ=255)

                                      gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct                                 >  1:1197432/1‑71 (MQ=255)

                                       aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:2463630/1‑47 (MQ=255)

                                        cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                               <  1:2438973/71‑3 (MQ=255)

                                        cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                               <  1:632030/71‑3 (MQ=255)

                                        cGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctgc                               >  1:652245/1‑71 (MQ=255)

                                          gCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                              <  1:203597/70‑1 (MQ=255)

                                          gCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGCtctgct                              <  1:30981/70‑1 (MQ=255)

                                           cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                              >  1:2786236/1‑69 (MQ=255)

                                           cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                            <  1:1314144/71‑1 (MQ=255)

                                            aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:3723872/58‑1 (MQ=255)

                                            aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:2312843/58‑1 (MQ=255)

                                             tttACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGt                                    >  1:384371/3‑61 (MQ=255)

                                             tttACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGt                                    >  1:1605840/3‑61 (MQ=255)

                                               tACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGAt                                             >  1:1667617/1‑50 (MQ=255)

                                               tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                         >  1:3258548/1‑38 (MQ=255)

                                               tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctg                                >  1:1551595/1‑63 (MQ=255)

                                                 cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGtt                                                     <  1:298228/40‑1 (MQ=255)

                                                 cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                        <  1:1875898/69‑1 (MQ=255)

                                                 cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATg                       <  1:453814/70‑1 (MQ=255)

                                                    ccAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                     >  1:3008480/1‑69 (MQ=255)

                                                    ccAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                     >  1:1196244/1‑69 (MQ=255)

                                                     cAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCaa                                       <  1:1407586/50‑1 (MQ=255)

                                                     cAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAg                  >  1:2751226/1‑71 (MQ=255)

                                                     cAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                           <  1:2138075/62‑1 (MQ=255)

                                                     cAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                           <  1:2218805/62‑1 (MQ=255)

                                                       ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                <  1:3030388/71‑1 (MQ=255)

                                                       ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                >  1:1054204/1‑71 (MQ=255)

                                                       ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCCGTGATGCgg                                         >  1:277263/1‑46 (MQ=255)

                                                          ggCGGTTTACTCTGATCATTTCAGCCATGTCTCgg                                                 >  1:3648951/1‑35 (MQ=25)

                                                          ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                              >  1:54484/1‑38 (MQ=255)

                                                              gTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                <  1:2082187/64‑1 (MQ=255)

                                                              gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                      <  1:2786197/42‑1 (MQ=255)

                                                              gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                      <  1:247820/42‑1 (MQ=255)

                                                               tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt        >  1:3660374/1‑71 (MQ=255)

                                                                ttACTCCGACCATTTCAGCGATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt        >  1:1965689/1‑70 (MQ=255)

                                                                ttACTCCGACCATTTCAGCCATGTTTCGGTGATGc                                           >  1:1230345/1‑35 (MQ=255)

                                                                ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                       >  1:579764/1‑39 (MQ=255)

                                                                ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa             >  1:3287485/1‑65 (MQ=255)

                                                                ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt        >  1:917966/1‑70 (MQ=255)

                                                                ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt        >  1:1093774/1‑70 (MQ=255)

                                                                 tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAg                  <  1:1790166/59‑1 (MQ=255)

                                                                 tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTa      >  1:1522456/1‑71 (MQ=255)

                                                                  aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTAt     >  1:1135165/1‑71 (MQ=255)

                                                                  aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTAt     >  1:2720778/1‑71 (MQ=255)

                                                                  aCTCCGAACATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg           >  1:1763574/1‑65 (MQ=255)

                                                                   cTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                     >  1:397751/1‑54 (MQ=255)

                                                                   cTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCac                  <  1:1950544/57‑2 (MQ=255)

                                                                   cTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa             >  1:216338/1‑62 (MQ=255)

                                                                   cTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATc    <  1:2534019/71‑1 (MQ=255)

                                                                    tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                         >  1:3063680/1‑49 (MQ=255)

                                                                    tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                          <  1:1765018/48‑1 (MQ=255)

                                                                    tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                          <  1:3016000/48‑1 (MQ=255)

                                                                    tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                          <  1:3938729/48‑1 (MQ=255)

                                                                    tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                          <  1:3054275/48‑1 (MQ=255)

                                                                      cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct                                 >  1:287753/1‑39 (MQ=255)

                                                                      cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa             >  1:911869/1‑59 (MQ=255)

TTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTT  >  minE/2571078‑2571217

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