BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	40.6%	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	40.6%	101.5 / 160.0	133	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (45/34); new base T (31/23); total (76/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 133 total aligned reads displayed.

CGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGT  >  minE/2359817‑2359944

cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                           <  1:1330886/71‑1 (MQ=255)

     taacGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                       >  1:887778/4‑70 (MQ=11)

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

     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                       >  1:1418373/1‑70 (MQ=255)

       tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACATCACACCTTCACAGGCTTa                                                     >  1:3335971/1‑70 (MQ=25)

        cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                     >  1:323696/1‑69 (MQ=255)

         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                     >  1:4008091/1‑68 (MQ=255)

         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                     >  1:1434570/1‑68 (MQ=255)

         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                     >  1:2536669/1‑68 (MQ=255)

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

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

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

             cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                               <  1:2296680/70‑1 (MQ=255)

              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                  <  1:2075196/66‑1 (MQ=35)

              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                              >  1:3772624/1‑70 (MQ=35)

              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTGACAGTACACCTTCACAGGCTTACAg                                                  <  1:2979961/66‑1 (MQ=11)

              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                  <  1:2622659/66‑1 (MQ=255)

              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                  <  1:1578504/66‑1 (MQ=255)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                             <  1:4034623/69‑1 (MQ=255)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                             <  1:4037246/69‑1 (MQ=255)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                             <  1:3394110/69‑1 (MQ=255)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                             <  1:2865072/69‑1 (MQ=255)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                             <  1:1095512/69‑1 (MQ=255)

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

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

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

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

                      gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                     >  1:650659/1‑71 (MQ=34)

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

                         ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                  <  1:2973672/71‑1 (MQ=34)

                         ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                  <  1:1029874/71‑1 (MQ=34)

                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                   <  1:2986783/70‑1 (MQ=255)

                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                   <  1:672014/70‑1 (MQ=255)

                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                   <  1:1271121/70‑1 (MQ=255)

                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                   <  1:1502696/70‑1 (MQ=255)

                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                   <  1:1547257/70‑1 (MQ=255)

                          tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                              >  1:3547254/1‑58 (MQ=35)

                          tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                              >  1:2606740/1‑58 (MQ=35)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                             cACTTCTGATACCTCCAGAATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                <  1:459405/69‑1 (MQ=11)

                              aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                             <  1:1486207/71‑1 (MQ=34)

                               cTTCTGATACCTCCAGCATGTCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                            <  1:3998118/71‑1 (MQ=25)

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

                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                     >  1:3125457/1‑60 (MQ=35)

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

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

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

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

                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCt                                            >  1:1438513/1‑47 (MQ=38)

                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                           >  1:918568/1‑48 (MQ=38)

                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                           >  1:1970187/1‑48 (MQ=38)

                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                           >  1:3206190/1‑48 (MQ=38)

                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                           >  1:2742273/1‑48 (MQ=38)

                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                           >  1:2233623/1‑48 (MQ=38)

                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                           >  1:2317436/1‑48 (MQ=38)

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

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

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

                                         cTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                     >  1:1256572/1‑36 (MQ=37)

                                         cTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                     >  1:3811428/1‑36 (MQ=37)

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

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

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

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

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

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

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

                                                 tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc           <  1:1201955/70‑1 (MQ=255)

                                                   ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAg         >  1:2711819/1‑70 (MQ=255)

                                                   ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc        >  1:218495/1‑71 (MQ=255)

                                                   ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc        <  1:458544/71‑1 (MQ=255)

                                                   ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc        >  1:2173327/1‑71 (MQ=255)

                                                   ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc        >  1:372940/1‑71 (MQ=255)

                                                     tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg            >  1:3949218/1‑65 (MQ=255)

                                                     tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg            >  1:1278785/1‑65 (MQ=255)

                                                     tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg            >  1:2859021/1‑65 (MQ=255)

                                                       aaaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                              >  1:1103399/3‑45 (MQ=255)

                                                        caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCgg   <  1:2611431/71‑1 (MQ=255)

                                                        caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCgg   <  1:3349351/71‑1 (MQ=255)

                                                        caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCgg   <  1:877423/71‑1 (MQ=255)

                                                         aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                   <  1:2315367/54‑1 (MQ=255)

                                                         aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                   <  1:1657639/54‑1 (MQ=255)

                                                         aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                   <  1:597173/54‑1 (MQ=255)

                                                         aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGt  <  1:134639/71‑1 (MQ=255)

CGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGT  >  minE/2359817‑2359944

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`