Bedtools multicover how to submit a bam list
Hi,
I'm trying to use bedtools muticov to find out Copy Number Variations among my samples, using the CANOES soft. I got many bam samples and canno't figure out how to "impute" a list of bams to bedtools multicov.
Already tried the following
``
find *.bam | xargs bedtools multicov -bed probes.bed -q 20 > canoes.reads.txt
An other way around:
find *.bam > bam.list
bedtools multicov -bams bam.list -bed probes.bed -q 20 > canoes.reads.txt
CANOES tutorial suggests the following form, but I can't understand how a quote around bam list could enhance things.
bedtools multicov -bams `cat bamlist.txt` -bed probes.bed -q 20 > canoes.reads.txt
Thank you eveyone!
• 3,310 views
•
link
written
by
benoit.tessoulin •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Difference between mapped reads and counts
written by josmantorres •
-
Downstream analysis of VCF files obtained from VarScan2
written by Raheleh
-
Running Freebayes on 100+ Samples
written by tpaisie
-
Retain file names as sample names while making pileup
written by bioinfo8
-
MACS2 coverage calculations
written by ariel.balter
-
How to convert bam-readcount output to BED-VAF format
written by Haitao •
-
Extracting counts from RNA-Seq: trouble with bedtools + multicov
written by Kristin Muench
-
Identify genes in a tophat aligned bam/sam file
written by gchaves •
-
SNP zygosity from RNA-Seq data (BAM file)
written by ChIP
-
MeDip-seq data interpretation + analysis / MEDIPS, MACS2
written by florian.noack •
The ` (backtick) causes the system command to run with the output being passed to program. They are there for a reason.
Thanks I got it for the backtick, but still doesn't do the job here...
And what's wrong with just following the tutorial? :-P
hehe, it doesn't work when I follow the tutorial :-P Gonna ask to bedtools team if it's bedtools' version issue. But i'm still surprised that other way around don't make it work out.