Hello,
When I use the NGS Picard tool "SAM/BAM alignment summary metrics" with a BAM file produced with Tophat in Galaxy I get the following message:
INFO:root:## executing java -Xmx4g -jar /galaxy/home/g2main/galaxy_main/tool-data/shared/jars/CollectAlignmentSummaryMetrics.jar VALIDATION_STRINGENCY=LENIENT ASSUME_SORTED=true ADAPTER_SEQUENCE= IS_BISULFITE_SEQUENCED=false MAX_INSERT_SIZE=1000 OUTPUT=/galaxy/main_pool/pool5/tmp/job_working_directory/2368774/dataset_2672529_files/CollectAlignmentSummaryMetrics.metrics.txt R=/galaxy/main_pool/pool5/tmp/job_working_directory/2368774/dataset_2672529_files/mm9.fa_fake.fasta TMP_DIR=/tmp INPUT=/galaxy/main_database/files/002/666/dataset_2666199.dat returned status 1 and stderr: [Mon Jul 11 03:27:56 EDT 2011] net.sf.picard.analysis.CollectAlignmentSummaryMetrics MAX_INSERT_SIZE=1000 ADAPTER_SEQUENCE=[AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT, AGATCGGAAGAGCTCGTATGCCGTCTTCTGCTTG, AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT, AGATCGGAAGAGCGGTTCAGCAGGAATGCCGAGACCGATCTCGTATGCCGTCTTCTGCTTG, AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT, AGATCGGAAGAGCACACGTCTGAACTCCAGTCACNNNNNNNNATCTCGTATGCCGTCTTCTGCTTG, IS_BISULFITE_SEQUENCED=false] INPUT=/galaxy/main_database/files/002/666/dataset_2666199.dat OUTPUT=/galaxy/main_pool/pool5/tmp/job_working_directory/2368774/dataset_2672529_files/CollectAlignmentSummaryMetrics.metrics.txt REFERENCE_SEQUENCE=/galaxy/main_pool/pool5/tmp/job_working_directory/2368774/dataset_2672529_files/mm9.fa_fake.fasta ASSUME_SORTED=true TMP_DIR=/tmp VALIDATION_STRINGENCY=LENIENT IS_BISULFITE_SEQUENCED=false STOP_AFTER=0 VERBOSITY=INFO QUIET=false COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false [Mon Jul 11 03:27:57 EDT 2011] net.sf.picard.analysis.CollectAlignmentSummaryMetrics done. Runtime.totalMemory()=507379712 Exception in thread "main" java.lang.IllegalArgumentException: No enum const class net.sf.samtools.SAMFileHeader$SortOrder.sorted at java.lang.Enum.valueOf(Enum.java:196) at net.sf.samtools.SAMFileHeader$SortOrder.valueOf(SAMFileHeader.java:58) at net.sf.samtools.SAMFileHeader.getSortOrder(SAMFileHeader.java:239) at net.sf.picard.analysis.SinglePassSamProgram.makeItSo(SinglePassSamProgram.java:85) at net.sf.picard.analysis.SinglePassSamProgram.doWork(SinglePassSamProgram.java:54) at net.sf.picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:157) at net.sf.picard.cmdline.CommandLineProgram.instanceMainWithExit(CommandLineProgram.java:117) at net.sf.picard.analysis.CollectAlignmentSummaryMetrics.main(CollectAlignmentSummaryMetrics.java:106)
If I use the NGS: Picard BAM Index Statistics, the script is running without any failure so I think the BAM format should be okay. Any help would be gratefully received.
Thanks.
Andreas Huge