BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,556	T→C	33.6%	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	33.6%	128.6 / 102.1	110	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (29/44); new base C (11/26); total (40/70)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCC  >  minE/1763502‑1763618

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                >  1:2864187/1‑71 (MQ=17)

 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                               <  1:2782081/71‑1 (MQ=21)

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

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

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

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

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

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                                        >  1:2592439/1‑71 (MQ=34)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGGGAGcc                                         <  1:213908/70‑1 (MQ=25)

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

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

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

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCttt                                      <  1:2198610/71‑1 (MQ=11)

               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                  >  1:735155/1‑70 (MQ=35)

               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                  >  1:1334619/1‑70 (MQ=35)

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGCCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcccaa                                >  1:382529/1‑68 (MQ=255)

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

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

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                               >  1:1120202/1‑69 (MQ=35)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                <  1:3303659/68‑1 (MQ=35)

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

                         aGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                       >  1:2855628/1‑71 (MQ=255)

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

                           tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                      <  1:1389696/70‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                         <  1:3034766/51‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                         <  1:1500057/51‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                         <  1:2027074/51‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                         <  1:1693132/51‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                         <  1:3713203/51‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     <  1:32964/71‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     <  1:3435621/71‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     <  1:457165/71‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     <  1:458118/71‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     <  1:968673/71‑1 (MQ=255)

                            gtgGCTGGTCATCCTCTGAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                    >  1:1171459/1‑71 (MQ=255)

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

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

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

                            gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     <  1:1754889/70‑1 (MQ=255)

                            gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                     <  1:2677027/70‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                   <  1:1046117/71‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                   <  1:2788200/71‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                   <  1:3541034/71‑1 (MQ=255)

                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                      <  1:1347038/66‑1 (MQ=255)

                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                      <  1:1785165/66‑1 (MQ=255)

                               gCTGGTCATCCTCT‑AGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                      <  1:2533530/65‑1 (MQ=255)

                                cgggTCATCCTCTCAGACCAGCTAGGGATCGTCTCCTTGGTGAGcc                                         <  1:2112674/44‑1 (MQ=11)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:1307101/1‑50 (MQ=38)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:2164671/1‑50 (MQ=38)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:63771/1‑50 (MQ=38)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:1865097/1‑50 (MQ=38)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATtcc                <  1:1033401/71‑1 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:1393914/71‑1 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:3155937/71‑1 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:2543370/71‑1 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:3300630/71‑1 (MQ=255)

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

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

                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:295219/70‑1 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGGGAGCCGTTACCCCACCAACAAGCTAATccc                <  1:948882/70‑1 (MQ=255)

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

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

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

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

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

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

                                     cATCCTCTCAGACCAGCCAGGGAGCGTCGCCTTGGTGAGc                                          <  1:2034552/40‑1 (MQ=14)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                       <  1:1649689/41‑1 (MQ=38)

                                        ccTCTCAGACCAGCTAGGGATCATCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg         >  1:966940/1‑70 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                    >  1:3239172/1‑58 (MQ=255)

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

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

                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                          >  1:3574496/1‑36 (MQ=255)

                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                     >  1:753533/1‑41 (MQ=255)

                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca       >  1:1932286/1‑71 (MQ=255)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                              <  1:1699565/46‑1 (MQ=38)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                              <  1:3468902/46‑1 (MQ=38)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                              <  1:4090671/46‑1 (MQ=38)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                              <  1:1093736/46‑1 (MQ=38)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                              <  1:1770767/46‑1 (MQ=38)

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

                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca     >  1:3736701/1‑71 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                       >  1:1470048/1‑35 (MQ=40)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:3116627/47‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:432221/47‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:2050892/47‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:2616085/47‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:673031/47‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                          <  1:3511004/47‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                               gACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc  <  1:1828893/70‑1 (MQ=255)

                                                aCCAGCCAGGGATCGTCGCCTTTGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc  <  1:2772757/69‑1 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCC  >  minE/1763502‑1763618

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