Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Compare VCF Files

    Can anyone tell me how to compare the vcf files to find the differences of two or more? Use Existing tools or souce code?

  • #2
    Originally posted by lyz1030 View Post
    Can anyone tell me how to compare the vcf files to find the differences of two or more? Use Existing tools or souce code?

    Comment


    • #3
      Common snps

      which vcftool i can use to find snps that are common between two vcf files.
      If there is any other software please recommend
      Regards

      Comment


      • #4
        Originally posted by huma Asif View Post
        which vcftool i can use to find snps that are common between two vcf files.
        If there is any other software please recommend
        Regards
        vcf-isec

        (padding to 10 characters to post)

        Comment


        • #5
          vcf-isec

          i used vcf-isec to filter snps found in one file but not in other using
          ./vcf-isec -c -f /home/fisa/orignal_vcf/A.vcf.gz /home/fisa/database_vcf/1kgenome.vcf.gz >/home/huma/lifescope_vcf/A_filtered_1Kgenome.vcf

          what changes i need in command to get common

          Comment


          • #6
            You can drop the -c argument.

            Comment


            • #7
              thank you

              let me try

              Comment


              • #8
                There's also VAAST if you want to look at more complicated interactions (e.g. SNPs in group 1, but not group 2, unless they're in group 3).

                Comment


                • #9
                  vcf-isec error

                  i used vcf-isec to find common snps between my two files using command
                  ./vcf-isec -f /home/asi/orignal_vcf/A.vcf.gz /home/asi/database_vcf/1Kgenome.vcf.gz >/home/asi/lifescope_vcf/A_common_1Kgenome.vcf
                  i got this error
                  Use of uninitialized value $fltr in string ne at ./vcf-isec line 337, <$__ANONIO__> line 74872.
                  Use of uninitialized value $fltr in string ne at ./vcf-isec line 337, <$__ANONIO__> line 74872.

                  Any help would be highly appreciated

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Essential Discoveries and Tools in Epitranscriptomics
                    by seqadmin




                    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
                    04-22-2024, 07:01 AM
                  • seqadmin
                    Current Approaches to Protein Sequencing
                    by seqadmin


                    Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                    04-04-2024, 04:25 PM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, Today, 08:06 AM
                  0 responses
                  13 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 04-30-2024, 12:17 PM
                  0 responses
                  15 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 04-29-2024, 10:49 AM
                  0 responses
                  21 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 04-25-2024, 11:49 AM
                  0 responses
                  27 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X