Re: [Vcftools-help] AC=0

[Warren W. K. <wkr...@gm...>]

Worst-case scenario, this would do it:

 perl -ane 'if(F[7] != m/^AC=0;/){ print $_}'  > output.vcf < input.vcf

This'll likely be slow though and I have not tested it.


On Mon, May 26, 2014 at 1:10 PM, Joanna Kelley <joa...@gm...>wrote:

> Hi all,
>
> I have a vcf with non-variable positions, a small example is below. I'd
> like to filter out the non-variable sites. I've tried both the --mac 1 and
> the --non-ref-ac 1 flags, but neither removes those sites. Any advice on
> how to remove the non-variable sites?
>
> Thanks,
> Joanna
>
> #CHROM    POS    ID      REF    ALT    QUAL    FILTER    INFO    FORMAT
> Victoria    9732    ....
> chr1    958777    chr1_958777    C    .    .    .    AC=0;AN=24      GT
> 0/0    0/0    0/0
> chr1    958778    chr1_958778    T    .    .    .    AC=0;AN=24      GT
> 0/0    0/0    0/0
> chr1    958779    chr1_958779    T    .    .    .    AC=0;AN=24      GT
> 0/0    0/0    0/0
> chr1    958780    chr1_958780    A    G    .    .    AC=2;AN=24      GT
> 0/0    0/0    1/1
>
>
> ------------------------------------------------------------------------------
> The best possible search technologies are now affordable for all companies.
> Download your FREE open source Enterprise Search Engine today!
> Our experts will assist you in its installation for $59/mo, no commitment.
> Test it for FREE on our Cloud platform anytime!
>
> http://pubads.g.doubleclick.net/gampad/clk?id=145328191&iu=/4140/ostg.clktrk
> _______________________________________________
> Vcftools-help mailing list
> Vcf...@li...
> https://lists.sourceforge.net/lists/listinfo/vcftools-help
>
>