Question: Error in SAMTools mpileup
gravatar for julia92796
23 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 23 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 21 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: 1144 users visited in the last hour