BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→G	41.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	G	41.6%	156.7 / 186.1	192	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (87/7); new base G (65/10); total (160/32)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 192 total aligned reads displayed.

TATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTT  >  minE/2359822‑2359940

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:941792/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:1451251/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:3360407/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:2791972/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:1267829/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:2395271/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:1125299/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:3084566/1‑68 (MQ=11)

tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAgg                                                     >  1:1131114/1‑68 (MQ=11)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:2604245/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:569405/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:1217575/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:2427552/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:735322/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:149483/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:400037/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:1682582/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:1371674/1‑68 (MQ=35)

  tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:644065/68‑1 (MQ=255)

  tCGTTACTTATGTCAGCATTCGCACTTCTGAGACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                   >  1:3145594/1‑68 (MQ=25)

    gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTa                                                 <  1:340087/68‑1 (MQ=11)

    gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTa                                                 <  1:2862248/68‑1 (MQ=11)

    gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTa                                                 <  1:2288187/68‑1 (MQ=11)

     ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTAc                                                >  1:2077053/1‑68 (MQ=11)

     ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCAACCCTCACAGGCCACCTTCACAGGCTTAc                                                >  1:2599404/1‑68 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:1397196/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:1366094/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:609084/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:1231974/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:1173309/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:3029088/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:3424697/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:3046725/60‑1 (MQ=11)

          tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                   <  1:103759/60‑1 (MQ=11)

           aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                         >  1:2633202/1‑69 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      >  1:2698356/1‑69 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      >  1:2564163/1‑69 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      <  1:1012203/69‑1 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      <  1:3261700/69‑1 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      >  1:154547/1‑69 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      >  1:1333756/1‑69 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      >  1:878108/1‑69 (MQ=35)

              tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                      <  1:1145529/69‑1 (MQ=35)

                aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                    >  1:3245896/1‑69 (MQ=11)

                aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACATTCACAGGCTTACAGAACGCTcacc                                    >  1:1307533/1‑69 (MQ=11)

                    ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                 >  1:2557774/1‑68 (MQ=255)

                    ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                 >  1:1468303/1‑68 (MQ=255)

                      cGCACTTCTGATACCTCCAGCATGCCTCACTGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:1617826/1‑68 (MQ=255)

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

                      cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:465190/1‑68 (MQ=35)

                      cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:485687/1‑68 (MQ=35)

                      cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:3332045/1‑68 (MQ=35)

                      cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:3266157/1‑68 (MQ=35)

                      cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:714963/1‑68 (MQ=35)

                      cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:648698/1‑68 (MQ=35)

                      cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                               >  1:3218946/1‑68 (MQ=35)

                             cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                        >  1:3316545/1‑67 (MQ=11)

                             cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                        >  1:2037778/1‑67 (MQ=11)

                                    cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                    >  1:2510337/1‑65 (MQ=255)

                                    cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                >  1:1556064/1‑69 (MQ=255)

                                    cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                >  1:3108664/1‑69 (MQ=255)

                                    cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                >  1:3289554/1‑69 (MQ=255)

                                    cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                >  1:1757099/1‑69 (MQ=255)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:2807150/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:1795919/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:330571/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:1627691/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:782238/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:1893489/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:2643662/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:89436/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:2190161/1‑42 (MQ=14)

                                      ccAGCAACCCTCACAGGCCACCTTCACAGGCTTACAGAACGc                                         >  1:2151649/1‑42 (MQ=14)

                                          cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:2477373/1‑69 (MQ=255)

                                          cATGCCTCACAGCACACCTTCACAGGCGTACAGCACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:2218878/1‑69 (MQ=255)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:385253/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:2548108/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:10643/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:294518/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:1117140/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:1988576/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:3184260/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:1810867/3‑67 (MQ=17)

                                            acccTCACAGGCCACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc          >  1:1765904/3‑67 (MQ=17)

                                             gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc       <  1:1854423/69‑1 (MQ=255)

                                                 cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                          >  1:1162133/1‑46 (MQ=38)

                                                 cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                          >  1:3384706/1‑46 (MQ=38)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                >  1:2351360/1‑39 (MQ=40)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:1958958/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:1025369/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:1284873/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:135390/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:1437157/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:1570986/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:2042850/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:2715672/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:2763092/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:2858436/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:2978748/1‑69 (MQ=255)

                                                  acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt  >  1:3018668/1‑69 (MQ=255)

TATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTT  >  minE/2359822‑2359940

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