BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,007,247	T→C	pseudogene (75/393 nt)	yhiK →	ECK3474:JW5674:b3489; hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,007,247	0	T	C	100.0%	212.6 / NA	62	pseudogene (75/393 nt)	yhiK	ECK3474:JW5674:b3489; hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (38/24); total (38/24)

ATGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACATAAAGAACACATGGATGATCCTCTATTCCATCTAATCC  >  W3110S.gb/4007194‑4007284

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAg                                   >  1:1143976/1‑58 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:2888794/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:2771879/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:2717907/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:2714497/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:2494448/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:2224895/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:2019747/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAaca                               >  1:1915523/1‑62 (MQ=255)

aTGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAac                                >  1:159530/1‑61 (MQ=255)

       aTCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGgat                        <  1:112750/62‑1 (MQ=255)

       aTCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGgat                        <  1:657777/62‑1 (MQ=255)

       aTCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGgat                        <  1:1946349/62‑1 (MQ=255)

       aTCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGgat                        <  1:519/62‑1 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGga                         >  1:384961/1‑46 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTcc           >  1:995349/1‑60 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTcc           >  1:768281/1‑60 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:442168/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2881103/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2812476/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2745890/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:70149/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:692973/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2896573/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2641564/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2496558/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:359485/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:1235691/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:1005561/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2385971/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:123614/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:381759/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:2056995/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:1387072/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:535399/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:1947821/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:1490766/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:1492697/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:433124/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:184489/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:438456/1‑61 (MQ=255)

                      aGGAATTAGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCa          >  1:1596518/1‑61 (MQ=255)

                           ttATCCAATATTTTTAATAAAAAACACAAAGAACACATGgat                        <  1:2193816/42‑1 (MQ=38)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2396717/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2848076/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2817115/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1058122/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1086255/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1173011/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1204234/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1389540/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1411129/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1517576/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1747233/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1876382/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2036282/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2175175/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2494479/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2811697/62‑1 (MQ=255)

                             aGCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2841806/62‑1 (MQ=255)

                              gCCAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:2674651/61‑1 (MQ=255)

                               ccAATATTTTTAATAAAAAACACAAAGAACACATGGATGATCCTCTATTCCATCTAATcc  <  1:1786944/60‑1 (MQ=255)

ATGAAAAATCCTACTGATCCTGAGGAATTAGCCAATATTTTTAATAAAAAACATAAAGAACACATGGATGATCCTCTATTCCATCTAATCC  >  W3110S.gb/4007194‑4007284

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