Samtools Mpileup Error In Making The Consensus Sequences
Dear All
I am trying to make consensus sequences using command
samtools mpileup -uf ref.fa aln.bam | bcftools view -cg - | vcfutils.pl vcf2fq > cns.fq.
But i am getting this error:
[bcf_sync] incorrect number of fields (4 != 5) at 15:16294867
[afs] 0:49470.732 1:50.052 2:33.216
Can some one help me to sort out this problem. What exactly is going wrong.
Regards
1 answer
Log in to answer this question.
More posts like this
-
How to make consensus sequence from aligned reads ?
written by Info.shi •
-
get consensus sequence from bam file with insertions
written by Vvedenskiia •
-
How to do consensus calling from bam files
written by AW
-
strange output when trying to get consensus from a sorted bam file
written by lait
-
Error in SAMTools mpileup
written by julia92796 •
-
Generating wrong consensus sequence using samtools + bcftools
written by zhuhenan0514 •
-
Samtools Mpileup Defined Regions
written by Vova Naumov
-
Pileup Haploid Base Calling, Creating Consensus Sequence
written by Lee Katz
-
Consensus With Mpileup
written by Dunand •
-
To Avoid The 'Argument "F" Isn'T Numeric In Numeric
written by Zhshqzyc •
