BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,556	T→C	20.9%	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	20.9%	217.6 / 55.5	110	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (38/49); new base C (11/12); total (49/61)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 114 total aligned reads displayed.

AGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763504‑1763624

aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    <  1:2380503/71‑1 (MQ=21)

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

  tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                  <  1:558882/71‑1 (MQ=34)

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

      gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                 <  1:891315/68‑1 (MQ=35)

      gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                 <  1:145518/68‑1 (MQ=35)

       aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:487884/1‑62 (MQ=33)

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

       aCCGTGTCTCAGTTACAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:1456122/71‑1 (MQ=11)

          gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:2430026/71‑1 (MQ=21)

          gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          >  1:2111815/1‑71 (MQ=21)

          gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          >  1:1972860/1‑71 (MQ=21)

          gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          >  1:1508757/1‑71 (MQ=21)

           tgtCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcc                                         >  1:2255707/1‑71 (MQ=11)

           tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                            >  1:2208221/1‑68 (MQ=25)

             tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:2281983/66‑1 (MQ=35)

             tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcccc                                       <  1:647968/71‑1 (MQ=34)

              ctcAGTTCAAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            >  1:3649597/1‑65 (MQ=25)

               tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                               <  1:3679460/61‑1 (MQ=35)

               tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:758477/64‑1 (MQ=35)

               tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     >  1:2111049/1‑71 (MQ=34)

                cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:784932/55‑1 (MQ=35)

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

                 aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:1779182/71‑1 (MQ=35)

                  gTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                  <  1:427387/71‑1 (MQ=38)

                   ttCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                      <  1:1671640/50‑1 (MQ=11)

                   ttCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                      <  1:1750824/50‑1 (MQ=11)

                   ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 <  1:1056026/71‑1 (MQ=39)

                   ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                <  1:2743661/56‑1 (MQ=255)

                   ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                <  1:558555/56‑1 (MQ=255)

                    tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 <  1:1266359/70‑1 (MQ=40)

                     ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:111776/70‑1 (MQ=255)

                     ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:2044260/1‑69 (MQ=40)

                     ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:3006784/1‑69 (MQ=40)

                      cAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                                   >  1:1947051/1‑50 (MQ=38)

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

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

                         tttgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:533928/42‑1 (MQ=25)

                         tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGAGCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           <  1:1511741/71‑1 (MQ=255)

                          gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     >  1:2755377/1‑44 (MQ=38)

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

                           tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          <  1:326166/70‑1 (MQ=255)

                           tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          <  1:599074/70‑1 (MQ=255)

                           tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACGAACAAGCTaa                          <  1:3525364/70‑1 (MQ=255)

                           tgGCTGGTCATCCT‑TCAGACCAGCTAGGGATCGTCGCCTTgg                                                     >  1:80830/1‑42 (MQ=14)

                             gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:40463/42‑1 (MQ=38)

                                ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     >  1:2497126/1‑54 (MQ=35)

                                ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          <  1:805875/65‑1 (MQ=255)

                                 gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                      >  1:1034269/1‑52 (MQ=35)

                                 gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                >  1:2653829/1‑58 (MQ=255)

                                 gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                    <  1:3674005/70‑1 (MQ=255)

                                 gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   >  1:3004412/1‑71 (MQ=255)

                                 gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   >  1:2158517/1‑71 (MQ=255)

                                  tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:1152383/35‑1 (MQ=37)

                                  tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                                   >  1:2193491/1‑38 (MQ=40)

                                  tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                  <  1:1040736/71‑1 (MQ=255)

                                  tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                  <  1:730055/71‑1 (MQ=255)

                                  tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                  <  1:2836870/71‑1 (MQ=255)

                                  tCATCCTCTCAGACCAGCTAGGGATCGGCGCCTTGGTGAGCCg                                              <  1:1939004/43‑1 (MQ=255)

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

                                  tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACAAGCTAATCCCATCt                  <  1:3011193/71‑1 (MQ=14)

                                   cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           <  1:703833/61‑1 (MQ=255)

                                   cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 >  1:2751326/1‑71 (MQ=255)

                                   cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                  <  1:1007611/70‑1 (MQ=255)

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

                                    aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:793509/1‑71 (MQ=255)

                                    aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:1055405/1‑71 (MQ=255)

                                     tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         <  1:3071162/61‑1 (MQ=255)

                                     tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         >  1:2096065/1‑61 (MQ=255)

                                     tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         <  1:2504072/61‑1 (MQ=255)

                                     tCCTCTCAGACCAGCTAGGGATCCTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:166610/1‑70 (MQ=255)

                                        tctcAGACCAGGTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:3530331/39‑1 (MQ=255)

                                        tctcAGACCAGCTAGGGCTCGTCGCCTTGGTGAGCCGtt                                            <  1:3525807/39‑1 (MQ=255)

                                        tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:1068761/48‑1 (MQ=255)

                                        tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca             >  1:1135996/1‑70 (MQ=255)

                                        tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca             >  1:1628578/1‑70 (MQ=255)

                                         ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:1094660/67‑1 (MQ=255)

                                         ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:1346736/67‑1 (MQ=255)

                                          tcAGACCAGTTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt          >  1:1239224/1‑71 (MQ=255)

                                          tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:1652261/36‑1 (MQ=40)

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

                                          tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt          >  1:481156/1‑71 (MQ=255)

                                          tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt          >  1:1990575/1‑71 (MQ=255)

                                          tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt          >  1:601127/1‑71 (MQ=255)

                                           cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         >  1:1434136/1‑55 (MQ=255)

                                           cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc         >  1:2647512/1‑71 (MQ=255)

                                            aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca           >  1:1275123/1‑68 (MQ=255)

                                            aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        >  1:2865699/1‑71 (MQ=255)

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

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

                                             gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          <  1:358859/52‑1 (MQ=255)

                                             gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          <  1:2677510/52‑1 (MQ=255)

                                              aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa      <  1:2619304/71‑1 (MQ=255)

                                               ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt     <  1:2578281/71‑1 (MQ=255)

                                               ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt     <  1:1740104/71‑1 (MQ=255)

                                                cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg    >  1:2522687/1‑71 (MQ=255)

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

                                                   cTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  >  1:2971014/1‑70 (MQ=255)

                                                    tAGGGATCGTCGCCTTGTTGAGCCGTTACCCCACCAACAAGCTAAt                         >  1:3435490/1‑46 (MQ=255)

                                                    tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:3812435/1‑55 (MQ=255)

AGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763504‑1763624

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