Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • A problem of using Samtools to make VCF file

    Hi!

    I used the samtools commands and 1000 genome data(phase 3, low-coverage sequenced data) to generate the VCF file for both snp and indel, and below are the two command lines to get the vcf file for 5 people. The final "5genomes.var.flt.vcf" file has around 163432 snp and indel sites. However, when I used the same way to generate the vcf file for 100 people, the final vcf file includes only 1465 snp and indel, which is too few. I checked the last locus for both files, it looks like both of them processed the whole chromosome 20.

    Do you have any idea what's wrong with it?

    Thanks a lot. I look forward to your reply.

    ————— Here are the commands I used to generate VCF file——————

    ../samtools-0.1.19/samtools mpileup -uf ../human_g1k_v37.fasta HG00096.chrom20.ILLUMINA.bwa.GBR.low_coverage.20120522.bam HG 00268.chrom20.ILLUMINA.bwa.FIN.low_coverage.20130415.bam HG00419.chrom20.ILLUMINA.bwa.CHS.low_coverage.20130415.bam HG00759. chrom20.ILLUMINA.bwa.CDX.low_coverage.20130415.bam HG01112.chrom2 0.ILLUMINA.bwa.CLM.low_coverage.20120522.bam | ../samtools-0.1.19/bcftools/bcftools view -bvc - > 5genomes.var.raw.bcf


    ../samtools-0.1.19/bcftools/bcftools view 5genomes.var.raw.bcf | ../samtools-0.1.19/bcftools/vcfutils.pl varFilter -D100 > 5genomes.var.flt.vcf

Latest Articles

Collapse

  • seqadmin
    Techniques and Challenges in Conservation Genomics
    by seqadmin



    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

    Avian Conservation
    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
    03-08-2024, 10:41 AM
  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-14-2024, 06:13 AM
0 responses
32 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-08-2024, 08:03 AM
0 responses
71 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-07-2024, 08:13 AM
0 responses
80 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-06-2024, 09:51 AM
0 responses
68 views
0 likes
Last Post seqadmin  
Working...
X