BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,523,255	(T)_6→7	coding (24/2190 nt)	yfeA ←	predicted diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,523,249	1	.	T	100.0%	133.2 / NA	41	N10K (AAT→AAA)	yfeA	predicted diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base T (30/11); total (30/11)

GGATAGCTGAATAAGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATA‑TTTTTTATCAGGTTATGCTCCACAAACATGAGTAC  >  W3110S.gb/2523195‑2523284

ggATAGCTGAATAAGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATA‑ttttttt                              <  1:2139635/62‑2 (MQ=255)

           tAAGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATg                    >  1:887974/1‑62 (MQ=255)

           tAAGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATg                    >  1:4438272/1‑62 (MQ=255)

             aGTACCACGGTGAGCGTAAACTCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCt                  <  1:1327679/62‑1 (MQ=255)

             aGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCt                  <  1:2945899/62‑1 (MQ=255)

             aGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCt                  <  1:4086305/62‑1 (MQ=255)

             aGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCt                  <  1:3826427/62‑1 (MQ=255)

             aGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCt                  <  1:2919210/62‑1 (MQ=255)

             aGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCt                  <  1:258510/62‑1 (MQ=255)

             aGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCt                  <  1:197057/62‑1 (MQ=255)

                aCCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCca               <  1:4375576/62‑1 (MQ=255)

                aCCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCca               <  1:3741994/62‑1 (MQ=255)

                aCCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCca               <  1:2021681/62‑1 (MQ=255)

                aCCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCca               <  1:2812966/62‑1 (MQ=255)

                      gTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAAc          >  1:2284783/1‑61 (MQ=255)

                      gTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAAc          >  1:1782631/1‑61 (MQ=255)

                      gTGAGCGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAAc          >  1:3067956/1‑61 (MQ=255)

                      gTGAGAGTAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAAc          >  1:2307557/1‑61 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTcc                >  1:4353486/1‑48 (MQ=39)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:931000/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:1430040/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:608744/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:4473439/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:4348020/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:2356333/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:1941273/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:2383947/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACaaa           >  1:1097722/1‑53 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:3666753/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:2473793/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:2341611/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:4438187/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:2504409/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:618147/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:640817/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:693088/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTac  >  1:2654273/1‑62 (MQ=255)

                             tAAACGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTa   >  1:3859444/1‑61 (MQ=255)

                                aCGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTa   >  1:3666360/1‑58 (MQ=255)

                                aCGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTa   >  1:148539/1‑58 (MQ=255)

                                aCGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTCCACAAACATGAGTa   >  1:2671780/1‑58 (MQ=255)

                                aCGCTAGTGTGAATATCTTGATATTTTTTTATCAGGTTATGCTACACAAACATGAGTa   >  1:4422304/1‑58 (MQ=255)

GGATAGCTGAATAAGTACCACGGTGAGCGTAAACGCTAGTGTGAATATCTTGATA‑TTTTTTATCAGGTTATGCTCCACAAACATGAGTAC  >  W3110S.gb/2523195‑2523284

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