Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to view only the variants that are present in multiple VCF files?

    Suppose I have 10 .vcf files and I want to generate a .vcf file that only contains the common variants that are present in all of the 10 files. Can I use GATK CombineVariants followed by SelectVariants intersection to achieve this? Is there another way?

    Thank you

    "If you want to extract just the records in common between two VCFs, you would first run CombineVariants on the two files to generate a single VCF and then run SelectVariants to extract the common records with -select 'set == "Intersection"', as worked out in the detailed example in the documentation guide."


  • #2
    VCFtools: http://vcftools.sourceforge.net/perl....html#vcf-isec

    Comment


    • #3
      I tried vcf-isec and it did not seem to work.

      I was eventually able to install and setup vcftools, including sorting, indexing, and compressing the vcf files with tabix and bgzip.

      The vcf-isec command gave a warning that column names do not match (i.e. 1-Normal and 1-Tumor). The command ran, but the output vcf file was 28 bytes of unreadable characters. Each of the 11 input files are around 80kb.

      These are vcf files generated using MuTect (for comparing tumor to normal samples).

      I am not sure if vcf-isec will work with MuTect vcf files or if there is something I am doing wrong. Maybe I can process the files ahead of time to get them to work.

      The idea is that MuTect gives a list of somatic mutations in cancer samples by comparing to matched normal samples. What I am trying to do is take several MuTect vcf files, and see which variants are present in multiple vcf files or the 'intersection.'

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Latest Developments in Precision Medicine
        by seqadmin



        Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

        Somatic Genomics
        “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
        05-24-2024, 01:16 PM
      • seqadmin
        Recent Advances in Sequencing Analysis Tools
        by seqadmin


        The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
        05-06-2024, 07:48 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 05-24-2024, 07:15 AM
      0 responses
      13 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-23-2024, 10:28 AM
      0 responses
      17 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-23-2024, 07:35 AM
      0 responses
      20 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-22-2024, 02:06 PM
      0 responses
      10 views
      0 likes
      Last Post seqadmin  
      Working...
      X