BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	327,825	A→G	38.7%	R206R (CGT→CGC)	aes ←	acetyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	327,825	0	A	G	38.7%	26.0 / 40.7	62	R206R (CGT→CGC)	aes	acetyl esterase
Reads supporting (aligned to +/- strand): ref base A (24/14); new base G (9/15); total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAC  >  minE/327762‑327872

tGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGt                                            >  1:238689/1‑69 (MQ=255)

              tCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCTCCCAACAGACGGCGAGTCACGGAAt                                    >  1:2244661/1‑63 (MQ=255)

              tCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGt                               >  1:3306770/1‑68 (MQ=255)

              tCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTa                              >  1:227057/1‑69 (MQ=255)

              tCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTa                              >  1:2957301/1‑69 (MQ=255)

              tCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTa                              >  1:1310076/1‑69 (MQ=255)

               cTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGt                               <  1:2323698/67‑1 (MQ=255)

                tGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTaa                             >  1:2825954/1‑68 (MQ=255)

                tGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAAt                            >  1:647537/1‑69 (MQ=255)

                  cAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATc                           >  1:2119298/1‑68 (MQ=255)

                  cAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATc                           >  1:1975900/1‑68 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGAGGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:2020727/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:1763987/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCTTAAAGcc                  <  1:2966965/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:1016632/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:1726121/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:2250901/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:798670/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:814792/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:817862/69‑1 (MQ=255)

                          gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGcc                  <  1:960979/69‑1 (MQ=255)

                             gCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                <  1:490263/68‑1 (MQ=255)

                             gCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                <  1:543592/68‑1 (MQ=255)

                             gCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                <  1:1306863/68‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:3167626/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:3188300/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:2899617/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:610286/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:750203/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:2770935/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:757101/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:2371907/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:1220697/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:2119672/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:1999181/69‑1 (MQ=255)

                                       aTCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:1301139/69‑1 (MQ=255)

                                        tCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:797432/68‑1 (MQ=255)

                                        tCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa     <  1:581721/68‑1 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaaaa   >  1:923447/1‑68 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:1427511/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:693220/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:1269027/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:1189830/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:1168095/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:447197/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:1454683/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:862056/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:1619077/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:2944622/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:2240148/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:2625196/1‑69 (MQ=255)

                                          ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:927735/1‑69 (MQ=255)

                                          acAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAc  >  1:644387/2‑69 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATc                           >  1:2547935/1‑35 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTaaa                     >  1:1380480/1‑41 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                >  1:1266852/1‑46 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                >  1:1572908/1‑46 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCAt               >  1:1415277/1‑47 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCAt               >  1:1056618/1‑47 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCAt               >  1:1461084/1‑47 (MQ=255)

                                                   gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCAt               >  1:912104/1‑47 (MQ=255)

                                                         cAGACGACGAGTCACGGAATCCCGTAATCCGTATAGcc                  <  1:1026343/38‑1 (MQ=255)

                                                              gacgaGTCACGGAATCCCGTAATCCGTAAAGCCCa                <  1:2367591/35‑1 (MQ=255)

TGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAC  >  minE/327762‑327872

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