BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,556	T→C	61.2%	noncoding (288/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,556	0	T	C	61.2%	69.7 / 115.1	103	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (23/17); new base C (34/29); total (57/46)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 103 total aligned reads displayed.

AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  Exported/1763501‑1763624

aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCATACCAGCTACTGATCGTCGCCtt                                                       >  1:2489581/1‑71 (MQ=11)

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

 ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                       <  1:2856721/70‑1 (MQ=14)

 ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                       <  1:1513451/70‑1 (MQ=14)

 ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                       <  1:3200670/70‑1 (MQ=14)

 ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                       <  1:1704874/70‑1 (MQ=14)

 ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                       <  1:3541036/70‑1 (MQ=14)

  gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     <  1:1069199/71‑1 (MQ=34)

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

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

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

    gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                       >  1:635793/1‑67 (MQ=14)

      cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATGGTCGCCTTGGTGAg                                                 >  1:2489957/1‑71 (MQ=21)

      cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCTTTgg                                                     >  1:2075699/1‑67 (MQ=255)

      cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     >  1:2812842/1‑67 (MQ=25)

         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                                   <  1:2885743/66‑1 (MQ=35)

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

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

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

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

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

            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:2131255/1‑71 (MQ=21)

              tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:319754/1‑71 (MQ=34)

              tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:326635/1‑71 (MQ=21)

                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                        >  1:718334/1‑70 (MQ=25)

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

                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:2347824/71‑1 (MQ=25)

                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:3160557/1‑70 (MQ=25)

                      ttCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACcaac                                 <  1:1432169/71‑1 (MQ=14)

                      ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                  >  1:139546/1‑70 (MQ=255)

                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:2856661/68‑1 (MQ=255)

                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:325792/68‑1 (MQ=255)

                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:2099635/68‑1 (MQ=255)

                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:3065684/1‑70 (MQ=255)

                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:2779802/1‑70 (MQ=255)

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

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

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

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

                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:1796663/1‑68 (MQ=255)

                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                              >  1:2968764/1‑71 (MQ=255)

                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                              >  1:2933418/1‑71 (MQ=255)

                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                              >  1:1508276/1‑71 (MQ=255)

                          aGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                             <  1:2986740/71‑1 (MQ=255)

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

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:20227/1‑55 (MQ=35)

                             gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:2404789/1‑52 (MQ=255)

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

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

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

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

                              tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         <  1:3462107/71‑1 (MQ=255)

                               ggCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                        >  1:662784/1‑71 (MQ=255)

                               ggCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                        >  1:2204479/1‑71 (MQ=255)

                                 cTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:2556846/58‑1 (MQ=255)

                                 cTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:3311226/58‑1 (MQ=255)

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

                                  tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                              >  1:825626/1‑62 (MQ=255)

                                  tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                              >  1:2252877/1‑62 (MQ=255)

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

                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:1801142/45‑1 (MQ=255)

                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:3186557/45‑1 (MQ=255)

                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:3190860/45‑1 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                               >  1:1208789/1‑40 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                               >  1:3442614/1‑40 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                               >  1:3188040/1‑40 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                               >  1:3120799/1‑40 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:1920564/42‑1 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                <  1:2151476/71‑1 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:3656878/1‑71 (MQ=255)

                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                <  1:2444683/71‑1 (MQ=255)

                                        tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:2036044/71‑1 (MQ=255)

                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:2682/71‑1 (MQ=255)

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

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

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

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

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

                                          ctctCAGACAAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca             >  1:2449872/1‑71 (MQ=255)

                                           tctcAGACCAGCCAGGGATCGTCGCCTTGGTTAGCCGTTAc                                          <  1:3094605/41‑1 (MQ=14)

                                           tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:577151/41‑1 (MQ=255)

                                           tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:501423/41‑1 (MQ=255)

                                           tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:2036842/1‑47 (MQ=255)

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

                                             tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca           >  1:561373/1‑70 (MQ=255)

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

                                              cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:3293774/1‑37 (MQ=255)

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

                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                        >  1:1646973/1‑39 (MQ=255)

                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                        >  1:3366957/1‑39 (MQ=255)

                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:508183/1‑58 (MQ=255)

                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:764650/1‑58 (MQ=255)

                                                gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:34246/46‑1 (MQ=255)

                                                gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:331541/46‑1 (MQ=255)

                                                gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg       <  1:2750547/71‑1 (MQ=255)

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

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

                                                   cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:2947195/67‑1 (MQ=255)

                                                   cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg    <  1:2723869/71‑1 (MQ=255)

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

AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  Exported/1763501‑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: `‑`