BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,402	G→A	53.5%	A188A (GCC→GCT)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,402	0	G	A	53.5%	7.2 / 86.3	99	A188A (GCC→GCT)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (21/25); new base A (38/15); total (59/40)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 101 total aligned reads displayed.

GCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTG  >  minE/125344‑125470

gctgctTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAAtt                                                             <  1:1368783/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTa                                                                >  1:128164/1‑46 (MQ=255)

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

                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAgt                                                    >  1:2902718/1‑57 (MQ=255)

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

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

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

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

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

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

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

                   aCGGAACACCCAGGCTTCTTTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCg                                         >  1:460851/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:1197319/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                    cttTCAGCTTGTGGACGTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGAttg                         <  1:3214108/68‑2 (MQ=255)

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

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

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

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

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

                                      ttCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTAc                        <  1:1650719/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:880417/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:79244/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                    gTTATCAGCGTTAATTTGCTTCAGTGGTGAAAAGCGTTGACCTTCCTGATTGCGACCATAGGCAGGCCa        >  1:2472213/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:2741277/69‑1 (MQ=255)

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

GCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCAGATTATGGACGTTATCGGCGTTAATTTGTTTCAGCGGCGAAAAGCGTTGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTG  >  minE/125344‑125470

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