BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	27.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	27.5%	227.2 / 96.6	127	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (61/31); new base A (25/10); total (86/41)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 132 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC  >  minE/2359815‑2359935

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    >  1:1726744/1‑71 (MQ=9)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:1118004/1‑70 (MQ=31)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:1935819/1‑70 (MQ=31)

  cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACCCCTTca                                                    >  1:1734077/1‑69 (MQ=11)

  cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    >  1:4051093/1‑69 (MQ=255)

    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                <  1:3648981/71‑1 (MQ=255)

     tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACTGGc                                               >  1:3240280/1‑71 (MQ=21)

     tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                >  1:3569664/1‑70 (MQ=35)

     tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                               >  1:3434232/1‑71 (MQ=34)

     tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                               >  1:960219/1‑71 (MQ=34)

     tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCAAAg                                                 >  1:2638166/1‑69 (MQ=14)

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACAGCTTCACAGGCt                                              >  1:3930823/1‑71 (MQ=21)

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              >  1:2566357/1‑71 (MQ=34)

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              >  1:3355606/1‑71 (MQ=34)

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              >  1:3360308/1‑71 (MQ=34)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                               <  1:1958852/69‑1 (MQ=35)

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

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

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                            >  1:701525/1‑71 (MQ=34)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                            >  1:3795828/1‑71 (MQ=34)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                            >  1:1159514/1‑71 (MQ=34)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                            >  1:3016827/1‑71 (MQ=34)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                            >  1:795577/1‑71 (MQ=34)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                            >  1:833078/1‑71 (MQ=34)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    >  1:2629704/1‑62 (MQ=255)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                            >  1:469311/1‑68 (MQ=35)

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

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

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                          <  1:3151965/70‑1 (MQ=255)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:320507/1‑59 (MQ=31)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:67762/1‑59 (MQ=31)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                 <  1:1612642/62‑1 (MQ=35)

            ttACTTATGTCAGAATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                 <  1:1310302/62‑1 (MQ=25)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                               >  1:3683982/1‑62 (MQ=35)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            >  1:3622071/1‑65 (MQ=255)

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

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

               cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                     <  1:1536827/71‑1 (MQ=255)

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

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

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

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

                  aTGTCAGAATTCGCACTTCTGATAGCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              <  1:1436536/59‑1 (MQ=11)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                             <  1:2445142/56‑1 (MQ=35)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                               >  1:3093076/1‑70 (MQ=35)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGAtt                                             <  1:2240645/56‑1 (MQ=255)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                              <  1:512093/70‑1 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                              <  1:497198/70‑1 (MQ=255)

                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                            >  1:2578589/1‑71 (MQ=255)

                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                           >  1:535871/1‑71 (MQ=34)

                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                           >  1:3631865/1‑71 (MQ=34)

                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                           >  1:1438667/1‑71 (MQ=34)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          <  1:1849008/71‑1 (MQ=34)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          <  1:3669979/71‑1 (MQ=34)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          <  1:4020747/71‑1 (MQ=34)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          <  1:2025525/71‑1 (MQ=34)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                <  1:244666/48‑1 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                             >  1:3279903/1‑49 (MQ=38)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                             >  1:1587839/1‑49 (MQ=38)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                         >  1:3154535/1‑53 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    >  1:2653493/1‑42 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                              >  1:1726614/1‑48 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                              >  1:3507741/1‑48 (MQ=11)

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

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

                              gCACTTCTGATACCTCTAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                       <  1:3804444/70‑1 (MQ=11)

                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                       <  1:483242/70‑1 (MQ=35)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                         <  1:1284106/50‑1 (MQ=38)

                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                   >  1:1310802/1‑39 (MQ=37)

                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                   <  1:1396406/71‑1 (MQ=35)

                                   tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                        >  1:1091642/1‑64 (MQ=255)

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

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

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCAAAGGCTTACAGAACGCt                                   >  1:3550474/1‑52 (MQ=255)

                                      gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg              >  1:3956241/1‑71 (MQ=40)

                                      gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg              >  1:968737/1‑71 (MQ=40)

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt             >  1:3293890/1‑71 (MQ=255)

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt             <  1:2496302/71‑1 (MQ=255)

                                         aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                   >  1:1956228/1‑62 (MQ=11)

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

                                           cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg         >  1:357564/1‑71 (MQ=255)

                                           cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg         >  1:408283/1‑71 (MQ=255)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc        >  1:2154726/1‑71 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                >  1:3952570/1‑61 (MQ=40)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                >  1:675944/1‑61 (MQ=40)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt           >  1:3478906/1‑66 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt           >  1:3050707/1‑66 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt       <  1:68014/70‑1 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                  <  1:2875922/57‑1 (MQ=40)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc    <  1:1875040/71‑1 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                             >  1:1338726/1‑44 (MQ=38)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                             >  1:895311/1‑44 (MQ=38)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                             >  1:2557123/1‑44 (MQ=38)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg            <  1:2770168/61‑1 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc  <  1:1534568/71‑1 (MQ=255)

                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg   >  1:1602477/1‑69 (MQ=255)

                                                    gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                >  1:2377558/1‑39 (MQ=40)

                                                    gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                >  1:2838494/1‑39 (MQ=40)

                                                    gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                >  1:2912351/1‑39 (MQ=40)

                                                    gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                >  1:3146564/1‑39 (MQ=40)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC  >  minE/2359815‑2359935

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