Problem about Getting Consensus Sequence Of Mapping Against Reference Genome Using Samtools Mpileup

Hi,

I have a set of reads mapped against a reference genome and need take consensus sequence.So I used the Samtools-Bcftools-Vcfutils.pl pipeline ,I used the command:

samtools mpileup -uf <reference.fasta> <file.sorted.bam> | \
  bcftools view -cg - | \
  vcfutils.pl vcf2fq | \
  perl fastq2fasta.pl cns.fasta

In the vcf file,there actually exists indels, but in the final consensus sequence, there is no indels.

For example, in the vcf

Supercontig_1.1 2953301 .       T       TCGTGGCCT       178     .       INDEL;IS=2,0.014493;DP=138;VDB=2.191453e-01;AF1=1;AC1=2;DP4=0,0,75,55;MQ=20;FQ=-290     GT:PL:GQ        1/1:219,255,0:99

Someone have any suggestions?

bcftools samtools consensus

• 2,904 views • link

updated 4.4 years ago by Ram

• written 11.7 years ago by wangyang703092

0 answers

No answers yet.

More posts like this

How to differentiate Homozygous reference calls from No calls in VCF
written 9.6 years ago by second_exon

I'm using samtools multisample SNP calling to call SNPs from different samples. The command I'm using is, samtools mpileup -uf ref.fa *.bam | bcftools view …
Consensus from mpileup error
written 10.0 years ago by mks002

Dear all, I am trying to generate the consensus sequence for whole genome sequencing. The reference has 11 segments. The issue is when i am …
VCF file REF column only contains N
written 10.2 years ago by genie66 •

I am using samtools mpileup to call variant. Command used: samtools mpileup -u -f "/gpfs/home/g/o/Homo_sapiens.GRCh37.75.dna.primary_assembly.fa" "/gopisiva/Mason/Single_cell/SC_Stressed/*.bam" | bcftools view -bvcg - > SC_Stressed.raw.bcf" bcftools view …
SNPs only output from bcftools
written 11.9 years ago by biogirl

Hi there, I'm trying to call SNPs from an alignment (BWA) to a reference genome using bcftools (0.1.18). However, my output contains INDELs too. I'm …
Allele frequency of vcf file
written 12.1 years ago by juncheng

I know it is a stupid question, but I'm a beginner of snp calling and really confused. I run samtools and bcftools get vcf files …
Why Is The Variant File Produced By Samtools So Large? It Is About ~34G.
written 13.0 years ago by Jordan

Hi, I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file of size 62G. When I do the variant analysis by <a href='http://samtools.sourceforge.net/'>samtools</a>, it gives a huge bcf file …
What Better Way To Get Consensus Sequence Of Mapping Against Reference Genome Using Samtools Mpileu…
written 13.3 years ago by thiagomafra

Hi people, I have a set of reads mapped against a reference genome and need take consensus sequence using <a href='http://samtools.sourceforge.net/'>samtools</a> mpileup. I used the …
I Did A Samtools Pileup And Ended In A Vcf Format Like This
written 14.3 years ago by satish.nair •

Hi Everyone I was doing the <a href='http://samtools.sourceforge.net/'>samtools</a> mpileup piping with bcftools using this command ("samtools mpileup -uf hg19.fa SingleCell5.chrM.bam | bcftools view -cg -v …
How To Filter The Snps/Indel Results From Mpileup?
written 15.0 years ago by Haiping •

Hi. I trid to use mpileup to identified the SNP/Indels. The commands I used were: <pre><code>samtools fillmd –bAr sample.sorted.bam ref.fa > sample.sorted.baq.bam samtools mpileup -uf …
How To Filter The Snps/Indel Results From Mpileup?
written 15.0 years ago by User 7551 •

Hi. I trid to use mpileup to identified the SNP/Indels. The commands I used were: samtools fillmd –bAr sample.sorted.bam ref.fa > sample.sorted.baq.bam samtools mpileup -ugf …