BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,900	C→T	43.0%	G225G (GGC→GGT)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,900	0	C	T	43.0%	17.7 / 40.7	58	G225G (GGC→GGT)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand): ref base C (20/13); new base T (16/9); total (36/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGA  >  minE/2086836‑2086941

aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAg                                        >  1:2862762/1‑68 (MQ=255)

aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:1092562/1‑69 (MQ=255)

aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:876641/1‑69 (MQ=255)

aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:254596/1‑69 (MQ=255)

aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:1836495/1‑69 (MQ=255)

aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGc                                       >  1:591809/1‑69 (MQ=255)

    ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2570840/1‑68 (MQ=255)

    ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:358267/1‑68 (MQ=255)

    ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2482359/1‑68 (MQ=255)

    ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2268366/1‑68 (MQ=255)

    ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:2263092/1‑68 (MQ=255)

    ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                    >  1:1800903/1‑68 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:2038609/68‑1 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:2107965/68‑1 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                <  1:1824154/68‑1 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1800644/1‑69 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1692510/1‑69 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:3240772/1‑69 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:123991/1‑69 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1551505/1‑69 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1547892/1‑69 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:1504772/1‑69 (MQ=255)

        gacgaGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTAc                               >  1:724597/1‑69 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:639723/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:337538/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2574738/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:783654/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:1393137/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:24715/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2391327/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:2239759/69‑1 (MQ=255)

                 gTGGCGAAAGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGc                      <  1:366354/69‑1 (MQ=255)

                   ggCGAAAGTGTTTGAATCGCTGGCACGACTCGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCaa                    <  1:3276427/69‑1 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCa                 >  1:845534/1‑67 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:351842/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3404085/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2541444/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:366492/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3357949/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:519403/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3072137/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2595199/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:2588373/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1010921/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1118078/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:3352757/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:180907/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1766405/1‑68 (MQ=255)

                        aaGTGTTTAAATCGCTGGCGCGACTGGAACAGATTTTAGGTCAGCACCGTTACCTGACCGGCAACCAg                >  1:1327848/1‑68 (MQ=255)

                                 aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:930283/68‑1 (MQ=255)

                                 aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:2104850/68‑1 (MQ=255)

                                 aaTCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa       <  1:851755/68‑1 (MQ=255)

                                     gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1931762/69‑1 (MQ=255)

                                     gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1309153/69‑1 (MQ=255)

                                     gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:1394879/69‑1 (MQ=255)

                                     gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:473901/69‑1 (MQ=255)

                                     gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa  <  1:2030927/69‑1 (MQ=255)

                                     gCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCCGCAACCAGCTAACCGAAGCCGa  <  1:3099656/69‑1 (MQ=255)

AAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGA  >  minE/2086836‑2086941

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`