BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,151	C→T	13.9%	D352D (GAC→GAT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,151	0	C	T	13.9%	315.0 / 34.6	122	D352D (GAC→GAT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (58/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.89e-01

Only 100 of 125 total aligned reads displayed.

ACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATT  >  minE/2571082‑2571221

aCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAcc                                                                       <  1:1789316/71‑1 (MQ=255)

aCTGGACACTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc                                                                        >  1:801584/1‑70 (MQ=255)

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

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

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

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

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

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

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

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

           aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                            >  1:2296564/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)

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

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

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

                   atCGATGCCTTACACGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    >  1:2142560/1‑71 (MQ=255)

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

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

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

                     cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                                  >  1:3917161/1‑71 (MQ=255)

                      gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                                  <  1:1608881/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)

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

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

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

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

                             tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                            <  1:3022781/69‑3 (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)

                                 cGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc                                      >  1:2823678/1‑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:632030/71‑3 (MQ=255)

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

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

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

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

                                        aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:3723872/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)

                                              gCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGAc                         >  1:3320366/1‑71 (MQ=255)

                                               cccAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                               <  1:2781918/64‑1 (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:2218805/62‑1 (MQ=255)

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

                                                  aGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAg                      >  1:456662/1‑70 (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)

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

                                                          gTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGc              <  1:1996898/70‑1 (MQ=255)

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

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

                                                           tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:3454220/39‑1 (MQ=255)

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

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

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

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

                                                            ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                           >  1:579764/1‑39 (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)

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

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

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

                                                               cTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTAt         <  1:3107350/70‑1 (MQ=255)

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

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

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

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

                                                                tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                              <  1:3016000/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)

                                                                   gACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa    <  1:3872986/71‑1 (MQ=255)

                                                                    aCCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATg                           <  1:1630614/47‑1 (MQ=255)

                                                                    aCCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                         >  1:3422157/1‑49 (MQ=255)

                                                                    aCCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                         >  1:1563447/1‑49 (MQ=255)

                                                                     ccATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt  >  1:3805910/1‑71 (MQ=255)

ACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATT  >  minE/2571082‑2571221

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