BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,414	T→C	47.7%	K184K (AAA→AAG)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,414	0	T	C	47.7%	2.6 / 98.8	109	K184K (AAA→AAG)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (23/34); new base C (37/15); total (60/49)
Fisher's exact test for biased strand distribution p-value = 1.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTA  >  minE/125361‑125478

gTACGGAACACCCAGGCTTCTTTCAGATTATTGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:2765523/69‑1 (MQ=255)

gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:1970540/69‑1 (MQ=255)

gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:237929/69‑1 (MQ=255)

gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:536265/69‑1 (MQ=255)

gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:1448472/69‑1 (MQ=255)

gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:358125/69‑1 (MQ=255)

gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:2740856/69‑1 (MQ=255)

gTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAg                                                   <  1:3212579/69‑1 (MQ=255)

  aCGGAACCCCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:882678/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCGGTGGTGAAAAGCg                                                 >  1:2630930/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGc                                                  >  1:859148/1‑68 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:28759/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2826107/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:3065593/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:3380948/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2576589/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:460851/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2494485/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:2377523/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:1056445/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:1045577/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                                 >  1:1119219/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCa                                                              >  1:1259497/1‑56 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:792122/1‑69 (MQ=255)

  aCGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:2426253/1‑69 (MQ=255)

   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1904489/1‑68 (MQ=255)

   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1197319/1‑68 (MQ=255)

   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:2983539/1‑68 (MQ=255)

   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:228711/1‑68 (MQ=255)

   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:2305292/1‑68 (MQ=255)

   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1387996/1‑68 (MQ=255)

   cGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCg                                                 >  1:1499356/1‑68 (MQ=255)

      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:2625045/69‑1 (MQ=255)

      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:3146023/69‑1 (MQ=255)

      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:66677/69‑1 (MQ=255)

      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:2680212/69‑1 (MQ=255)

      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:810253/69‑1 (MQ=255)

      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:2358052/69‑1 (MQ=255)

      aaCACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGa                                             <  1:1971037/69‑1 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGt                                                                  >  1:2565070/1‑38 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                                               >  1:1220988/1‑41 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                                               >  1:1076942/1‑41 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                                               >  1:2825725/1‑41 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                                               >  1:2947691/1‑41 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                                               >  1:1251358/1‑41 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                                               >  1:1477887/1‑41 (MQ=255)

                cttcttTCAGATTATGGACGTTATCGGCGTTAATTTGTTTc                                                               >  1:68826/1‑41 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1942959/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1521373/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:299409/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:3221688/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:325195/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2758796/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:689254/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:547953/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2434310/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2499633/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2592203/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:2588433/69‑2 (MQ=255)

                  tcttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1402113/69‑2 (MQ=255)

                     ttCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                                 <  1:1353458/66‑2 (MQ=255)

                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:2330705/67‑1 (MQ=255)

                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:3369575/67‑1 (MQ=255)

                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:1650719/67‑1 (MQ=255)

                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:295906/67‑1 (MQ=255)

                     ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                <  1:2620457/67‑1 (MQ=255)

                                   gTTATCAGCGTTGATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2418676/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATaggcaggcga                >  1:1377747/1‑67 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATag                        >  1:3277944/1‑61 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:3315230/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:675752/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:371806/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:880417/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2905413/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:37225/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:79244/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2783464/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:200541/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2219600/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2472213/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:261269/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:2695279/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:3023737/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa                >  1:3174308/1‑69 (MQ=255)

                                   gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGAGACCATAGGCAGGCCa                >  1:499470/1‑69 (MQ=255)

                                      aTCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATa                         >  1:3383527/1‑57 (MQ=255)

                                      aTCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTc             >  1:3268000/1‑69 (MQ=255)

                                      aTCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTc             >  1:2107416/1‑69 (MQ=255)

                                         ggCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTg          <  1:2118272/69‑1 (MQ=255)

                                         ggCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTg          <  1:2741277/69‑1 (MQ=255)

                                          gCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTg          <  1:1379588/68‑1 (MQ=255)

                                               aaTTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc    <  1:12607/69‑1 (MQ=255)

                                                 tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa  >  1:1073343/1‑69 (MQ=255)

                                                 tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa  >  1:151714/1‑69 (MQ=255)

                                                 tttGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTa  >  1:397131/1‑69 (MQ=255)

                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:3310131/49‑1 (MQ=255)

                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:1764743/49‑1 (MQ=255)

                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:863758/49‑1 (MQ=255)

                                                   tGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATaggcag                    <  1:791661/49‑1 (MQ=255)

                                                   tGTTTCAGCGGCGAAAAGCGTTGAACTTCCTGATTACGACCATagggag                    <  1:3292707/49‑1 (MQ=38)

GTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTA  >  minE/125361‑125478

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