BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,068,568	T→A	G22G (GGT→GGA)	ymdF →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,068,568	0	T	A	100.0%	186.6 / NA	56	G22G (GGT→GGA)	ymdF	conserved hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (34/22); total (34/22)

CAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGA  >  W3110S.gb/1068537‑1068627

cAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAAttt                             <  1:1013323/64‑1 (MQ=255)

 agaagaCCGCGAAAGCGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:2016920/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:17741/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:585503/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:590136/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:583417/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:2111924/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:616912/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:664310/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:1934394/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:498426/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:2255823/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:1806152/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:2303231/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:776839/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:2426427/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:82444/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:2455281/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:276256/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:473523/64‑1 (MQ=255)

 agaagaCCGCGAAAGAGCATCAGAAACAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:1651586/64‑1 (MQ=255)

  gaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTc                            <  1:1303377/63‑1 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGTGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:455609/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTg   >  1:1871543/1‑64 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:601575/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:671344/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:487235/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:773314/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:789562/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:812811/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:480468/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:845245/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:857849/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:975619/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2306058/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2018920/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1077368/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1103268/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1174516/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1188058/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1306683/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1318497/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1391201/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1504575/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1923742/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:1973201/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:421350/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:212736/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2156586/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2219659/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2234835/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2241704/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2278219/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:2450666/1‑65 (MQ=255)

                          gCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:34013/1‑65 (MQ=255)

                          gCAGGAAAAAAAGCTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa  >  1:285913/1‑65 (MQ=255)

CAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGA  >  W3110S.gb/1068537‑1068627

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