BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,559	G→T	18.1%	noncoding (285/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,763,559	0	G	T	18.1%	387.1 / 46.7	138	noncoding (285/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (50/63); new base T (7/18); total (57/81)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 143 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAG  >  Exported/1763502‑1763629

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                           >  1:3700193/1‑71 (MQ=32)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCATGGATCGTCGCCtt                                                            <  1:3224739/70‑1 (MQ=255)

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGa                                                       <  1:4498578/71‑1 (MQ=11)

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGa                                                       <  1:426074/71‑1 (MQ=11)

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                           <  1:2004363/67‑1 (MQ=18)

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                         >  1:2543939/1‑69 (MQ=25)

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                                       <  1:3281157/71‑1 (MQ=21)

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                                       <  1:3699791/71‑1 (MQ=21)

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

     cTGGACAGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                      >  1:663053/1‑71 (MQ=255)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                          >  1:4960723/1‑66 (MQ=25)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                          >  1:2143655/1‑66 (MQ=25)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                          >  1:483060/1‑66 (MQ=25)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                          >  1:2142630/1‑66 (MQ=25)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                          >  1:1349280/1‑66 (MQ=25)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                           >  1:742252/1‑65 (MQ=18)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                                    <  1:4658648/70‑1 (MQ=25)

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

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

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                              <  1:5154498/71‑1 (MQ=34)

                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGc                                                     <  1:207018/60‑1 (MQ=11)

                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGc                                                     <  1:4878817/60‑1 (MQ=11)

                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGc                                                     <  1:891750/60‑1 (MQ=11)

                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCCTCGCCTTGGTGAGc                                                     <  1:4588196/60‑1 (MQ=11)

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

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

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

                  cAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcac                                          >  1:1022699/1‑70 (MQ=11)

                  cAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtca                                           >  1:3746126/1‑69 (MQ=11)

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

                      tCCAGTGTGGCTGGTCATCCTCTAAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                      >  1:566425/1‑70 (MQ=255)

                          gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                                  >  1:2036092/1‑70 (MQ=255)

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

                             tgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                              >  1:4628818/1‑71 (MQ=11)

                             tgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                              >  1:1421740/1‑71 (MQ=11)

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

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

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

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                <  1:1336918/51‑1 (MQ=11)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                <  1:212854/51‑1 (MQ=11)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                <  1:3657948/51‑1 (MQ=11)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                <  1:4992091/51‑1 (MQ=11)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                <  1:2505864/51‑1 (MQ=11)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                <  1:4379064/51‑1 (MQ=11)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                <  1:1834253/51‑1 (MQ=11)

                                cTGGTCATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                            >  1:4443353/1‑70 (MQ=255)

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

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

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

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

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

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

                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                        <  1:1085530/71‑1 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATc                        <  1:3401115/67‑1 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                        <  1:1050380/67‑1 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                        <  1:4180275/67‑1 (MQ=255)

                                        ccTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                        <  1:3524311/66‑1 (MQ=255)

                                          tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                 <  1:2355233/55‑1 (MQ=255)

                                          tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                 <  1:1320844/55‑1 (MQ=255)

                                          tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                 <  1:1479530/55‑1 (MQ=255)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                         >  1:4973576/1‑62 (MQ=255)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                         >  1:3195243/1‑62 (MQ=255)

                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                         <  1:2857380/62‑1 (MQ=255)

                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                         <  1:3149788/62‑1 (MQ=255)

                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                         <  1:1723317/62‑1 (MQ=255)

                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                         <  1:4049123/62‑1 (MQ=255)

                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                         <  1:3928392/62‑1 (MQ=255)

                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac                 >  1:4233531/1‑70 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                       >  1:133872/1‑47 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                       >  1:1800905/1‑47 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                       >  1:375083/1‑47 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                       >  1:4438580/1‑47 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                       >  1:2826484/1‑47 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACaaa                                       >  1:4117983/1‑47 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca                  >  1:3872538/1‑68 (MQ=255)

                                             cAGACCATCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc              >  1:4784450/1‑71 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc              >  1:1234723/1‑71 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc              >  1:3560832/1‑71 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc              >  1:4550254/1‑71 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc              >  1:4398980/1‑71 (MQ=255)

                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                       >  1:3015142/1‑62 (MQ=255)

                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                       >  1:836957/1‑62 (MQ=255)

                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                       >  1:1866600/1‑62 (MQ=255)

                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGACCCGTTACCCCACCAACAAGCAAATCCCATCt                       >  1:4051324/1‑62 (MQ=17)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                <  1:4866292/68‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                <  1:190786/68‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                <  1:2888818/68‑1 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                <  1:4346734/68‑1 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             >  1:3721379/1‑71 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             >  1:4733271/1‑71 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             >  1:336297/1‑71 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             <  1:1960879/65‑1 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             <  1:2737834/65‑1 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             <  1:1363365/65‑1 (MQ=255)

                                                      caGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa     <  1:4147964/70‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg            <  1:1744350/62‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg            <  1:4387291/62‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg            <  1:3173839/62‑1 (MQ=255)

                                                         ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAagag  >  1:5191429/1‑71 (MQ=255)

                                                         ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAagag  >  1:3070544/1‑71 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAG  >  Exported/1763502‑1763629

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