Question: Error in SAMTools mpileup
gravatar for julia92796
10 months ago by
julia927960 wrote:

I have a BAM file containing an alignment of the neanderthal y chromosome to the human y chromosome. I'm trying to extract the consensus sequence of the neanderthal reads so that I can isolate genes. I found this command online:

samtools mpileup -uf ref.fa aln.bam | bcftools view -cg - | vcf2fq > cns.fq

and am trying to use it to output the consensus sequence. However, I keep getting an error:

Error: Could not parse --min-ac g Use of uninitialized value $l in numeric lt (<) at line 566.
Use of uninitialized value $l in numeric lt (<) at line 566.
<mpileup> Set max per-file depth to 8000

Does anyone know how to fix this? If not, is there an alternate method to get a consensus sequence from reads?

ADD COMMENTlink written 10 months ago by julia927960

I also encounter the similar problem, and haven't resolved this yet.

Is there anyone who knows what's wrong with this error?

Thanks a lot. : )

ADD REPLYlink written 9 months ago by galanhsnu0
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: 966 users visited in the last hour