BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	115,740	C→T	37.4%	L87L (CTG→CTA)	yacH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	115,740	0	C	T	37.4%	37.2 / 37.8	64	L87L (CTG→CTA)	yacH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (17/23); new base T (4/20); total (21/43)
Fisher's exact test for biased strand distribution p-value = 5.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TGCTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTG  >  minE/115674‑115807

tGCTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAg                                                                   <  1:1882584/69‑1 (MQ=255)

tGCTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAg                                                                   <  1:2697579/69‑1 (MQ=255)

tGCTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAg                                                                   <  1:170088/69‑1 (MQ=255)

tGCTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAg                                                                   <  1:1698732/69‑1 (MQ=255)

tGCTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAg                                                                   <  1:496453/69‑1 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:2529408/1‑68 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:2936513/1‑68 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:1493506/1‑68 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:2599771/1‑68 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:197103/1‑68 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:2348529/1‑68 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:911949/1‑68 (MQ=255)

  cTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACATGCGAAAGTAGGGCGTCGGGATACAgc                                                                  >  1:2131685/1‑68 (MQ=255)

               tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                     >  1:1185824/1‑68 (MQ=255)

               tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                     >  1:2380867/1‑68 (MQ=255)

               tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                     >  1:1913621/1‑68 (MQ=255)

               tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                     >  1:2019721/1‑68 (MQ=255)

               tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                     >  1:2296664/1‑68 (MQ=255)

               tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                     >  1:126973/1‑68 (MQ=255)

               tGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGa                                                     >  1:880100/1‑68 (MQ=255)

                               gCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAg                                                                   <  1:3283716/38‑1 (MQ=255)

                                  aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                  <  1:1927186/68‑1 (MQ=255)

                                  aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                  <  1:4942/68‑1 (MQ=255)

                                  aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                  <  1:2292634/68‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATATCTCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:2668343/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCATTGCTGAATg                       <  1:3387176/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:288219/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:2819282/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:3051442/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:3081135/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:3182131/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:654405/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:7026/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:2282317/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:2187023/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:1539974/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:2178574/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:2169835/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:1766218/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:1839019/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:1840732/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:1974296/69‑1 (MQ=255)

                                            gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:1945015/69‑1 (MQ=255)

                                              gAAAGTAGGGCGTCGGGGTATAGCGCGAGAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                       <  1:973161/67‑1 (MQ=255)

                                                   tAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGAttt                >  1:225097/1‑69 (MQ=255)

                                                    aGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTg               >  1:2525818/1‑69 (MQ=255)

                                                            ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc       >  1:2646863/1‑69 (MQ=255)

                                                            ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc       >  1:2349546/1‑69 (MQ=255)

                                                            ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc       >  1:1752683/1‑69 (MQ=255)

                                                            ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGAGCTGAATGCGGATTTGACGAcggc       >  1:1924490/1‑69 (MQ=255)

                                                            gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:2646368/41‑1 (MQ=255)

                                                            gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:911778/41‑1 (MQ=255)

                                                            gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:1161205/41‑1 (MQ=255)

                                                            gggATACAGCGCGA‑GGGCGCGACCCATTGATCTATTtgtg                                   <  1:3398690/40‑1 (MQ=39)

                                                              gATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:2672119/39‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:144692/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:1349907/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:551754/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:1437560/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:14412/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:2362157/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:1878231/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:1914786/69‑1 (MQ=255)

                                                                 aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:2446344/69‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:2086310/35‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:3407067/35‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:393848/35‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:1902935/35‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:1219932/35‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:1208716/35‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                   <  1:1193323/35‑1 (MQ=255)

                                                                  cAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg  <  1:1400908/68‑1 (MQ=255)

TGCTTGAGCAACGTTTGTCGGATAGGTTGATGCCATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTG  >  minE/115674‑115807

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