BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,011,433	T→C	33.0%	L259L (CTT→CTC)	yghZ →	aldo‑keto reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,011,433	0	T	C	33.0%	92.5 / 65.3	118	L259L (CTT→CTC)	yghZ	aldo‑keto reductase
Reads supporting (aligned to +/- strand): ref base T (56/23); new base C (30/9); total (86/32)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

Only 100 of 118 total aligned reads displayed.

ATTCCGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAA  >  minE/2011368‑2011497

aTTCCGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTa                                                                 >  1:256746/1‑67 (MQ=255)

aTTCCGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTAc                                                                >  1:3264314/1‑68 (MQ=255)

aTTCCGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTAc                                                                >  1:633157/1‑68 (MQ=255)

aTTCCGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTAc                                                                >  1:2607857/1‑68 (MQ=255)

aTTCCGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTAc                                                                >  1:1577463/1‑68 (MQ=255)

   ccGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGaa                                                            >  1:3425495/1‑69 (MQ=255)

   ccGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGaa                                                            >  1:1978969/1‑69 (MQ=255)

   ccGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGaa                                                            >  1:2814180/1‑69 (MQ=255)

   ccGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGaa                                                            >  1:1809636/1‑69 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:2378692/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:2437997/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:567893/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:1859420/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:2610859/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:741361/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:1778379/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:1689318/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:3205016/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:354249/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:796133/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:2884317/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:932823/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:1214510/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAcc                                                     >  1:3336963/1‑68 (MQ=255)

           ttCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAAc                                                      >  1:2896303/1‑67 (MQ=255)

                        cGTGAAGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCTAAGCCa                                                        <  1:442816/52‑1 (MQ=255)

                        cGTGAAGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCa                                                        <  1:2687842/52‑1 (MQ=255)

                        cGTGAAGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCa                                                        <  1:2147168/52‑1 (MQ=255)

                        cGTGAAGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCa                                                        <  1:489182/52‑1 (MQ=255)

                        cGTGAAGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCa                                                        <  1:1413499/52‑1 (MQ=255)

                        cGTGAAGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCa                                                        <  1:2859606/52‑1 (MQ=255)

                        cGTGAAGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCa                                                        <  1:2657809/52‑1 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTgcgc                                       >  1:2335650/1‑68 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:1410713/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:237940/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:3434474/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:488653/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:3179957/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:493902/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:3418277/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:1903247/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCATGCGCt                                      >  1:981327/1‑69 (MQ=255)

                         gTGAAGGGAATAAAGTTCGTGGTCTGAAACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      >  1:484196/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAATCTCAACAGCCTGCGCTTAt                                   >  1:259879/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:675053/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:1804864/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:1788245/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:1117654/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:43559/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:1257725/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:2796357/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:1356157/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:900897/1‑69 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      <  1:3060245/66‑1 (MQ=255)

                            aaGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      <  1:3130118/66‑1 (MQ=255)

                            aaGGGAATAAAATTCGTGGTCTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTAt                                   >  1:2277304/1‑69 (MQ=255)

                                   tAAAGTTCGTTGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCt                                      <  1:212780/59‑1 (MQ=255)

                                    aaaGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATtgaatgaa                           <  1:3076479/69‑1 (MQ=255)

                                      aGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAAt                         >  1:271000/1‑69 (MQ=255)

                                        ttCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATgg                       <  1:1573385/69‑1 (MQ=255)

                                          cGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGGGCTTATTGAATGAAATgg                       >  1:842024/1‑67 (MQ=255)

                                          cGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATgg                       >  1:2742934/1‑67 (MQ=255)

                                          cGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATgg                       >  1:1171590/1‑67 (MQ=255)

                                          cGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATgg                       >  1:426218/1‑67 (MQ=255)

                                              gTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCAcagc                 >  1:2781872/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATgg                       >  1:818583/1‑61 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:404734/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:2042054/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:2138041/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:1815250/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:455867/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:83838/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:860993/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:1307410/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:1027932/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:2938047/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:2941778/1‑69 (MQ=255)

                                                cTGACGCCGAAAATGCTCACCGAAGCCAACCTCAACAGCCTGCGCTTATTGCATGAAATGGCAcagcag               >  1:1049078/1‑69 (MQ=255)

                                                cTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCAcagcag               <  1:2953671/69‑1 (MQ=255)

                                                    cacCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCAcagcag               <  1:154030/65‑1 (MQ=255)

                                                        gAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAca       >  1:536542/1‑69 (MQ=255)

                                                        gAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAca       >  1:962964/1‑69 (MQ=255)

                                                        gAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAca       >  1:2350831/1‑69 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:2530651/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:3309865/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:340319/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:950475/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:2712681/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:2704697/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:1443006/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:811646/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:2454164/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:1841531/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:669442/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatc    <  1:1943260/69‑1 (MQ=255)

                                                           aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaata    <  1:649563/69‑2 (MQ=255)

                                                             tGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatcaa  >  1:3234731/1‑69 (MQ=255)

                                                             tGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatcaa  >  1:1427539/1‑69 (MQ=255)

                                                             tGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGAcaatcaa  >  1:1998855/1‑69 (MQ=255)

                                                               cTTACCGACGCCAACCTCAACAGCCTGCGCTTATTGAATGAACTGGCACAGCAGCGTGGAcaatc    <  1:2568919/65‑1 (MQ=255)

ATTCCGCAAGATTCACGGATGCATCGTGAAGGGAATAAAGTTCGTGGTCTGACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAA  >  minE/2011368‑2011497

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