Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • htetre
    Member
    • Jul 2013
    • 28

    #1

    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
  • htetre
    Member
    • Jul 2013
    • 28

    #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

    • SEQadmin2
      Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
      by SEQadmin2



      Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

      Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
      ...
      10-14-2026, 09:55 AM
    • SEQadmin2
      New Genomics Technologies Take Aim at Long-Standing Limits
      by SEQadmin2


      Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

      We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
      ...
      09-28-2026, 10:25 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, Yesterday, 11:13 AM
    0 responses
    16 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, Yesterday, 10:48 AM
    0 responses
    16 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-15-2026, 10:52 AM
    0 responses
    35 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-14-2026, 05:39 AM
    0 responses
    42 views
    0 reactions
    Last Post SEQadmin2  
    Working...