BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	810,589	T→G	54.6%	M275R (ATG→AGG)	bioB →	biotin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	810,589	0	T	G	54.6%	11.2 / 21.0	22	M275R (ATG→AGG)	bioB	biotin synthase
Reads supporting (aligned to +/- strand): ref base T (4/6); new base G (6/6); total (10/12)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

CATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGC  >  W3110S.gb/810521‑810645

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:2507148/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:980632/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:556339/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:3157612/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:2354198/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:2206568/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:1941118/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:1575569/1‑71 (MQ=255)

cATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAggt                                                        >  1:1381615/1‑71 (MQ=255)

    atgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                    <  1:1087997/71‑1 (MQ=255)

    atgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                    <  1:2131115/71‑1 (MQ=255)

    atgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                    <  1:3258972/71‑1 (MQ=255)

     tgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                    <  1:2862244/70‑1 (MQ=255)

     tgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                    <  1:1778494/70‑1 (MQ=255)

             tctTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAg                                            >  1:12377/1‑70 (MQ=255)

                             cTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                    <  1:3070371/46‑1 (MQ=255)

                                     cgcgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCg                                 <  1:2997707/57‑1 (MQ=255)

                                     cgcgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCg                                 <  1:559374/57‑1 (MQ=255)

                                     cgcgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCg                                 <  1:643397/57‑1 (MQ=255)

                                     cgcgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCg                                 <  1:740654/57‑1 (MQ=255)

                                     cgcgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCg                                 <  1:954214/57‑1 (MQ=255)

                                       cgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGa                                >  1:1530009/1‑56 (MQ=255)

                                       cgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAt                               >  1:2603567/1‑57 (MQ=255)

                                       cgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAt                               >  1:2754254/1‑57 (MQ=255)

                                       cgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAt                               >  1:2773106/1‑57 (MQ=255)

                                       cgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAt                               >  1:1073553/1‑57 (MQ=255)

                                       cgAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAt                               >  1:2365234/1‑57 (MQ=255)

                                                      acAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGc  >  1:1606425/1‑71 (MQ=255)

                                                      acAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGc  >  1:3043572/1‑71 (MQ=255)

                                                      acAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGc  >  1:1391224/1‑71 (MQ=255)

                                                      acAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCg                                 <  1:2278417/40‑1 (MQ=255)

CATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGC  >  W3110S.gb/810521‑810645

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