Entering edit mode
8 weeks ago

I have a vcf file with samples from one species called against the reference genome of an other species. I am trying to remove sites any sites fixed among my samples for non reference alleles.

I have constructed this code below for the 11 samples in my vcf file. I am not sure if this is the correct logic for what I am attempting to filter from my vcf file (I included an example for just three samples for simplicity of the post!). I have not been able to find any direction on other posts and would very much appreciate any suggestions or guidance. Thank you!!

zcat |\
SnpSift filter \

"isHom( GEN[0] ) & isRef( GEN[0] ) | isHet( GEN[0] ) && \ isHom( GEN[1] ) & isRef( GEN[1] ) | isHet( GEN[1] ) && \ isHom( GEN[2] ) & isRef( GEN[2] ) | isHet( GEN[2] ) && \"

alleles filter SNPsift nonrefence • 125 views

Login before adding your answer.

Traffic: 2692 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6