BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,055,632	C→T	A9A (GCG→GCA)	yrhB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,055,632	0	C	T	100.0%	115.5 / NA	36	A9A (GCG→GCA)	yrhB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (10/26); total (10/26)

CATGTAGAGTAATGACGACCGGGAGATCTGCATCATCAAGGTAATGGTTCGCTTTCGCGAATGCGTCGTGATAAGTAATCATAATTATAATAAATATTATTGTTGAGTG  >  W3110S.gb/4055577‑4055685

cATGTAGAGTAATGACGACCGGGAGATCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAAt                                                 <  1:2125788/62‑1 (MQ=255)

       agTAATGACGACCGGGAGATCTGCATCATCAAGGTAATGGTTCGCTTTtgcg                                                    >  1:366878/1‑52 (MQ=255)

         tAATGACGACCGGGAGATCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAAt                                                 <  1:2863703/53‑1 (MQ=255)

              acgacCGGGAGATCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGt                                   <  1:929382/62‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2081641/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2593519/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2950604/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2372995/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2350409/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:3023065/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2178625/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2139446/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:1875028/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:186859/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:843144/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:893283/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:1752931/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:1745954/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:1742119/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:163538/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:160374/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:1484770/61‑1 (MQ=255)

                       tgaTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGGCGTGATAAGTAATCATAAt                          <  1:2578113/61‑1 (MQ=255)

                        gaTCTGCATCATCAAGGTAATGGTTCTCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2808197/61‑1 (MQ=255)

                        gaTCTGCATCATCAAGGGAATGGATCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:152088/61‑1 (MQ=255)

                         aTCTGCATCATCAAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAAt                          <  1:2304430/60‑1 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATta                        >  1:996404/1‑53 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATaataaat                >  1:2396071/1‑61 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATaataaat                >  1:136452/1‑61 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATaataaat                >  1:2245778/1‑61 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATaataaat                >  1:394067/1‑61 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATaataaat                >  1:812908/1‑61 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATaataaat                >  1:1823060/1‑61 (MQ=255)

                                  atcaAGGTAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATaataaa                 >  1:179262/1‑60 (MQ=255)

                                         tAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATAATAAATATTAttg         >  1:484623/1‑61 (MQ=255)

                                         tAATGGTTCGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATAATAAATATTAttg         >  1:963057/1‑61 (MQ=255)

                                                 cGCTTTTGCGAATGCGTCGTGATAAGTAATCATAATTATAATAAATATTATTGTTGAgtg  <  1:158895/60‑1 (MQ=255)

CATGTAGAGTAATGACGACCGGGAGATCTGCATCATCAAGGTAATGGTTCGCTTTCGCGAATGCGTCGTGATAAGTAATCATAATTATAATAAATATTATTGTTGAGTG  >  W3110S.gb/4055577‑4055685

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