Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: How to output all variant alleles at a given site, regardless of quality
0
gravatar for correlationmatrix
16 months ago by
correlationmatrix • 20
correlationmatrix • 20 wrote:

How can one use either samtools mpileup or any other tool to, for a given position, output all alternate alleles regardless of quality? Just "samtools mpileup" (and the parameters I've tried so far) gives <*> for a position for which I observe an alternate base in IGV. (The purpose of doing this is for quality control)
mpileup snp variant • 500 views
ADD COMMENT • link •
modified 16 months ago by Brian Bushnell ♦ 15k • written 16 months ago by correlationmatrix • 20
2
gravatar for Brian Bushnell
16 months ago by
Brian Bushnell ♦ 15k
Walnut Creek, USA
Brian Bushnell ♦ 15k wrote:

You can do this with BBMap's CallVariants tool:

callvariants.sh in=mapped.bam ref=ref.fasta ploidy=2 vcf=all.vcf clearfilters

"clearfilters" sets all filters to zero, so no variant will fail. Your vcf file might be pretty big, though.
ADD COMMENT • link written 16 months ago by Brian Bushnell ♦ 15k

Works! Although in this case I was only interested in looking up specific regions, but subsetting the BAM file using samtools first solves the problem. However, I got some Java errors when running callvariants.sh: "Invalid maximum heap size: -Xmx101593m. The specified size exceeds the maximum representable size.". Lowering the memory parameters to "-Xmx2048m -Xms2048m" still gave errors about "java.lang.OutOfMemoryError: Java heap space", but the program finished anyway and the output VCF seemed to look fine.

ADD REPLY • link written 16 months ago by correlationmatrix • 20

Oh... I'm guessing that you are using an older 32-bit version of Java? I suggest you upgrade to Java 8 64-bit version.

ADD REPLY • link written 16 months ago by Brian Bushnell ♦ 15k

You're right, that was the issue!

ADD REPLY • link written 16 months ago by correlationmatrix • 20
Please log in to add an answer.

Similar posts • Search »

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