BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,161,857	T→C	37.5%	C17C (TGT→TGC)	ydjY →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,161,857	0	T	C	37.5%	27.2 / 33.6	56	C17C (TGT→TGC)	ydjY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (26/9); new base C (8/13); total (34/22)
Fisher's exact test for biased strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTGGAA  >  minE/1161790‑1161920

aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                                <  1:1862192/69‑1 (MQ=255)

aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                                <  1:1644731/69‑1 (MQ=255)

aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                                <  1:3351372/69‑1 (MQ=255)

aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                                <  1:36486/69‑1 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTa                                                              >  1:3193100/1‑69 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTa                                                              >  1:2098896/1‑69 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTa                                                              >  1:2591298/1‑69 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTa                                                              >  1:636546/1‑69 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTa                                                              >  1:1163787/1‑69 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTa                                                              >  1:1053667/1‑69 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTa                                                              >  1:926406/1‑69 (MQ=255)

  aTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGACTGCACTCTGCTTa                                                              >  1:2400387/1‑69 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGTCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:2547734/69‑1 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:3404512/69‑1 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:606314/69‑1 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:961968/69‑1 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:3123782/69‑1 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:2874917/69‑1 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:2314296/69‑1 (MQ=255)

               cGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:1487053/69‑1 (MQ=255)

                gATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:3400200/68‑1 (MQ=255)

                gATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:2206877/68‑1 (MQ=255)

                gATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:1879386/68‑1 (MQ=255)

                gATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:1247305/68‑1 (MQ=255)

                gATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGCTTACTGGCTGTTGCTg                                                 <  1:495335/68‑1 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:2612133/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:3404462/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:1381981/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:336666/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:1400376/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:1429641/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:2974154/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:1336323/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCgg                                       >  1:2199044/1‑69 (MQ=255)

                         aCATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCg                                        >  1:1974284/1‑68 (MQ=255)

                                      tCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaacaa                          <  1:125556/69‑1 (MQ=255)

                                      tCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaacaa                          <  1:348080/69‑1 (MQ=255)

                                          ggAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAagag                       <  1:1722345/68‑1 (MQ=255)

                                          ggAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAagag                       <  1:2303480/68‑1 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAagag                       >  1:1633848/1‑67 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:565198/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:1006512/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:819046/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:866454/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:1018670/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:1986369/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:3379468/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:3118158/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:277288/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:274472/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:1639610/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:2471261/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:2279841/1‑69 (MQ=255)

                                           gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa                     >  1:2123376/1‑69 (MQ=255)

                                                              ctctGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTgg    <  1:40601/67‑1 (MQ=255)

                                                              ctctGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTgg    <  1:663067/67‑1 (MQ=255)

                                                              ctctGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTgg    <  1:663731/67‑1 (MQ=255)

                                                              ctctGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTgg    <  1:2400848/67‑1 (MQ=255)

                                                              ctctGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTgg    <  1:2229195/67‑1 (MQ=255)

                                                                 tGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTGGaa  >  1:2641879/1‑66 (MQ=255)

AAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGCGGCAAAAGTGGAA  >  minE/1161790‑1161920

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