BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	65,135	C→T	D216N (GAT→AAT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	65,135	0	C	T	89.2%	103.0 / 4.8	37	D216N (GAT→AAT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand): ref base C (4/0); new base T (19/14); total (23/14)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACG  >  W3110S.gb/65076‑65190

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCaa                                              <  1:482740/71‑1 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCaa                                              <  1:3306107/71‑1 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCaa                                              <  1:3686420/71‑1 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCaa                                              <  1:2768736/71‑1 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCaa                                              <  1:260373/71‑1 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCaa                                              <  1:504076/71‑1 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCa                                               >  1:3910740/1‑70 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCa                                               >  1:2662863/1‑70 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCa                                               >  1:534199/1‑70 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCa                                               >  1:1478861/1‑70 (MQ=255)

tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCa                                               >  1:3808578/1‑70 (MQ=255)

             cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAg                                 >  1:213444/1‑71 (MQ=255)

             cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAg                                 >  1:2447087/1‑71 (MQ=255)

             cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAg                                 >  1:658345/1‑71 (MQ=255)

             cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAg                                 >  1:1531535/1‑71 (MQ=255)

             cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAg                                 >  1:1156590/1‑71 (MQ=255)

             cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATAAGGATTGTAGTTGGCAAACCAGGCGTTGa                                  >  1:1341673/1‑70 (MQ=255)

                    cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGttttt                            >  1:2613197/1‑69 (MQ=255)

                    cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTc                           >  1:3650980/1‑70 (MQ=255)

                    cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTc                           >  1:4093651/1‑70 (MQ=255)

                    cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTc                           >  1:651226/1‑70 (MQ=255)

                       aCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       <  1:916121/71‑1 (MQ=255)

                        cTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCa                      <  1:2819688/71‑1 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:495094/1‑37 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:542835/1‑37 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:565458/1‑37 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:351012/1‑37 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:3487462/1‑37 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:3201035/1‑37 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:2176516/1‑37 (MQ=255)

                          gCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                      >  1:1889078/1‑37 (MQ=255)

                                                    tCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:854815/63‑1 (MQ=255)

                                                    tCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:348394/63‑1 (MQ=255)

                                                    tCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:4115071/63‑1 (MQ=255)

                                                    tCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:2537969/63‑1 (MQ=255)

                                                    tCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:1361172/63‑1 (MQ=255)

                                                    tCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:1247971/63‑1 (MQ=255)

TTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACG  >  W3110S.gb/65076‑65190

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