BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	861,111	G→A	L640L (CTG→TTG)	ybiW ←	predicted pyruvate formate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	861,111	0	G	A	100.0%	145.7 / NA	44	L640L (CTG→TTG)	ybiW	predicted pyruvate formate lyase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (20/24); total (20/24)

CGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCG  >  W3110S.gb/861052‑861149

cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAAtt                                      <  1:3957189/62‑1 (MQ=255)

cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAAtt                                      <  1:80179/62‑1 (MQ=255)

cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAAtt                                      <  1:2303406/62‑1 (MQ=255)

cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAAtt                                      <  1:674902/62‑1 (MQ=255)

 gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:852049/62‑1 (MQ=255)

 gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:3756777/62‑1 (MQ=255)

 gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:2335802/62‑1 (MQ=255)

 gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:4322763/62‑1 (MQ=255)

 gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:475099/62‑1 (MQ=255)

 gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:2753878/62‑1 (MQ=255)

 gcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:359441/62‑1 (MQ=255)

 gagTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:2535583/60‑1 (MQ=255)

           tGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAg                           >  1:282178/1‑62 (MQ=255)

            gccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:1565741/51‑1 (MQ=255)

            gccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:412324/51‑1 (MQ=255)

             ccgccAACCGGTCCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:1288680/1‑62 (MQ=255)

             ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:4560966/1‑62 (MQ=255)

             ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:4554217/1‑62 (MQ=255)

             ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:3565820/1‑62 (MQ=255)

             ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:2444616/1‑62 (MQ=255)

             ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:3723058/1‑62 (MQ=255)

             ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:4264603/1‑62 (MQ=255)

             ccgccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:4026921/1‑62 (MQ=255)

                ccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                         >  1:2939279/1‑59 (MQ=255)

                ccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCt                       >  1:1928688/1‑61 (MQ=255)

                ccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGACt                       >  1:1861562/1‑61 (MQ=255)

                  aaCCGGAGCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtcgtc                                  >  1:994379/1‑48 (MQ=38)

                  aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtcgtc                                  >  1:3933102/1‑48 (MQ=255)

                  aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATtcgtc                                  >  1:762320/1‑48 (MQ=255)

                         ccACGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:4603427/38‑1 (MQ=255)

                           aCGACCGTAGCGCGGATTATGGTACTGTTTCAATtc                                     <  1:1568834/36‑1 (MQ=255)

                            cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGa            >  1:4283294/1‑60 (MQ=255)

                            cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGa            >  1:361406/1‑60 (MQ=255)

                            cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGa            >  1:2957753/1‑60 (MQ=255)

                            cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGa            >  1:2683219/1‑60 (MQ=255)

                            cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGa            >  1:569464/1‑60 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:1473909/62‑1 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:4104952/62‑1 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:2213167/62‑1 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:2596491/62‑1 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:4594265/62‑1 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:3245686/62‑1 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:4505267/62‑1 (MQ=255)

                                    gcgcGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg  <  1:3884201/62‑1 (MQ=255)

CGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCG  >  W3110S.gb/861052‑861149

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