BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,234,526	T→C	T28A (ACG→GCG)	dsbB ←	oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,234,526	0	T	C	100.0%	113.1 / NA	31	T28A (ACG→GCG)	dsbB	oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (9/22); total (9/22)

AATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTG  >  W3110S.gb/1234473‑1234582

aaTACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCa                                                  <  1:4593752/62‑1 (MQ=255)

 aTACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCa                                                  <  1:1158831/61‑1 (MQ=255)

 aTACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCa                                                  <  1:3070744/61‑1 (MQ=255)

   aCAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGt                                                >  1:3618998/1‑61 (MQ=255)

       agCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCag                                           <  1:355692/62‑1 (MQ=255)

                ggTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTaaa                                  <  1:3387007/62‑1 (MQ=255)

                      aGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCa                             >  1:2228260/1‑61 (MQ=255)

                      aGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCa                             >  1:2388118/1‑61 (MQ=255)

                      aGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCa                             >  1:497150/1‑61 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGGGCAGTAAACGCCATCa                          <  1:2665236/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:2606357/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:1009365/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:2995358/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:284098/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:4106969/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:2177385/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:2125411/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:4443570/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:1882531/62‑1 (MQ=255)

                        tAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCa                          <  1:4538400/62‑1 (MQ=255)

                                   gCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCAcgcg               <  1:3659761/62‑1 (MQ=255)

                                      ggAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGccc            >  1:4451891/1‑62 (MQ=255)

                                      ggAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGccc            >  1:117021/1‑62 (MQ=255)

                                      ggAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGccc            >  1:4706859/1‑62 (MQ=255)

                                      ggAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGcc             >  1:4017775/1‑61 (MQ=255)

                                         aCCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCgg         >  1:2068436/1‑62 (MQ=255)

                                           cacaGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc       <  1:515992/62‑1 (MQ=255)

                                           cacaGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc       <  1:3107771/62‑1 (MQ=255)

                                           cacaGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc       <  1:4467801/62‑1 (MQ=255)

                                           cacaGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc       <  1:1164008/62‑1 (MQ=255)

                                           cacaGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc       <  1:2451910/62‑1 (MQ=255)

                                                cgccgcCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAAc                        >  1:444619/1‑40 (MQ=255)

                                                cgccgcCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTtgtg  >  1:2800425/1‑62 (MQ=255)

                                                cgccgcCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTtgtg  >  1:1595625/1‑62 (MQ=255)

AATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTG  >  W3110S.gb/1234473‑1234582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`