Question: Consensus from mpileup error
gravatar for mks002
2.6 years ago by
Bangalore. India
mks002160 wrote:

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 seeing the consensus output, the indels are not incorporated i.e for the insertion the bases should be added at that particular position and deleted in the other case. Only the SNPs are showing up properly in consensus.

So why indels are coming or what should be done to get it correctly done.

   samtools-0.1.18/samtools mpileup -E -uf reference.fna sorted.bam | samtools-0.1.18/bcftools/bcftools view -cg - | samtools-0.1.18/bcftools/ vcf2fq > cns.fastq

seqtk seq -a cns.fastq > cns.fasta

One more query is that the start 1-3 position of the consensus is also changed, which i checked in the mpileup is same as reference without and SNPs or indels. why is this issue coming.

What would be the possible way to do this. Awaiting for kind replies .:)

sequencing snp alignment genome • 943 views
ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by mks002160
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 770 users visited in the last hour