BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	38.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	38.7%	166.8 / 158.8	142	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (43/44); new base T (36/19); total (79/63)
Fisher's exact test for biased strand distribution p-value = 8.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 143 total aligned reads displayed.

CGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359817‑2359945

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

  cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                          >  1:3837865/1‑71 (MQ=255)

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

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

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

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

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

              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                             <  1:3421310/56‑1 (MQ=33)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACCCCTTCACAGGCTTACAGAACg                                               >  1:3011218/1‑68 (MQ=11)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:3757847/1‑58 (MQ=11)

                aTGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCTTACAGAACGCTc                                            >  1:316016/1‑71 (MQ=11)

                 tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                               >  1:162628/1‑67 (MQ=255)

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

                   tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:3066581/58‑1 (MQ=35)

                    cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:916242/60‑1 (MQ=35)

                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    <  1:3683391/58‑1 (MQ=35)

                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    <  1:2020419/58‑1 (MQ=35)

                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    <  1:1036876/58‑1 (MQ=35)

                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:3461064/59‑1 (MQ=35)

                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:352980/59‑1 (MQ=35)

                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 >  1:222178/1‑61 (MQ=35)

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

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

                      gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      <  1:649386/71‑1 (MQ=34)

                      gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCCTACAGAACGCTCCCCTAc                                      <  1:2371334/71‑1 (MQ=21)

                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:707724/58‑1 (MQ=11)

                       cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      <  1:3009409/70‑1 (MQ=255)

                         ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                     >  1:3792520/1‑53 (MQ=35)

                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    >  1:1699124/1‑52 (MQ=35)

                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:29054/53‑1 (MQ=35)

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

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

                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            >  1:2527005/1‑60 (MQ=11)

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

                            gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:765274/49‑1 (MQ=38)

                            gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCt                                             <  1:1315852/58‑1 (MQ=35)

                             cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    <  1:3704037/66‑1 (MQ=255)

                              aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      >  1:1748979/1‑47 (MQ=38)

                              aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:1951112/1‑68 (MQ=35)

                              aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        >  1:43282/1‑45 (MQ=11)

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

                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                             >  1:306075/1‑71 (MQ=35)

                               cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      <  1:3904080/62‑1 (MQ=255)

                               cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                             >  1:3269095/1‑70 (MQ=255)

                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:2941641/48‑1 (MQ=11)

                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:2628916/48‑1 (MQ=11)

                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:559555/48‑1 (MQ=11)

                                 tttGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:2456192/45‑1 (MQ=11)

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

                                 tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:281792/1‑49 (MQ=11)

                                 tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:2368829/1‑49 (MQ=11)

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

                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1015562/1‑63 (MQ=35)

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

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

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

                                  cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:1051294/1‑48 (MQ=11)

                                  cTGATACATCCAGCATGCCTCACATCACACCTTCACAGGCTTACAGAACGCTc                                            >  1:2666345/1‑53 (MQ=11)

                                  cTGATACATCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1505053/1‑63 (MQ=255)

                                   tGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:875362/1‑69 (MQ=40)

                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                          <  1:1854277/68‑1 (MQ=40)

                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                          <  1:809365/68‑1 (MQ=40)

                                        ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    <  1:2243065/39‑1 (MQ=40)

                                         cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               <  1:1086609/59‑1 (MQ=35)

                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    <  1:2123321/53‑1 (MQ=35)

                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               <  1:2897072/58‑1 (MQ=35)

                                          tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            <  1:2225354/45‑1 (MQ=11)

                                           ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   >  1:2030117/1‑37 (MQ=40)

                                           ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:617099/1‑61 (MQ=40)

                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGTACGCt                                             <  1:4005586/42‑1 (MQ=255)

                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:2454654/1‑49 (MQ=38)

                                            cAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                             >  1:171729/1‑42 (MQ=14)

                                            cAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                             >  1:2680535/1‑42 (MQ=14)

                                            cAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                             >  1:1544702/1‑42 (MQ=14)

                                            cAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                           <  1:3648493/44‑1 (MQ=11)

                                            cAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:662250/1‑56 (MQ=11)

                                            aaGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                             >  1:384816/2‑42 (MQ=14)

                                             aGCATACCTCACAGTACACCTTCACAGGGTTACAg                                                   <  1:2545623/35‑1 (MQ=25)

                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg                        >  1:1825035/1‑60 (MQ=255)

                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg                        >  1:3375560/1‑60 (MQ=255)

                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg                   <  1:113118/65‑1 (MQ=255)

                                                 tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            <  1:2711933/38‑1 (MQ=40)

                                                 tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt                 <  1:3456776/65‑1 (MQ=255)

                                                 tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCa           <  1:1628835/71‑1 (MQ=255)

                                                   ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:4018070/1‑53 (MQ=255)

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

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

                                                     tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:2657063/45‑1 (MQ=38)

                                                     tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                            <  1:3432506/50‑1 (MQ=255)

                                                     tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                            <  1:1460469/50‑1 (MQ=255)

                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:444702/1‑49 (MQ=255)

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

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

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

                                                       acaGCACACCTTAACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:2179401/43‑1 (MQ=255)

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

                                                         aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc               <  1:2587804/59‑1 (MQ=255)

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

                                                          gCACACCTTCACATGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg  >  1:3431880/1‑71 (MQ=255)

                                                           cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                    <  1:3780511/52‑1 (MQ=255)

CGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359817‑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: `‑`