SNP分析

老师,我想问一下您,做群体遗传时,分析SNP 变异率时,可以用zcat clean.vcf.gz|grep -v "#"|wc -l得到总的SNP位点数除以参考基因组的大小来计算么?

另外对于vcftools --gzvcf all.varFilter.vcf.gz --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这里的参数调整以后,直接影响结果的大小比较,比较迷茫,如何根据自己的数据去调整合适的参数大小

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

不同的过滤条件得到的SNP数量肯定不一样的,

你那个计算方法只是得到基因组上SNP的密度,你再看看文献吧,看看那个指标比较适合你想表达的意思;

请先 登录 后评论
  • 2 关注
  • 0 收藏,1025 浏览
  • 叶雪儿 提出于 2023-12-27 22:07

相似问题