BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,690,717	G→A	32.1%	T202T (ACC→ACT)	hslO ←	heat shock protein Hsp33

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,690,717	0	G	A	32.1%	101.1 / 58.3	115	T202T (ACC→ACT)	hslO	heat shock protein Hsp33
Reads supporting (aligned to +/- strand): ref base G (64/14); new base A (10/27); total (74/41)
Fisher's exact test for biased strand distribution p-value = 1.55e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 115 total aligned reads displayed.

GATCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCCTGCTG  >  minE/2690656‑2690782

gATCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagca                                                            >  1:3069569/1‑69 (MQ=255)

gATCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagca                                                            >  1:296699/1‑69 (MQ=255)

  tCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGt                                                          >  1:551505/1‑69 (MQ=255)

  tCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGt                                                          >  1:2383891/1‑69 (MQ=255)

  tCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGt                                                          >  1:3397616/1‑69 (MQ=255)

  tCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGt                                                          >  1:2589930/1‑69 (MQ=255)

  tCGAAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagcag                                                           >  1:2089644/1‑68 (MQ=255)

   cGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGtt                                                         >  1:3214642/1‑69 (MQ=255)

   cGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGtt                                                         >  1:875993/1‑69 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2399500/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:3042610/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2869905/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2792377/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:330132/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2312693/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2248144/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2101573/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:1205761/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:1007177/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:872119/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:5819/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:1105330/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:417498/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:338449/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:354611/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:1316858/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:1320833/69‑1 (MQ=255)

      aaaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:3386130/69‑1 (MQ=255)

       aaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:284572/68‑1 (MQ=255)

       aaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2727062/68‑1 (MQ=255)

       aaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:2105536/68‑1 (MQ=255)

       aaCCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAAGTCAGCAGttctt                                                      <  1:3382178/68‑1 (MQ=255)

                 tcttctTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATg                                           >  1:1574782/1‑69 (MQ=255)

                 tcttctTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATg                                           >  1:1157112/1‑69 (MQ=255)

                 tcttctTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATg                                           >  1:449983/1‑69 (MQ=255)

                 tcttctTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATg                                           >  1:2167976/1‑69 (MQ=255)

                 tcttctTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGGTTTGATg                                           >  1:1744010/1‑69 (MQ=255)

                  cttcttCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATgg                                          >  1:341623/1‑69 (MQ=255)

                  cttcttCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATgg                                          >  1:94510/1‑69 (MQ=255)

                  cttcttCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATgg                                          >  1:177711/1‑69 (MQ=255)

                           gATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGtt                                  <  1:1052882/68‑1 (MQ=255)

                           gATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGtt                                  <  1:182012/68‑1 (MQ=255)

                           gATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGtt                                  <  1:2493391/68‑1 (MQ=255)

                           gATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGtt                                  <  1:3420414/68‑1 (MQ=255)

                           gATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGtt                                  <  1:897533/68‑1 (MQ=255)

                           gATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTa                                 >  1:1484650/1‑69 (MQ=255)

                           gATACAAACGCCACAACACTTCGTTGGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTa                                 >  1:3107856/1‑69 (MQ=255)

                              aCAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGc                               <  1:136952/68‑1 (MQ=255)

                              aCAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGc                               <  1:898047/68‑1 (MQ=255)

                               cAAACGTCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:435749/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTgagc                                                             >  1:789040/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:779679/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:623733/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:3286537/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:1666952/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:1675030/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:176848/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:3212231/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:2220777/1‑37 (MQ=255)

                               cAAACGCCACAACACTTCGTTTGCCGGTAAGGTcagc                                                             >  1:2224650/1‑37 (MQ=255)

                                       acaacaCTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggt                     >  1:353273/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCTGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:1394279/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:2733981/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:2182077/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:2937277/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:1354197/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:115580/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:1504099/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:587647/1‑69 (MQ=255)

                                        caacaCTTCGTTTGCCGGTAAAGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtg                    >  1:3280858/1‑69 (MQ=255)

                                          acacTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtgg                   >  1:252566/1‑68 (MQ=255)

                                           cacTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:1838154/1‑68 (MQ=255)

                                           cacTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:1603158/1‑68 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGcc                         >  1:147931/1‑60 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:1255330/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:2038786/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:1983275/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:1082603/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:6057/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:2235590/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:655734/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:810993/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:2463484/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAggtggt                  >  1:2987037/1‑67 (MQ=255)

                                            acTTCGTTTGCCGGTAAGGTCAGCAGATCTTCGGTTTTGATGGttt                                       >  1:1859457/1‑46 (MQ=255)

                                             cTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGc                               <  1:1156230/53‑1 (MQ=255)

                                              ttCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGtt                                  <  1:1914901/49‑1 (MQ=255)

                                              ttCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCaaa              >  1:390906/1‑69 (MQ=255)

                                              ttCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCaaa              >  1:2350606/1‑69 (MQ=255)

                                                   ttGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAgtcgt         <  1:209802/69‑1 (MQ=255)

                                                   ttGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAgtcgt         <  1:1082105/69‑1 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGt                                   >  1:1457205/1‑35 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAgtcgtc        >  1:304210/1‑62 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCctgctg  >  1:1126647/1‑68 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCctgctg  >  1:1255596/1‑68 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCctgctg  >  1:1312597/1‑68 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCctgctg  >  1:1058030/1‑68 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCctgctg  >  1:2546564/1‑68 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCctgctg  >  1:2602434/1‑68 (MQ=255)

                                                           aaGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGATAGCGTCGCCAGGTGGTCAAAGTCGTCctgctg  >  1:2224197/1‑68 (MQ=255)

GATCGTAAACCGTCACCTCTTCTTCGTGATACAAACGCCACAACACTTCGTTTGCCGGTAAGGTCAGCAGTTCTTCGGTTTTGATGGTTTCGGTTAGCGTCGCCAGGTGGTCAAAGTCGTCCTGCTG  >  minE/2690656‑2690782

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