Question: Consensus from mpileup error
1
gravatar for mks002
16 months ago by
mks002100
Bangalore. India
mks002100 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/vcfutils.pl 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 • 502 views
ADD COMMENTlink modified 16 months ago • written 16 months ago by mks002100
Please log in to add an answer.

Help
Access

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