BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,559	G→T	14.0%	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	14.0%	255.4 / 23.1	100	noncoding (285/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (34/52); new base T (2/12); total (36/64)
Fisher's exact test for biased strand distribution p-value = 7.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

Only 100 of 108 total aligned reads displayed.

AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAG  >  Exported/1763501‑1763629

aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                            >  1:423641/1‑71 (MQ=14)

aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                            >  1:3045702/1‑71 (MQ=14)

aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                            >  1:2080024/1‑71 (MQ=14)

aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                            >  1:449594/1‑71 (MQ=14)

aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                            >  1:2701320/1‑71 (MQ=14)

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

 ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                           <  1:2274331/71‑1 (MQ=17)

  gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTgg                                                          <  1:490577/71‑1 (MQ=11)

  gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTgg                                                          <  1:2481990/71‑1 (MQ=11)

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

    gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTg                                                        <  1:566110/71‑1 (MQ=11)

    gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                         >  1:3175186/1‑70 (MQ=25)

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

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

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

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

       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                     >  1:468084/1‑71 (MQ=21)

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

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

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

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

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

         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                                                   >  1:268553/1‑71 (MQ=21)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:2343415/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:344131/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:1759088/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:1094390/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:2628856/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:1755242/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:158616/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:1250629/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:3369470/71‑2 (MQ=11)

         gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                                                   <  1:593331/71‑2 (MQ=11)

          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                                        >  1:3395954/1‑65 (MQ=35)

          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                                        >  1:1622022/1‑65 (MQ=35)

          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                                        >  1:2592988/1‑65 (MQ=35)

          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                  <  1:2657554/71‑1 (MQ=21)

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

          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                  <  1:2984759/71‑1 (MQ=21)

           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                 <  1:2269088/71‑1 (MQ=34)

           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                 <  1:1534099/71‑1 (MQ=34)

            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                 <  1:1459344/70‑1 (MQ=35)

            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                 >  1:2770499/1‑70 (MQ=35)

            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                                >  1:1172289/1‑71 (MQ=34)

                tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGcc                                                    >  1:3054889/1‑63 (MQ=11)

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

                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                                 <  1:2906679/62‑1 (MQ=25)

                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                                >  1:2429005/1‑63 (MQ=25)

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

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

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

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

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                <  1:3391515/71‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                <  1:203697/71‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                <  1:123096/71‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                <  1:2716738/71‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACAAGCTa                                <  1:1654433/71‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                           <  1:1004419/44‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                           <  1:2456802/44‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                           <  1:2172058/44‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                           <  1:1141152/44‑1 (MQ=255)

                              tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                               <  1:2333982/70‑1 (MQ=255)

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

                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                               <  1:3071650/52‑1 (MQ=255)

                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                               <  1:1472666/52‑1 (MQ=255)

                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                             <  1:2296588/70‑1 (MQ=255)

                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                             <  1:736120/70‑1 (MQ=255)

                                  tGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                              >  1:42192/1‑67 (MQ=11)

                                    atCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATccga                          <  1:1856091/68‑3 (MQ=11)

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

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

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

                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                     <  1:3165672/70‑1 (MQ=255)

                                           tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca                  >  1:257865/1‑70 (MQ=255)

                                            ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                                  >  1:2983573/1‑53 (MQ=255)

                                             tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                              >  1:1039970/1‑40 (MQ=40)

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

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

                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt               <  1:3278459/69‑1 (MQ=255)

                                               aGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTaa                               <  1:662917/53‑1 (MQ=11)

                                               aGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTaa                               <  1:2262613/53‑1 (MQ=11)

                                               aGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTaa                               <  1:2266679/53‑1 (MQ=11)

                                                 aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa           <  1:2777904/71‑1 (MQ=255)

                                                  ccAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt          >  1:298506/1‑71 (MQ=255)

                                                    aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                        <  1:2272994/39‑1 (MQ=25)

                                                    aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                        <  1:1166168/39‑1 (MQ=25)

                                                    aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                        <  1:1545529/39‑1 (MQ=25)

                                                    aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                        <  1:2374791/39‑1 (MQ=25)

                                                    aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                        <  1:3302956/39‑1 (MQ=25)

                                                    aGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg        <  1:1970495/71‑1 (MQ=255)

                                                     gccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc       <  1:2416041/68‑1 (MQ=255)

                                                     gccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc       <  1:873395/68‑1 (MQ=255)

                                                      cTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCa      <  1:2564884/71‑1 (MQ=255)

                                                        acggATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCt                                 >  1:1614449/3‑42 (MQ=14)

                                                         gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                 <  1:1278351/41‑1 (MQ=255)

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

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

                                                          ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             <  1:2549819/60‑1 (MQ=255)

                                                          ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc             <  1:2654497/60‑1 (MQ=255)

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

AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAG  >  Exported/1763501‑1763629

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