Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • mpileup only calling snps for chr1 not entire genome

    Hello,

    I am having a problem when trying to call SNPS from two alignment files, and this also happened when I tried calling SNPs from only 1 alignment file. My problem is I am only getting SNPs for Chr01, and not for the entire genome (total Chr10).

    Here is the code I am using:
    samtools mpileup -B -Q 20 -P Illumina -C50 -u -f genome.fa wt_align_sorted.bam mut_align_sorted.bam | bcftools view -vcg - >var.vcf

    I also have a second question in terms of the genome fasta to use. I have run the samtools faidx command on my genome and get 4 columns with 10 rows. I have tried using that output ".fai" file in the mpileup command but it does not work. Do I use the faidx file I mentioned above or the fasta file I used for mapping? When I use the fasta file I used for mapping without indexing I get SNPs but only for Chr01.

    I am using Samtools v 1.3 and bcftools v 1.3

    Thank you for your help

  • #2
    issue resolved

    Issue was resolved for this problem. Turns out the bcftools version 1.3 uses -call instead of the older version that uses -view. And the reason it crashed after chr01 was because I did not allot enough space for it to finish.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Best Practices for Single-Cell Sequencing Analysis
      by seqadmin



      While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
      Yesterday, 07:15 AM
    • 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

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 08:18 AM
    0 responses
    13 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 08:04 AM
    0 responses
    12 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 06-03-2024, 06:55 AM
    0 responses
    13 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-30-2024, 03:16 PM
    0 responses
    27 views
    0 likes
    Last Post seqadmin  
    Working...
    X