BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,556	T→C	24.8%	noncoding (288/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,763,556	0	T	C	24.8%	230.1 / 90.5	134	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (44/56); new base C (20/13); total (65/69)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 135 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763502‑1763624

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       >  1:3666093/1‑70 (MQ=31)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:1038247/70‑1 (MQ=31)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:459250/70‑1 (MQ=31)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:2600822/70‑1 (MQ=31)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:3267333/71‑1 (MQ=32)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATACTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:1387041/1‑71 (MQ=17)

  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:2304749/71‑1 (MQ=34)

  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:787008/71‑1 (MQ=34)

   gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    >  1:1851696/1‑70 (MQ=35)

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:2911336/66‑1 (MQ=31)

     cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                  <  1:1552146/70‑1 (MQ=35)

     cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                 >  1:3440050/1‑71 (MQ=21)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     >  1:1838709/1‑66 (MQ=35)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     >  1:89084/1‑66 (MQ=35)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                >  1:1064937/1‑71 (MQ=34)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                                  >  1:1974647/1‑69 (MQ=25)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                >  1:1872987/1‑71 (MQ=21)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                               <  1:892409/70‑1 (MQ=35)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                <  1:1983882/69‑1 (MQ=35)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                <  1:2349530/69‑1 (MQ=35)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                <  1:2387772/69‑1 (MQ=25)

         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:1111757/71‑1 (MQ=34)

         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:1321794/71‑1 (MQ=34)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:1681976/1‑70 (MQ=35)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:3459014/62‑1 (MQ=35)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                               >  1:2943825/1‑67 (MQ=25)

            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:3483634/71‑1 (MQ=34)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTggc                                                    >  1:2113468/1‑59 (MQ=25)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          >  1:2886690/1‑70 (MQ=35)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          >  1:1914604/1‑70 (MQ=25)

              gtCTCAGTTCCAGTGCGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:1699843/59‑1 (MQ=255)

                ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                      >  1:638504/1‑71 (MQ=21)

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                <  1:2033865/60‑1 (MQ=35)

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     >  1:1886031/1‑71 (MQ=21)

                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     >  1:2773888/1‑70 (MQ=35)

                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     >  1:1492610/1‑70 (MQ=25)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:2208610/1‑70 (MQ=35)

                    gTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:2433348/50‑1 (MQ=32)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:2616194/50‑1 (MQ=35)

                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:1398573/1‑70 (MQ=40)

                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:690178/1‑70 (MQ=40)

                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:3080610/1‑70 (MQ=40)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                  <  1:326245/52‑1 (MQ=35)

                        cAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                               <  1:196842/70‑1 (MQ=255)

                         aGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:3247073/45‑1 (MQ=32)

                          gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:3125497/71‑1 (MQ=255)

                          gtgtGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                      >  1:860322/1‑61 (MQ=255)

                          gtgtGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                      >  1:1572820/1‑61 (MQ=255)

                          gtgtGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:3432653/71‑1 (MQ=255)

                          gtgtGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCcca                                      <  1:431112/61‑1 (MQ=11)

                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          >  1:1130007/1‑71 (MQ=255)

                            gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          >  1:738911/1‑71 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     >  1:102039/1‑43 (MQ=38)

                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     >  1:1150785/1‑43 (MQ=38)

                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         >  1:1999865/1‑71 (MQ=255)

                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                       >  1:2768415/1‑71 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     <  1:1596719/40‑1 (MQ=40)

                                cTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                      <  1:255662/71‑1 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:42057/1‑71 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                <  1:617659/43‑1 (MQ=38)

                                  ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                    >  1:3906406/1‑71 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                    >  1:499085/1‑71 (MQ=255)

                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:3587879/69‑1 (MQ=255)

                                    tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    <  1:3166487/53‑1 (MQ=35)

                                    tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                  <  1:3773239/71‑1 (MQ=255)

                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 >  1:3582361/1‑71 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 <  1:188649/69‑1 (MQ=255)

                                       tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:2813571/58‑1 (MQ=255)

                                        ccTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:1012564/1‑44 (MQ=38)

                                        ccTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATACCATCTggg               >  1:1365011/1‑70 (MQ=255)

                                        ccTCTCAGAACAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                      >  1:932708/1‑47 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca             >  1:2261621/1‑71 (MQ=255)

                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          >  1:932486/1‑42 (MQ=255)

                                         ctctCAGACCAGCAAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:3767508/1‑41 (MQ=255)

                                          tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:2705588/51‑1 (MQ=255)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                        >  1:1300953/1‑42 (MQ=38)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca           >  1:2094946/1‑71 (MQ=255)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca           >  1:3405011/1‑71 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGCGCCGTTACCCCACCAACAAGCTAATCCCATc                   >  1:3987749/1‑62 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          >  1:3725527/1‑39 (MQ=40)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:3746897/60‑1 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 <  1:1336465/64‑1 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac            <  1:2837129/69‑1 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt          >  1:437235/1‑71 (MQ=255)

                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         <  1:3498215/55‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:3073409/71‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc         <  1:3021600/70‑1 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          >  1:2164523/1‑53 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:4132595/71‑1 (MQ=255)

                                               gACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        >  1:2264567/1‑70 (MQ=255)

                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   >  1:1759390/1‑57 (MQ=255)

                                                 ccAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt     <  1:1316207/71‑1 (MQ=255)

                                                 ccAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt     >  1:802607/1‑71 (MQ=255)

                                                  cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                              >  1:245132/1‑45 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  <  1:617076/71‑1 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  <  1:3755854/71‑1 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  <  1:2967163/71‑1 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  <  1:3791551/71‑1 (MQ=255)

                                                     cTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:561178/64‑1 (MQ=255)

                                                     cTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAtaga  >  1:667188/1‑67 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763502‑1763624

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