SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
picard-tools "Unrecognized tag type" error jtackney Bioinformatics 2 05-19-2013 03:55 PM
samtools installation error "cannot find -lcurses" ihelbig Bioinformatics 5 02-01-2013 11:48 AM
Picard Tools "loaded 0 genes" ckrna Bioinformatics 3 02-08-2012 11:26 AM
Error "RG ID on SAMRecord not found in header" from Picard's MarkDuplicates.jar‏ cliff Bioinformatics 4 11-10-2011 03:27 AM
picard: "No enum const clas" error mrfox Bioinformatics 1 04-17-2011 01:36 PM

Reply
 
Thread Tools
Old 08-14-2013, 10:11 AM   #1
serenaliao
Member
 
Location: Pittsburgh

Join Date: Jan 2013
Posts: 22
Default Using Picard "CalculateHsMetrics" to find coverage stat

Hi all,

I have been trying to use CalculateHsMetrics to get coverage information. I have two questions:

1. When I run the code:
java -Xmx10g -jar /home/liao/picard-tools-1.96/CalculateHsMetrics.jar BAIT_INTERVALS=target.bed TARGET_INTERVALS=target.bed INPUT=/data2/Pregnent_1_time_6_month/SAM/V8_GCCAAT_L003_R1_001.fastq.gz_bismark_pe.sam OUTPUT=sample.hybrid.stats.txt
Everything went well, and it took 20 mins to finish. However after I added PER_TARGET_COVERAGE options, and reference genome(fa extension), after a few seconds there's nothing, and no longer runnning:
java -Xmx10g -jar /home/liao/picard-tools-1.96/CalculateHsMetrics.jar BAIT_INTERVALS=target.bed TARGET_INTERVALS=target.bed INPUT=/data2/Pregnent_1_time_6_month/SAM/V8_GCCAAT_L003_R1_001.fastq.gz_bismark_pe.sam OUTPUT=sample.hybrid.stats.txt PER_TARGET_COVERAGE=Picard_PER_TARGET_COVERAGE.txt R=/home/liao/mm9/genome.fa
Does anyone know how to solve it?

2. So the near bait reads(according to their code) use 250 +/- off baits regions as criteria. Is there any convenient way to adjust the number, say to 10K,20K...?


Thanks in advance!!
serenaliao is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:57 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO