BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,359,876	C→T	32.2%	noncoding (1219/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,359,876	0	C	T	32.2%	160.6 / 113.5	121	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (50/32); new base T (16/23); total (66/55)
Fisher's exact test for biased strand distribution p-value = 5.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

Only 100 of 130 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGC  >  Exported/2359815‑2359946

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

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

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

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

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCtt                                                        >  1:2057340/2‑70 (MQ=11)

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCtt                                                        >  1:3217660/2‑70 (MQ=11)

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCtt                                                        >  1:4022429/2‑70 (MQ=11)

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCt                                                         >  1:1578302/2‑69 (MQ=11)

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

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

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                       <  1:1470147/69‑1 (MQ=255)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                     <  1:1590780/70‑1 (MQ=35)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                     <  1:4146891/70‑1 (MQ=35)

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

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                    <  1:612577/71‑1 (MQ=255)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                    <  1:4355287/71‑1 (MQ=255)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                                   <  1:3441075/71‑1 (MQ=34)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                    <  1:1411050/70‑1 (MQ=255)

             tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               >  1:2958641/1‑58 (MQ=31)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAAc                                                 >  1:1037015/1‑71 (MQ=34)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                  <  1:3019401/70‑1 (MQ=255)

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

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

                     tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                           >  1:1655544/1‑70 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                      >  1:1654707/1‑57 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                        >  1:1111928/1‑71 (MQ=34)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                        >  1:960316/1‑71 (MQ=34)

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

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

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               <  1:949872/47‑1 (MQ=32)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               <  1:2038990/47‑1 (MQ=32)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               <  1:2185006/47‑1 (MQ=32)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                             <  1:1020027/65‑1 (MQ=35)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                           >  1:2956437/1‑65 (MQ=35)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                     >  1:1925689/1‑70 (MQ=35)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                                     <  1:883135/54‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                                     <  1:705658/54‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                         <  1:3457662/66‑1 (MQ=255)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                                   <  1:2806292/55‑1 (MQ=35)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                                   <  1:1996104/55‑1 (MQ=35)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                   >  1:3927596/1‑71 (MQ=34)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                   >  1:1133134/1‑71 (MQ=34)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                   >  1:2815496/1‑71 (MQ=34)

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

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                       >  1:406779/1‑65 (MQ=255)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                       >  1:3984578/1‑65 (MQ=255)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                       >  1:2186018/1‑65 (MQ=255)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                <  1:2427504/71‑1 (MQ=34)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                <  1:2292556/71‑1 (MQ=34)

                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                        >  1:2854243/1‑61 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               <  1:156259/69‑1 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               <  1:3665446/69‑1 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               <  1:2249574/69‑1 (MQ=35)

                                   tCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                              >  1:1493806/1‑37 (MQ=11)

                                   tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                             >  1:3198516/1‑68 (MQ=11)

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

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                             <  1:1753502/69‑1 (MQ=39)

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

                                    cTGATACCTCCACCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                              <  1:2057722/68‑1 (MQ=255)

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

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

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg                         <  1:2685644/70‑1 (MQ=255)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                       <  1:2498028/51‑1 (MQ=11)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACTCTccc                                           >  1:3006355/1‑45 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                            >  1:699718/1‑60 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                            >  1:1618810/1‑60 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                            >  1:3673547/1‑60 (MQ=255)

                                               aGCATGCCTCACAGCACACCTTCACATGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc                >  1:2651377/1‑71 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                                <  1:552288/38‑1 (MQ=40)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                <  1:4282818/54‑1 (MQ=35)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc               >  1:1661703/1‑71 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                              >  1:4442921/1‑54 (MQ=38)

                                                     ccTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                          >  1:943124/1‑39 (MQ=255)

                                                     ccTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                          >  1:2650613/1‑39 (MQ=255)

                                                     ccTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                          >  1:1786483/1‑39 (MQ=255)

                                                     ccTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                          >  1:1259068/1‑39 (MQ=255)

                                                     ccTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                          >  1:65393/1‑39 (MQ=255)

                                                      cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               <  1:161238/49‑1 (MQ=255)

                                                      cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               <  1:4223839/49‑1 (MQ=255)

                                                      cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               <  1:897553/49‑1 (MQ=255)

                                                      cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCt         >  1:2351794/1‑71 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGACCGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt        <  1:2525004/69‑1 (MQ=255)

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

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

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

                                                          caGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                         <  1:2966485/35‑1 (MQ=255)

                                                          caGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                         <  1:3495687/35‑1 (MQ=255)

                                                          caGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                         <  1:124268/35‑1 (MQ=255)

                                                          caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                 >  1:857458/1‑59 (MQ=255)

                                                          caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                 >  1:182661/1‑59 (MQ=255)

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

                                                           aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGt    >  1:1790895/1‑71 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc                >  1:3565594/1‑57 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc                >  1:897781/1‑57 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc                >  1:3262194/1‑57 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg   >  1:3903094/1‑70 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg   >  1:343825/1‑70 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg   >  1:2578109/1‑70 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGc  <  1:18617/71‑1 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTA‑CCAACAACACATAGTGTCGCTGc                >  1:2699000/1‑56 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGC  >  Exported/2359815‑2359946

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