BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	22.7%	noncoding (1219/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,876	0	C	T	22.7%	355.9 / 99.7	163	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (91/35); new base T (25/12); total (116/47)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

Only 100 of 164 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359815‑2359945

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:398376/70‑1 (MQ=35)

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

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

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

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

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

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

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

              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)

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

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

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

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

                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                                       >  1:654355/1‑71 (MQ=255)

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

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

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

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

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

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

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    >  1:3281033/1‑70 (MQ=255)

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

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

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

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

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

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

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:2580725/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)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                                <  1:2084923/67‑1 (MQ=255)

                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:3107433/1‑71 (MQ=39)

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

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

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

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

                                         aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                   <  1:3508259/57‑1 (MQ=11)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                   taccTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:2218841/47‑1 (MQ=11)

                                                   taccTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:209900/47‑1 (MQ=11)

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

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

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

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

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

                                                     ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1797157/1‑46 (MQ=38)

                                                       tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCTCAGCtt       >  1:2884499/1‑71 (MQ=255)

                                                        cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:327475/1‑50 (MQ=255)

                                                        cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:2233203/1‑71 (MQ=255)

                                                        cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:2228416/1‑71 (MQ=255)

                                                        cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:296001/1‑71 (MQ=255)

                                                        cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:3193771/1‑71 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                       >  1:907942/1‑37 (MQ=40)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:2774547/1‑38 (MQ=40)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      <  1:1174947/70‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      <  1:1853080/70‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:4086759/1‑70 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     >  1:3603628/1‑71 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTAACCAACAACACATAGTGTCGCTGCCGCAGCTTc      <  1:601285/70‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCATACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     >  1:958113/1‑71 (MQ=255)

                                                          caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCgg    >  1:2299927/1‑71 (MQ=255)

                                                            gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                         >  1:2061349/1‑48 (MQ=255)

                                                            gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                         >  1:3465377/1‑48 (MQ=255)

                                                            gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                         >  1:2930303/1‑48 (MQ=255)

                                                            gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                       <  1:3000382/50‑1 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359815‑2359945

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