Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SNP base calling for multiple samples

    My goal is to find SNPs from multiple Sanger sequences, which cover varied regions, not necessary the same fragments. When a SNP is found, I want to see the base (or hetero/homo-zygote) of any samples which cover that position. What kind of software should I use? I would prefer a command-line software.


    I tried to align sequences via bwasw and pileup all via samtools. However, when I pileuped all samples together, the output did NOT clearly distinguish between when a SNP location not covered in a sample and when a sample had hetero in that SNP location. My commands below.


    samtools mpileup -uf Sorbil.fasta a1.bam a2.bam a3.bam | bcftools/bcftools view -bvcg - > raw.bcf

    bcftools/bcftools view raw.bcf | vcfutils.pl varFilter -D100 > flt.vcf

  • #2
    I'm not sure what you mean by "clearly", but you should see a 1/1, 0/1 or 0/0 in each entry. That tells you if the SNP is predicted to be homozygous alternate, heterozygous, or homozygous reference. Unfortunately, what you really want is the DP4 for each sample separately, and doing a multi-vcf with samtools won't do that. You'll need the individual vcf file to see that.

    Comment


    • #3
      In my case, some sequences showed 1/0 for a SNP position, however, some of those sequences did not cover that position.

      I want to know whether a SNP base is a polymorphic, or wide type, or uncovered in a given sequence/file.

      Comment

      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
      34 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-08-2024, 08:03 AM
      0 responses
      72 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-07-2024, 08:13 AM
      0 responses
      82 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