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
        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, Yesterday, 06:55 AM
      0 responses
      12 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-30-2024, 03:16 PM
      0 responses
      24 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-29-2024, 01:32 PM
      0 responses
      29 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 05-24-2024, 07:15 AM
      0 responses
      215 views
      0 likes
      Last Post seqadmin  
      Working...
      X