BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,447	A→G	46.4%	R173R (CGT→CGC)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,447	0	A	G	46.4%	6.7 / 57.6	69	R173R (CGT→CGC)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (22/15); new base G (30/2); total (52/17)
Fisher's exact test for biased strand distribution p-value = 1.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

TTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGC  >  minE/125382‑125514

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

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

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

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

ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                                                                    <  1:295906/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:3174308/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                 aTCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTc                                                 >  1:3268000/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:397131/1‑69 (MQ=255)

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

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

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

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

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

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

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

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCgg                                         >  1:1370818/1‑43 (MQ=255)

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGca                                       >  1:51718/1‑44 (MQ=255)

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc               >  1:2395529/1‑68 (MQ=255)

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc               >  1:2802582/1‑68 (MQ=255)

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc               >  1:1891082/1‑68 (MQ=255)

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc               >  1:1884412/1‑68 (MQ=255)

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc               >  1:776173/1‑68 (MQ=255)

                                                   tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTgg                >  1:1444255/1‑68 (MQ=255)

                                                    gACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggcg              <  1:949902/69‑3 (MQ=255)

                                                    gACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggcg              <  1:3258584/69‑3 (MQ=255)

                                                    gACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc                                        >  1:1554824/1‑43 (MQ=255)

                                                    gACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc                                        >  1:2590529/1‑43 (MQ=255)

                                                    gACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc                                        >  1:1038259/1‑43 (MQ=255)

                                                    gACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc                                        >  1:227815/1‑43 (MQ=255)

                                                    gACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc                                        >  1:1750026/1‑43 (MQ=255)

                                                    gACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc                                        >  1:442886/1‑43 (MQ=255)

                                                    gACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGc                                        >  1:1423864/1‑43 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTgg           >  1:1693892/1‑60 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:631306/1‑68 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:1333462/1‑68 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:703967/1‑68 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:749228/1‑68 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:1917902/1‑68 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:1146608/1‑68 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:2446655/1‑68 (MQ=255)

                                                                tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg   >  1:1532830/1‑68 (MQ=255)

                                                                 aCGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcgc  >  1:1745596/1‑68 (MQ=255)

TTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGC  >  minE/125382‑125514

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