BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,230,986	C→T	30.3%	E590E (GAG→GAA)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,230,986	0	C	T	30.3%	66.0 / 29.3	69	E590E (GAG→GAA)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base C (21/27); new base T (8/13); total (29/40)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGA  >  minE/1230920‑1231051

cGTTACTAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:426566/69‑1 (MQ=255)

cGTTACCAGTAGATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:1129413/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:2215334/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:2713461/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:1259863/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:1333612/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:2810470/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:1804705/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:3007211/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTCTAGCTTTTCATATATTc                                                                 <  1:2536628/69‑1 (MQ=255)

cGTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGGCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:3363092/69‑1 (MQ=255)

 gTTACCAGTAAATGCGGATAAACCTGTGCAGTGACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:645242/68‑1 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:3212283/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:1031354/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:29029/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:2962133/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:2964264/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:300809/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:3355271/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:36888/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:859086/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:936107/1‑69 (MQ=255)

           aaTGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCg                                                      >  1:1062149/1‑69 (MQ=255)

             tGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGa                                                    <  1:2876088/69‑1 (MQ=255)

             tGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGa                                                    <  1:2565576/69‑1 (MQ=255)

                 gATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATcc                                                 <  1:1351879/68‑1 (MQ=255)

                 gATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATcc                                                 <  1:1551021/68‑1 (MQ=255)

                     aGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCgg                                            >  1:792608/1‑69 (MQ=255)

                     aGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCgg                                            >  1:2913424/1‑69 (MQ=255)

                     aGCCTGAGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCgg                                            >  1:5591/1‑69 (MQ=255)

                                 gACGTTGTCATACGGGCTGTAGCTTTTCATATATTc                                                                 <  1:596762/36‑1 (MQ=255)

                                     ttGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                            <  1:979799/69‑1 (MQ=255)

                                     ttGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                            <  1:659586/69‑1 (MQ=255)

                                     ttGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                            <  1:157407/69‑1 (MQ=255)

                                     ttGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                            <  1:1710792/69‑1 (MQ=255)

                                     ttGTCATACGGGCTGTAGCCTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                            <  1:3335510/69‑1 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaa                           <  1:371354/69‑1 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaa                           <  1:171047/69‑1 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCCCCACTCTTCaa                           >  1:2883428/1‑69 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCCCCACTCTTCaa                           >  1:663033/1‑69 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCCCCACTCTTCaa                           >  1:1316868/1‑69 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCCCCACTCTTCaa                           >  1:385119/1‑69 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCCCCACTCTTCaa                           >  1:355499/1‑69 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCCCCACTCTTCaa                           >  1:3225705/1‑69 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCCCCACTCTTCaa                           >  1:2827709/1‑69 (MQ=255)

                                      tGTCATACGGGCTGTAGCTTTTCATATATTCGTAATATTGCGGATCCTGCGGATTCCACCACTCTTCaa                           >  1:2505832/1‑69 (MQ=255)

                                       gTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaaa                          >  1:1076156/1‑69 (MQ=255)

                                       gTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaaa                          >  1:893697/1‑69 (MQ=255)

                                       gTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaa                           >  1:2220853/1‑68 (MQ=255)

                                          aTACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTc                       <  1:3167965/69‑1 (MQ=255)

                                          aTACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTc                       <  1:814121/69‑1 (MQ=255)

                                          aTACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTc                       <  1:1945975/69‑1 (MQ=255)

                                          aTACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTc                       <  1:2635969/69‑1 (MQ=255)

                                          aTACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTc                       <  1:50021/69‑1 (MQ=255)

                                          aTACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTc                       <  1:1352640/69‑1 (MQ=255)

                                          aTACGGGCTGTAGCTTTTCATGCACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTc                       <  1:1946719/69‑1 (MQ=255)

                                            aCGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAc                     <  1:2761842/69‑1 (MQ=255)

                                            aCGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAc                     <  1:2134193/69‑1 (MQ=255)

                                                  tGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTg                >  1:2979416/1‑68 (MQ=255)

                                                  tGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTg                >  1:1667657/1‑68 (MQ=255)

                                                  tGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTg                >  1:2449088/1‑68 (MQ=255)

                                                  tGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTg                >  1:3205099/1‑68 (MQ=255)

                                                         tttCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAgg        <  1:1141346/69‑1 (MQ=255)

                                                         tttCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAgg        <  1:902689/69‑1 (MQ=255)

                                                               gTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCGAACTCACCAGTGGTAAGAGGAATTGa  <  1:1553238/69‑1 (MQ=255)

                                                               gTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGa  <  1:476276/69‑1 (MQ=255)

                                                               gTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGa  <  1:1730305/69‑1 (MQ=255)

                                                               gTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGa  <  1:195251/69‑1 (MQ=255)

                                                               gTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGa  <  1:2489702/69‑1 (MQ=255)

CGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGA  >  minE/1230920‑1231051

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