This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Picard MergeSam error : Cannot use index file with textual SAM file

I tried to run a simple mergesam using picard on two bam files using the following command

java -jar /opt/picard/picard.jar MergeSamFiles I=file1.bam I=file2.bam O=merged_output.bam USE_THREADING=TRUE

This is giving me the following error

Exception in thread "main" java.lang.RuntimeException: **Cannot use index file with textual SAM file**
    at htsjdk.samtools.SamReaderFactory$SamReaderFactoryImpl.open(SamReaderFactory.java:317)
    at htsjdk.samtools.SamReaderFactory$SamReaderFactoryImpl.open(SamReaderFactory.java:145)
    at picard.sam.MergeSamFiles.doWork(MergeSamFiles.java:142)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:208)
    at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:95)

What does it mean and how do I fix it ?

Thanks

picard mergesam
• 6,370 views • link
written by AB AB

Figured it out. I only had a partial bam file. Once i downloaded it again, MergeSam ran perfectly. Easy solution to a very silly problem.

Thanks for your time Pierre.

• 0 views • link
written by AB AB

if it answered your question, click on the green mark near my answer to close the post.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

1 answer

can your please run the following command:

file file1.bam file2.bam

it should return:

gzip compressed data, extra field

if not (eg: ASCII text) your files are SAM but not BAM.

if not, I suggest your update your version of picard because there is no such method SamReaderFactoryImpl.openat line 317. https://github.com/samtools/htsjdk/blob/master/src/main/java/htsjdk/samtools/SamReaderFactory.java#L317

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thank you. I tried the command and I got this

file1.bam: gzip compressed data, extra field

file2.bam: data

Does this mean my file2 is a SAM file ?

• 0 views • link
written by AB AB

mean my file2 is a SAM file

try head file2.bamor hexdump -c file2.bam | head to see if you can guess the format...

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this
  • sorting BAM file
    written by Mohammad Amin • Mohammad Amin
    I tried sorting the BAM file i created using gatk-package-4.2.5.0-local.jar SortSam the code is : java -jar gatk-package-4.2.5.0-local.jar SortSam -I {path.bam} -O {path.sorted.bam} -SO coordinate …
  • Problem in converting SAM to BAM file by DROP SEQ
    written by VANSHIKA • VANSHIKA
    Exception in thread "main" htsjdk.samtools.SAMFormatException: Error parsing SAM header. @RG line missing SM tag. Line: @RG ID:staralign1; ; Line number 10 at htsjdk.samtools.SAMTextHeaderCodec.reportErrorParsingLine(SAMTextHeaderCodec.java:258) at htsjdk.samtools.SAMTextHeaderCodec.access$200(SAMTextHeaderCodec.java:46) …
  • Add SM group STAR
    written by Lucy • Lucy
    Hi, I have samples sequenced over three lanes that I am mapping together (paired end reads, three forward and three reverse) using STAR as follows: …
  • Errors with Picard CollectRnaSeqMetrics
    written by bertb • bertb
    Hello, I'm trying to run a QC on my bam files for an RNAseq, but when I run the command: find *U*.bam -exec echo java …
  • BWA mem giving malformed SAM headers.
    written by Ismaelrymy • Ismaelrymy
    I'm running BWA to align my reads to a big reference. I started using samtools view to obtain directly the BAM file, but when I …
  • Problems with picard markduplicates
    written by Cecelia • Cecelia
    Hi, I was running picard markduplicates with a few bam files. java -jar /sw/bioinfo/picard/2.20.4/rackham/picard.jar MarkDuplicates INPUT=sorted.bam OUTPUT=md.bam METRICS_FILE=duplicate.txt READ_NAME_REGEX=null REMOVE_DUPLICATES=true CREATE_INDEX=true And I got no …
  • Dictionary size problem in Picard CalculateHsMetrics
    written by nicola.licheri • nicola.licheri
    Hi everyone, I'm following a Roche pipeline for the analysis of targeted bisulfite sequencing data and I'm in trouble with Picard CalculateHsMetrics. I create the …
  • Picard CollectHsMetrics malfunctioning
    written by erarroji • erarroji
    I tried to run CollectHsMetrics, this is my command: java -jar /picard.jar CollectHsMetrics \ I=190.sortedDeDup.bam \ O=190_hs_metrics.txt \ R=ucsc.hg19.fasta \ BAIT_INTERVALS=AgilentSSV6_bait_list.interval_list \ TARGET_INTERVALS=V6SureSelect/AgilentSSV6_targets_list.interval_list it works …
  • SortVCF in Picard
    written by Maheswari P • Maheswari P
    i am working on SORTVcf picard tool. **java -jar picard.jar SortVcf I=filename.vcf o=output.vcf** SEQUENCE_DICTIONARY=reference.dict I am getting error as Exception in thread "main" java.lang.IllegalStateException: Key …
  • error when running MarkDuplicates
    written by kk.mahsa kk.mahsa
    when i run Picard to mark duplicates i get this error: To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp Exception in thread "main" htsjdk.samtools.SAMException: /tmp/CSPI.61127736197136018.tmp/7237.tmpnot found at htsjdk.samtools.util.FileAppendStreamLRUCache$Functor.makeValue(FileAppendStreamLRUCache.java:63) …