BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,558	G→C	6.9%	noncoding (286/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,558	0	G	C	6.9%	460.8 / 15.1	144	noncoding (286/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (75/59); new base C (2/8); total (77/67)
Fisher's exact test for biased strand distribution p-value = 4.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGA  >  minE/1763502‑1763628

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

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

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

   gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                          >  1:1851735/1‑68 (MQ=18)

    tCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGtag                                                      <  1:3191511/71‑3 (MQ=11)

       ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                   >  1:3133810/1‑71 (MQ=34)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           >  1:360905/1‑62 (MQ=31)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           >  1:1824045/1‑62 (MQ=31)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTACGGATCGTCGCCtt                                                           >  1:4334884/1‑62 (MQ=14)

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

         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                 >  1:6318659/1‑71 (MQ=21)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGGGCCttt                                                <  1:3716574/71‑1 (MQ=11)

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

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

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

           cGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCttt                                                <  1:6856590/70‑1 (MQ=11)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                               <  1:5540168/71‑1 (MQ=34)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                               <  1:2758587/71‑1 (MQ=34)

             tgtCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTAcc                                             <  1:1863689/71‑1 (MQ=11)

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

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

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

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

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

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

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

                     ttCCAGTGTGGCTTGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                          <  1:6701931/50‑1 (MQ=11)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                          <  1:1234643/50‑1 (MQ=255)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCACCCAGGGATCGTCGCCTTg                                                          <  1:4796326/50‑1 (MQ=11)

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

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

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

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                        >  1:1002540/1‑50 (MQ=255)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                        >  1:6685291/1‑50 (MQ=255)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                     >  1:4125639/1‑69 (MQ=255)

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

                           tgtgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCATTTACCtcacca                                       <  1:4904356/63‑1 (MQ=11)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:3916158/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:1959860/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:3898897/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:1204011/43‑1 (MQ=255)

                              ggCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                             >  1:5636561/1‑70 (MQ=11)

                               tcTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                          >  1:5732985/2‑40 (MQ=37)

                               gCTTGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                          >  1:3705274/1‑40 (MQ=255)

                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                          >  1:4559303/1‑40 (MQ=37)

                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                          >  1:5830579/1‑40 (MQ=37)

                               gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:5252190/41‑1 (MQ=255)

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

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

                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                       >  1:4630944/1‑57 (MQ=38)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                                      >  1:1832984/1‑42 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                        <  1:4512604/56‑1 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                        <  1:2417011/56‑1 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                     <  1:1425607/42‑1 (MQ=255)

                                      aTCCTCTCAGACCAGCTAGGGATCTTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATACCATc                       >  1:5461882/1‑68 (MQ=255)

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

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

                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                     >  1:4600794/1‑70 (MQ=255)

                                      aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                                 >  1:2678150/1‑58 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                                 >  1:4182043/1‑41 (MQ=38)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   >  1:324328/1‑71 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   >  1:6387257/1‑71 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   >  1:2211759/1‑71 (MQ=255)

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

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

                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTGCCCCACCAACAAGCTAATCCCATCTggg                   >  1:436002/1‑70 (MQ=14)

                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   >  1:3798902/1‑70 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                  <  1:3768868/70‑1 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                  <  1:1995088/70‑1 (MQ=255)

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

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                               >  1:925043/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                               >  1:4780432/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                               >  1:4623624/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                               >  1:4159876/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                               >  1:5363663/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                               >  1:3859846/1‑38 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                          >  1:1572595/1‑42 (MQ=38)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                  >  1:3075915/1‑66 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                  >  1:1005194/1‑66 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                  >  1:3045631/1‑66 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                  >  1:5845895/1‑66 (MQ=255)

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

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

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

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGAACATcc            <  1:899174/71‑1 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc             >  1:1037817/1‑70 (MQ=255)

                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa          >  1:1523031/1‑71 (MQ=255)

                                                      caGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGcta    <  1:3504305/70‑3 (MQ=255)

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

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

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCccccca                                       >  1:6011880/1‑35 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt         <  1:4552700/65‑1 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc      >  1:1847715/1‑68 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc      >  1:2723595/1‑68 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag   <  1:2193679/71‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa    <  1:1051333/69‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa    <  1:3110735/69‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa    <  1:323108/69‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAaga  <  1:68243/71‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag   <  1:5895275/70‑1 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGA  >  minE/1763502‑1763628

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