This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bedtools Multicov Need A Bam Index File Specification Option

bedtools version 2.16.2 multicov used to compute the multiple sample coverage given a feature file(gtf bed).

format: bedtools multicov -bams alin1.bam aln2.bam .. -bed capturRegion.bed >out.coverage

official doc has mentioned that input bam files should be sorted and indexed, but it does not mention the details. suppose the bam file name is: sample1.bam, then the index file should be named: sample1.bam.bai(not sample1.bai) ,otherwise multicov will report an error: indexes not found.

I think it would be better to add an option which will allow the user to specify the bam index files or the suffix used for these index files.

bedtools coverage bam sort index
• 7,914 views • link
updated by Alex Reynolds Alex Reynolds • written by ewre ewre

1 answer

Like many things in bioinformatics this comes from a convention.

  • This is the default name that is generated from samtools index
  • many tools follow this convention (bedtools, IGV, IGB, UCSC etc...)
  • if you follow this convention like everybody, life will be easier
  • if at all this should be a feature request in bedtools not a "Tip"
• 62 views • link
written by Ido Tamir Ido Tamir

can not agree more

• 0 views • link
written by Tky Tky

thanks for the reply and the "feature request but tip" suggestion. well, i can not agree more with that life would be easier when we follow some convention. but there are always exceptions, picard output bam index file with the default form sample1.bai. and probably this is much better when the sample file name is long(LHX?_mix?_index?.recal.reduce.xxxx). I just mention it here to avoid someone who is newer like me in sequence file handling spending time googling the "indexes file not found" error.

• 0 views • link
written by ewre ewre

Log in to answer this question.

More posts like this
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • Using Gnu Parallel For Bedtools
    written by GouthamAtla GouthamAtla
    <p>I am trying to run <code>gnu:parallel</code> on <code>bedtools multicov</code> function where the original command is </p> <pre><code>bedtools multicov -bams bam1 bam2 bam3.. -bed anon.bed &gt; …
  • Conversion From Sam To Bed With Max Information
    written by roll roll
    <p>I am trying to convert my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file into bed file in such a way that i will keep some information from the <a …
  • Samtools View Of A Sam File Entire Chromosome
    written by dfernan dfernan
    <p>Hi, </p> <p>I have (what I think) is a valid <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file.</p> <p>When I do:</p> <pre><code>samtools view -S F1i1F.C57.2.Aligned.out.sam chr19_C57 </code></pre> <p>I get the …
  • Merging Bam Files
    written by rickyflintoff • rickyflintoff
    <p>I have two <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files that have a different sample in read group (@RG) section. I would like the sample to be the same. …
  • How To Rearrange Paired End Bam File?
    written by deepthithomaskannan deepthithomaskannan
    <p>Hello all,</p> <p>I have a paired end <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file and I want to use <a href='https://code.google.com/p/bedtools/'>bedtools</a> for them. After merging, the paired end read …
  • Converting Sam Files To Bam Files - Reproduce Results Nature Paper: Transcriptome Genetics Using Se…
    written by Sander Timmer Sander Timmer
    <p>I want to reproduce the results that people achieved in the following Nature paper: Transcriptome genetics using second generation sequencing in a Caucasian population <a …
  • Getting Unmapped Reads: Comparing Fastq To Bam
    written by User 9996 User 9996
    <p>given a FASTQ file and a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file of aligned reads, is there an efficient way to get all FASTQ reads that are in …
  • Converting Bam Or Sam Files Into Bed Files
    written by Nancy • Nancy
    <p>hi there,</p> <p>how can I convert <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> or <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> files into bed files?</p> <p>thanks so much for your help!!!!</p> <p>Nancy</p>