BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	41.5%	noncoding (1228/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,359,867	0	G	A	41.5%	116.9 / 160.9	141	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (41/36); new base A (33/25); total (78/63)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 141 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/2359815‑2359934

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   <  1:3987377/71‑1 (MQ=9)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   <  1:564710/70‑1 (MQ=31)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                  <  1:1355553/71‑1 (MQ=255)

  cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   >  1:2830976/1‑69 (MQ=31)

   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                <  1:653960/71‑1 (MQ=34)

   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAg                                                >  1:3522401/1‑71 (MQ=11)

   gCTTAATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcagg                                                 >  1:608013/1‑68 (MQ=9)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:3685340/70‑1 (MQ=35)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             >  1:929185/1‑70 (MQ=35)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:3176425/70‑1 (MQ=35)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                            <  1:952365/70‑1 (MQ=255)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           >  1:17193/1‑70 (MQ=35)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           >  1:2470281/1‑70 (MQ=35)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           <  1:2036201/70‑1 (MQ=255)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           <  1:209031/70‑1 (MQ=255)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:2644620/71‑1 (MQ=34)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   <  1:4585194/60‑1 (MQ=11)

           gTTACTTATGCCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:2962223/71‑1 (MQ=21)

             tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:1269440/1‑71 (MQ=255)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                      <  1:783686/70‑1 (MQ=35)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                 >  1:1662393/1‑59 (MQ=11)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:2072561/1‑70 (MQ=255)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                     >  1:465139/1‑71 (MQ=255)

                 tATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                            <  1:3342173/61‑1 (MQ=35)

                 tATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                            <  1:1140409/61‑1 (MQ=35)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                   <  1:215768/69‑1 (MQ=35)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                 >  1:3762593/1‑71 (MQ=34)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                 >  1:3104807/1‑71 (MQ=34)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                  <  1:3483555/70‑1 (MQ=255)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 >  1:2729010/1‑71 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                <  1:1973318/71‑1 (MQ=34)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                 <  1:274264/70‑1 (MQ=35)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 <  1:54783/70‑1 (MQ=255)

                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   >  1:423850/1‑50 (MQ=11)

                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           >  1:2022405/1‑58 (MQ=11)

                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                              <  1:2723060/71‑1 (MQ=255)

                     tCAGCACTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                  <  1:4392679/51‑1 (MQ=11)

                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             >  1:2879683/1‑52 (MQ=35)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                       <  1:2454631/56‑1 (MQ=35)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                  >  1:122100/1‑45 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                        <  1:1260279/55‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                        <  1:4136902/55‑1 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                    >  1:2282957/1‑58 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                    >  1:3054334/1‑58 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                    >  1:2054159/1‑58 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                             <  1:4017568/65‑1 (MQ=255)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:721464/1‑67 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                       >  1:3831472/1‑70 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                       >  1:3613028/1‑70 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                       >  1:427031/1‑70 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      <  1:3773015/55‑1 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      <  1:3581965/55‑1 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                       >  1:702745/1‑70 (MQ=255)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                       >  1:4319380/1‑70 (MQ=255)

                             agCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                        >  1:932603/2‑69 (MQ=25)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                         <  1:3279864/51‑1 (MQ=11)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:1073557/1‑54 (MQ=11)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      >  1:3601662/1‑70 (MQ=255)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      <  1:4186127/70‑1 (MQ=255)

                               cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           >  1:3952707/1‑48 (MQ=11)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   >  1:2807810/1‑39 (MQ=33)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   >  1:4637892/1‑39 (MQ=33)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           <  1:969170/47‑1 (MQ=38)

                                aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:4207532/1‑52 (MQ=11)

                                aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                       >  1:2917911/1‑51 (MQ=11)

                                aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                   >  1:1602628/1‑71 (MQ=255)

                                aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                   >  1:2326817/1‑71 (MQ=255)

                                aCTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      <  1:2089205/68‑1 (MQ=11)

                                aCTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      <  1:2326521/68‑1 (MQ=11)

                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                            <  1:3600746/61‑1 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:4268354/43‑1 (MQ=38)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           <  1:3715103/45‑1 (MQ=38)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           <  1:1444350/45‑1 (MQ=38)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                         >  1:1420254/1‑47 (MQ=38)

                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:4068785/47‑1 (MQ=38)

                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:2704013/1‑61 (MQ=35)

                                   tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 >  1:1611008/1‑54 (MQ=11)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   >  1:4519058/1‑35 (MQ=33)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:3557654/46‑1 (MQ=38)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                      >  1:1306436/1‑48 (MQ=38)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                <  1:3653045/70‑1 (MQ=40)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                >  1:3412105/1‑70 (MQ=40)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                <  1:2077903/70‑1 (MQ=40)

                                      gATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                  >  1:842277/1‑66 (MQ=11)

                                      gATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                  >  1:4176682/1‑66 (MQ=11)

                                        tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg           <  1:3799113/71‑1 (MQ=255)

                                        tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg           <  1:1372322/71‑1 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc         <  1:3965086/71‑1 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc         >  1:3164847/1‑71 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCCTACAGAACGCTCCCCTACCCAACAACACATAGTGTc         <  1:3627377/71‑1 (MQ=255)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGTTTACAg                                        <  1:4046631/38‑1 (MQ=25)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                        <  1:1539462/38‑1 (MQ=14)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                         >  1:2619381/1‑53 (MQ=11)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                         >  1:2140870/1‑53 (MQ=11)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                         >  1:2783897/1‑53 (MQ=11)

                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt      >  1:548387/1‑71 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    >  1:132964/1‑56 (MQ=35)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc       <  1:75550/69‑1 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg  >  1:4551920/1‑70 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg  >  1:1306844/1‑70 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/2359815‑2359934

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