BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,444	C→T	18.5%	R73R (CGC→CGT)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,444	0	C	T	18.5%	118.3 / 15.2	75	R73R (CGC→CGT)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand): ref base C (23/38); new base T (9/5); total (32/43)
Fisher's exact test for biased strand distribution p-value = 8.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

ATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGT  >  minE/2086383‑2086482

aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:2356178/69‑1 (MQ=255)

aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:736083/69‑1 (MQ=255)

aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:2695811/69‑1 (MQ=255)

aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:2088646/69‑1 (MQ=255)

aTCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCACACCGCACGCTGATCATGCGTAAGCTCa                                 <  1:3341617/69‑1 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1136321/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2307985/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2287064/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2672779/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:2862125/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3132640/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1961304/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:176701/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3286326/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1629300/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3294415/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3380941/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:428489/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1368484/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:902186/1‑69 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCa                                 >  1:911021/1‑68 (MQ=255)

 tCGTTATCATCTCTATGTTTCACTCACCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:3209422/1‑69 (MQ=255)

  cGTTATAATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCaa                                >  1:1129350/1‑68 (MQ=255)

   gTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAg                              <  1:205515/69‑1 (MQ=255)

   gTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAg                              <  1:2695212/69‑1 (MQ=255)

     tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAgg                             <  1:2826920/68‑1 (MQ=255)

     tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAgg                             <  1:136365/68‑1 (MQ=255)

     tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGa                            >  1:3267349/1‑69 (MQ=255)

     tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGa                            >  1:2733958/1‑69 (MQ=255)

     tATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGa                            >  1:2034952/1‑69 (MQ=255)

                  tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:1309249/67‑1 (MQ=255)

                  tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:886676/67‑1 (MQ=255)

                  tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:672507/67‑1 (MQ=255)

                  tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:2899456/67‑1 (MQ=255)

                  tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:1609419/67‑1 (MQ=255)

                  tttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGtt                 <  1:3280258/67‑1 (MQ=255)

                   ttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAt              >  1:2847203/1‑69 (MQ=255)

                   ttCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTCTAt              >  1:75260/1‑69 (MQ=255)

                    tCACTCGCCTGCCCGTGGGCGCACCGTACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:2361702/69‑1 (MQ=255)

                    tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:3033062/69‑1 (MQ=255)

                    tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:1804032/69‑1 (MQ=255)

                    tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:2718161/69‑1 (MQ=255)

                    tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:895914/69‑1 (MQ=255)

                    tCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtt             <  1:2420371/69‑1 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATgcgct                                       <  1:952597/38‑2 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:346644/38‑1 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:2165364/38‑1 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:700801/38‑1 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:56226/38‑1 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:1380019/38‑1 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:340838/38‑1 (MQ=255)

                         cGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCa                                       <  1:2009319/38‑1 (MQ=255)

                                 cGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATttccgtttccgt  <  1:2264929/67‑1 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:787209/1‑58 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:1395387/1‑58 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:2798603/1‑58 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:399766/1‑58 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:768848/1‑58 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:2326437/1‑58 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccgt        >  1:295853/1‑58 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTCAAAGGACTGGAACCGTTTAtttccg         >  1:2560376/1‑57 (MQ=255)

                                    gggCACACCGCACGCTGATCATGCGTAAGCTAAAAGGACt                          >  1:695683/1‑40 (MQ=25)

                                     ggCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:3013778/62‑1 (MQ=255)

                                        gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:585069/54‑1 (MQ=255)

                                        gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:747816/54‑1 (MQ=255)

                                        gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:562192/54‑1 (MQ=255)

                                        gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:1551110/54‑1 (MQ=255)

                                        gcaccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgt        <  1:845790/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:639891/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:2216952/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:461391/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:2642759/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:1765953/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:1361624/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:2512966/54‑1 (MQ=255)

                                             gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTAtttccgtttccg   <  1:2030563/54‑1 (MQ=255)

ATCGTTATCATCTCTATGTTTCACTCGCCTGCCCGTGGGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGT  >  minE/2086383‑2086482

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