BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,231,014	C→T	24.3%	G581D (GGT→GAT)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,231,014	0	C	T	24.3%	65.4 / 19.8	37	G581D (GGT→GAT)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base C (8/20); new base T (6/3); total (14/23)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTA  >  minE/1230964‑1231071

aCGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                                                <  1:451443/62‑1 (MQ=255)

aCGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                                                <  1:186658/62‑1 (MQ=255)

aCGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                                                <  1:187325/62‑1 (MQ=255)

aCGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                                                <  1:560403/62‑1 (MQ=255)

aCGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                                                <  1:575478/62‑1 (MQ=255)

 cGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaa                                               <  1:486849/62‑1 (MQ=255)

 cGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaa                                               <  1:554433/62‑1 (MQ=255)

 cGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCaa                                               <  1:50003/62‑1 (MQ=255)

 cGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCa                                                <  1:344765/61‑1 (MQ=255)

    gCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACt                                            >  1:152831/1‑62 (MQ=255)

    gCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACt                                            >  1:77089/1‑62 (MQ=255)

    gCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACt                                            >  1:247137/1‑62 (MQ=255)

                 aTGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAag                               >  1:226769/1‑62 (MQ=255)

                 aTGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAag                               >  1:504102/1‑62 (MQ=255)

                 aTGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAag                               >  1:570019/1‑62 (MQ=255)

                 aTGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAag                               >  1:173105/1‑62 (MQ=255)

                 aTGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAag                               >  1:454451/1‑62 (MQ=255)

                   gTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAaga                              <  1:135233/61‑1 (MQ=255)

                   gTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAaga                              <  1:147809/61‑1 (MQ=255)

                    ttcTCGTATTATTGCGGATCCTGCGGGTTACCCCACTCTTCa                                                <  1:168823/40‑1 (MQ=39)

                                       ccTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCaa              <  1:315805/57‑1 (MQ=255)

                                       ccTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCaa              <  1:462135/57‑1 (MQ=255)

                                        cTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCg        <  1:88699/62‑1 (MQ=255)

                                          gCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgtt      <  1:398126/62‑1 (MQ=255)

                                          gCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgtt      <  1:565884/62‑1 (MQ=255)

                                          gCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgtt      <  1:102934/62‑1 (MQ=255)

                                              gTTATCCCTCTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttgtt   >  1:242225/1‑61 (MQ=255)

                                              gTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttgtt   >  1:515479/1‑61 (MQ=255)

                                              gTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttgtt   >  1:362602/1‑61 (MQ=255)

                                              gTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttgtt   >  1:345353/1‑61 (MQ=255)

                                              gTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttgtt   >  1:579123/1‑61 (MQ=255)

                                              gTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttgtt   >  1:201226/1‑61 (MQ=255)

                                              gTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCGTCGTTGTTa  <  1:301677/62‑1 (MQ=255)

                                              gTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTa  <  1:440738/62‑1 (MQ=255)

                                              gTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTa  <  1:426667/62‑1 (MQ=255)

                                              gTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTa  <  1:506969/62‑1 (MQ=255)

                                              gTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTa  <  1:68590/62‑1 (MQ=255)

ACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTA  >  minE/1230964‑1231071

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