BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,081,481	T→C	27.4%	D392D (GAT→GAC)	exuT →	hexuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,081,481	0	T	C	27.4%	118.6 / 45.9	109	D392D (GAT→GAC)	exuT	hexuronate transporter
Reads supporting (aligned to +/- strand): ref base T (24/55); new base C (15/15); total (39/70)
Fisher's exact test for biased strand distribution p-value = 7.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GTGTATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCGCATG  >  minE/2081421‑2081537

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:1828961/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:1867880/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:757538/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:2071409/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:3107374/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:1309537/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:2220682/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:2808409/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:2616781/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCgg                                                  >  1:241163/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTgg                                                  >  1:125993/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTgg                                                  >  1:3331605/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTgg                                                  >  1:1342361/1‑69 (MQ=255)

gtgtATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCAGAtg                                                         >  1:1632348/1‑62 (MQ=255)

   tATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCGGTc                                                >  1:900059/1‑68 (MQ=255)

   tATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCGGTc                                                >  1:2118432/1‑68 (MQ=255)

   tATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCGGTc                                                >  1:2993862/1‑68 (MQ=255)

   tATCGGTGGTTTTGCTCACCAGGCACTGTCTGGCGCACTGATTACGCTCTCTTCTGACGTGTTCGGTc                                                >  1:8392/1‑68 (MQ=255)

   tATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTc                                                >  1:1526877/1‑68 (MQ=255)

   tATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTc                                                >  1:2648149/1‑68 (MQ=255)

   tATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTc                                                >  1:1686225/1‑68 (MQ=255)

       ggttgtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:3320665/1‑68 (MQ=255)

       ggtggtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:2477678/1‑68 (MQ=255)

       ggtggtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:2318608/1‑68 (MQ=255)

       ggtggtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:2939236/1‑68 (MQ=255)

       ggtggtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:1988075/1‑68 (MQ=255)

       ggtggtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:972131/1‑68 (MQ=255)

       ggtggtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:1566260/1‑68 (MQ=255)

       ggtggtTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTaa                                            >  1:1209940/1‑68 (MQ=255)

                      caGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            >  1:1168588/1‑69 (MQ=255)

                      caGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            >  1:710613/1‑69 (MQ=255)

                      caGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            >  1:1322873/1‑69 (MQ=255)

                      caGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            >  1:308785/1‑69 (MQ=255)

                      caGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            >  1:3097923/1‑69 (MQ=255)

                      aaGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            >  1:309837/2‑69 (MQ=255)

                             tGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            <  1:781661/62‑1 (MQ=255)

                             tGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            <  1:2450578/62‑1 (MQ=255)

                             tGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            <  1:2350468/62‑1 (MQ=255)

                             tGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCa                            <  1:3217643/62‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAAcctg              <  1:2064335/68‑3 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:1751294/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:3047598/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:235173/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:1913754/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:895870/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:526850/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:619446/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:2742644/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:1842366/68‑1 (MQ=255)

                                     gcgcTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:2687447/68‑1 (MQ=255)

                                        cTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATg          <  1:730705/69‑1 (MQ=255)

                                        cTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATg          <  1:36566/69‑1 (MQ=255)

                                        cTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATg          <  1:295550/69‑1 (MQ=255)

                                        cTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATg          <  1:2715293/69‑1 (MQ=255)

                                        cTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATg          <  1:1762067/69‑1 (MQ=255)

                                        cTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATg          <  1:2413410/69‑1 (MQ=255)

                                        cTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATg          <  1:157349/69‑1 (MQ=255)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:1298611/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:3196470/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:843209/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:130276/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:1711200/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:1961659/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:1652639/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:3247277/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:2427677/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:596988/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:609944/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:2521632/35‑1 (MQ=25)

                                             tACGCTCTCTTCTGACGTGTTCGGTCGTAACGAAg                                       <  1:1429069/35‑1 (MQ=25)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:1058772/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:3286756/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:3361979/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:414138/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:496727/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:1265906/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:3172421/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:1770658/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:2245233/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:1468695/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:266409/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:3250955/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:2816305/55‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:1394615/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:1771527/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:1979694/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:2117610/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:787002/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:2846206/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:2901415/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:1312523/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:308845/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:3114421/68‑1 (MQ=255)

                                             tACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:1267036/68‑1 (MQ=255)

                                             tACGCCCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:3433229/68‑1 (MQ=255)

                                              aCGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCg      <  1:3165976/67‑1 (MQ=255)

                                               cGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCgg              <  1:1422660/58‑1 (MQ=255)

                                                 ctctctTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCTCAtg  <  1:627973/68‑1 (MQ=255)

                                                        ccGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:453981/44‑1 (MQ=255)

                                                        ccGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTa                   <  1:1785027/44‑1 (MQ=255)

GTGTATCGGTGGTTTTGCCCACCAGGCACTGTCTGGCGCGCTGATTACGCTCTCTTCCGATGTGTTTGGTCGTAACGAAGTGGCAACGGCAAACGGTTTAACCGGGATGTCCGCATG  >  minE/2081421‑2081537

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