代码不适用

vcftools --gzvcf LP.vcf --recode --recode-INFO-all --stdout     --maf 0.05  --max-missing 0.8  --minDP 2  --maxDP 1000      \

>     --minQ 30 --minGQ 0 --min-alleles 2  --max-alleles 2 --remove-indels |gzip - > clean.vcf.gz

在用这个代码过滤vcf文件时,一秒就出结果,说代码不适用,我的物种是4倍体。

请先 登录 后评论

2 个回答

马杰宇

attachments-2023-12-2lWXnyYH657425e4a3747.pngVCFtools - UNKNOWN

(C) Adam Auton and Anthony Marcketta 2009


Parameters as interpreted:

--gzvcf LP.vcf

--recode-INFO-all

--maf 0.05

--max-alleles 2

--maxDP 1e+03

--min-alleles 2

--minDP 2

--minGQ 0

--minQ 30

--max-missing 0.8

--recode

--remove-indels

--stdout


Using zlib version: 1.2.7

Warning: Expected at least 2 parts in FORMAT entry: ID=PL,Number=.,Type=Float,Description="Normalized, Phred-scaled likelihoods for AA,AB,BB genotypes where A=ref and B=alt; not applicable if site is not biallelic">

Warning: Expected at least 2 parts in FORMAT entry: ID=PL,Number=.,Type=Float,Description="Normalized, Phred-scaled likelihoods for AA,AB,BB genotypes where A=ref and B=alt; not applicable if site is not biallelic">

After filtering, kept 258 out of 258 Individuals

Outputting VCF file...

After filtering, kept 0 out of a possible 105046 Sites

No data left for analysis!

Run Time = 0.00 seconds

这是结果文件
请先 登录 后评论
omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

你截图一下你的VCF,

如果不是二倍体的,你需要用专门的软件处理;

请先 登录 后评论
  • 2 关注
  • 0 收藏,1072 浏览
  • 马杰宇 提出于 2023-12-07 21:44

相似问题