On 7/11/14, 9:19 AM, Kanterakis, Efstathios wrote:

Hello,

Even though the usage message states "This option can be set to 'null' to clear the default value", when I run Picard (v1.115) CollectWgsMetrics.jar with MQ=null, I get

java.lang.IllegalArgumentException: Can not set int field picard.analysis.CollectWgsMetrics.MINIMUM_MAPPING_QUALITY to null value

same for CAP=null

I tried setting MQ=-1 which worked, but CAP=-1 gave me an java.lang.ArrayIndexOutOfBoundsException: -1

In the end I ran with MQ=-1, CAP=100000


Hi Stathis,

This is a bug in the code that generates the usage message. MQ is required, as you discovered.

And a feature request, would it be possible to report the number of mapped and unmapped bases so that one could calculate the number of bases remaining after filters?

Does CollectAlignmentSummaryMetrics not provide the information you need?

-Alec

Thank you very much.

Stathis



------------------------------------------------------------------------------
Open source business process management suite built on Java and Eclipse
Turn processes into business applications with Bonita BPM Community Edition
Quickly connect people, data, and systems into organized workflows
Winner of BOSSIE, CODIE, OW2 and Gartner awards
http://p.sf.net/sfu/Bonitasoft


_______________________________________________
Samtools-help mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/samtools-help

------------------------------------------------------------------------------
Want fast and easy access to all the code in your enterprise? Index and
search up to 200,000 lines of code with a free copy of Black Duck®
Code Sight™ - the same software that powers the world's largest code
search on Ohloh, the Black Duck Open Hub! Try it now.
http://p.sf.net/sfu/bds
_______________________________________________
Samtools-help mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/samtools-help

Reply via email to