BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	44.3%	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	44.3%	96.9 / 158.3	122	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (32/36); new base A (35/19); total (67/55)
Fisher's exact test for biased strand distribution p-value = 6.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 122 total aligned reads displayed.

CGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  minE/2359817‑2359936

cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                     >  1:904363/1‑69 (MQ=31)

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

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

     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                     >  1:2479075/1‑64 (MQ=31)

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

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

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

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

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

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

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

           tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:2572663/1‑71 (MQ=34)

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

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

                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)

                 tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                      >  1:353153/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:1036876/58‑1 (MQ=35)

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

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

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

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

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

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

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

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

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

                        aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                             >  1:218900/1‑69 (MQ=35)

                        aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                            >  1:1772510/1‑70 (MQ=255)

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

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

                         ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:1693522/55‑1 (MQ=35)

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

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

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

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

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

                           agCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                    >  1:2298846/2‑59 (MQ=11)

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

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

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

                             cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                      >  1:1412853/1‑71 (MQ=255)

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

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

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

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

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

                               cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                        <  1:2141356/67‑1 (MQ=255)

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

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

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

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

                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                          <  1:2628916/48‑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:2079094/70‑1 (MQ=40)

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

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

                                  cTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                   >  1:62419/1‑37 (MQ=14)

                                  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)

                                   tGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                        >  1:746094/1‑47 (MQ=11)

                                   tGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                         >  1:2419982/1‑46 (MQ=11)

                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTTCCCTACCCAACAACACATAgt              <  1:3719707/71‑1 (MQ=255)

                                     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)

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

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

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

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

                                           ccAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                        >  1:822748/1‑39 (MQ=14)

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

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

                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:2383528/1‑53 (MQ=35)

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

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

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

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

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

                                              gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt              <  1:3758762/62‑1 (MQ=255)

                                              gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt              <  1:286913/62‑1 (MQ=255)

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

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

                                               cATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                   <  1:3436196/56‑3 (MQ=11)

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

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

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

CGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  minE/2359817‑2359936

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