BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,490	C→A	5.8%	P89T (CCG→ACG)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,490	0	C	A	5.8%	484.7 / 8.6	138	P89T (CCG→ACG)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand): ref base C (50/80); new base A (3/5); total (53/85)
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.86e-01

Only 100 of 146 total aligned reads displayed.

GGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTC  >  minE/2086420‑2086555

 gcgcACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAccc                                                                  <  1:3071467/71‑1 (MQ=255)

 gcgcACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAccc                                                                  <  1:4715820/71‑1 (MQ=255)

    caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                               <  1:1096466/71‑1 (MQ=255)

      ccgcacGCTGATCATGCTCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCg                                                                 <  1:1838912/67‑1 (MQ=255)

      ccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCg                                                                 <  1:5528447/67‑1 (MQ=255)

      ccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCC‑TAGTGAACCCGCTGAt                                                            <  1:1682535/71‑1 (MQ=255)

       cgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGa                                                             >  1:3249511/1‑70 (MQ=255)

       cgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                            <  1:2694658/71‑1 (MQ=255)

       cgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                            <  1:3177949/71‑1 (MQ=255)

       cgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                            <  1:3518265/71‑1 (MQ=255)

       cgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                            >  1:4833595/1‑71 (MQ=255)

        gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                               >  1:4154122/1‑67 (MQ=255)

        gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATg                                                           <  1:2297165/71‑1 (MQ=255)

        gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATg                                                           <  1:6734877/71‑1 (MQ=255)

           cGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTg                                                        >  1:5029373/1‑71 (MQ=255)

           cGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTg                                                        >  1:1028991/1‑71 (MQ=255)

             cTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAcc                                                                   >  1:5049244/1‑58 (MQ=255)

              tGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTgg                                                       >  1:2050585/1‑69 (MQ=255)

               gATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaa                                                     >  1:1503769/1‑70 (MQ=255)

               gATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGaa                                                     >  1:397050/1‑70 (MQ=255)

                 tCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAAc                                                  >  1:5137224/1‑71 (MQ=255)

                 tCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCAGTAGTGAACCCGCTGATGCTGGAAAAc                                                  <  1:5066057/71‑1 (MQ=255)

                  cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCg                                                                 >  1:5257567/1‑55 (MQ=255)

                  cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCg                                                                 >  1:6107956/1‑55 (MQ=255)

                  cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCg                                                                 >  1:934001/1‑55 (MQ=255)

                  cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGa                                                             >  1:1090041/1‑59 (MQ=255)

                  cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                 >  1:5363982/1‑71 (MQ=255)

                  cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                 >  1:1831127/1‑71 (MQ=255)

                   aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                <  1:303117/71‑1 (MQ=255)

                         aaaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGAc                                          <  1:1429657/71‑1 (MQ=255)

                         aaaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGAc                                          <  1:6956836/71‑1 (MQ=255)

                         aaaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGAc                                          <  1:5718003/71‑1 (MQ=255)

                            cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                               <  1:4691789/47‑1 (MQ=255)

                            cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCtt                                       >  1:4250404/1‑71 (MQ=255)

                            cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCCCt                                                               <  1:4649805/47‑1 (MQ=255)

                               aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                      <  1:4965423/69‑1 (MQ=255)

                               aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                      <  1:4053271/69‑1 (MQ=255)

                               aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                      <  1:1912131/69‑1 (MQ=255)

                               aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                    >  1:2070487/1‑71 (MQ=255)

                               aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                    >  1:5218573/1‑71 (MQ=255)

                               aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                    >  1:1305843/1‑71 (MQ=255)

                                  ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:385220/69‑1 (MQ=255)

                                  ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:3411139/69‑1 (MQ=255)

                                  ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:6233434/69‑1 (MQ=255)

                                  ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:4414076/69‑1 (MQ=255)

                                    aCTGGAACCGTTTATTTCCGATTCCGTAGTGAACCCGCTGATGCTGGaa                                                     <  1:2968154/49‑1 (MQ=255)

                                      tGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaa                                                     <  1:4351229/47‑1 (MQ=255)

                                       ggAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGtt                            >  1:4774984/1‑71 (MQ=255)

                                        gAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGttt                           <  1:1362979/71‑1 (MQ=255)

                                         aaCCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAAc                                                  <  1:6062723/47‑1 (MQ=255)

                                          aCCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                         <  1:4378804/71‑1 (MQ=255)

                                          aCCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                         <  1:6124281/71‑1 (MQ=255)

                                            cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGt                                               <  1:796500/47‑1 (MQ=255)

                                             gTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGtttt                          <  1:4746955/67‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCTCTGATGCTGGAAAACGGTTGGACCtt                                       <  1:3042476/53‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                      <  1:572086/54‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCtt                                       <  1:6040275/53‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:5185889/57‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:2317299/57‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:777389/57‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:1690200/57‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:6582305/57‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:5466217/57‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                   <  1:59024/57‑1 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                      >  1:4023257/1‑70 (MQ=255)

                                              tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTAGACCttt                                      <  1:1403822/54‑1 (MQ=255)

                                               ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                <  1:1169125/43‑1 (MQ=255)

                                               ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                <  1:147592/43‑1 (MQ=255)

                                               ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                <  1:4626173/43‑1 (MQ=255)

                                               ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                <  1:1825607/43‑1 (MQ=255)

                                               ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                      >  1:1579984/1‑53 (MQ=255)

                                               ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                    >  1:1913538/1‑55 (MQ=255)

                                                tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                <  1:1991912/42‑1 (MQ=255)

                                                tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                   >  1:509383/1‑71 (MQ=255)

                                                 aTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGttt                           <  1:1397749/62‑1 (MQ=255)

                                                  tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                  >  1:2227868/1‑70 (MQ=255)

                                                  tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                  >  1:4345738/1‑70 (MQ=255)

                                                  tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                  <  1:2410748/70‑1 (MQ=255)

                                                  tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                  >  1:2000966/1‑70 (MQ=255)

                                                    tccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                <  1:6650770/70‑1 (MQ=255)

                                                    tccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                <  1:6428395/70‑1 (MQ=255)

                                                    tccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAAc               >  1:4988218/1‑71 (MQ=255)

                                                     ccgtttgcgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCt                                        >  1:885386/1‑45 (MQ=255)

                                                     ccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCt                                        >  1:3414909/1‑45 (MQ=255)

                                                      cgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg             <  1:4208835/71‑1 (MQ=255)

                                                      cgtttccgtAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg             <  1:1761916/71‑1 (MQ=255)

                                                       gtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCgg            <  1:3152398/71‑1 (MQ=255)

                                                       gtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg             >  1:3549591/1‑70 (MQ=255)

                                                       gtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg             >  1:213226/1‑70 (MQ=255)

                                                        tttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga           <  1:6431348/71‑1 (MQ=255)

                                                        tttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga           <  1:4167101/71‑1 (MQ=255)

                                                        tttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga           <  1:1632653/71‑1 (MQ=255)

                                                          tccgtAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgag          <  1:3763530/70‑1 (MQ=255)

                                                              tAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACg     >  1:2414992/1‑71 (MQ=255)

                                                                 tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                        <  1:1175807/49‑1 (MQ=255)

                                                                 tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGctc  <  1:5735486/71‑1 (MQ=255)

                                                                   aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgag          >  1:2696594/1‑61 (MQ=255)

                                                                   aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgag          >  1:53997/1‑61 (MQ=255)

                                                                   aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgag          >  1:4706061/1‑61 (MQ=255)

                                                                   aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgag          >  1:2116731/1‑61 (MQ=255)

GGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTC  >  minE/2086420‑2086555

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