Question: confused about how to run pipeline
0
gravatar for andrew.lamb
11 months ago by
andrew.lamb0 wrote:

I'm confused about cnvkit batch. I ahve one normal bam and one tumor bam file(Exome SEQ).

I have a bed file to use for --targets, and flat ref text to go with it. I don't have an access file, or anti-target file. Do I need to use the access command first to create an access file or will batch do that with this command:

cnvkit.py batch 
Tumor.bam \
--normal Normal.bam \
--targets targets.bed \
--fasta Homo_sapiens.GRCh38.dna_sm.primary_assembly.fa \
--output-reference my_reference.cnn \
--output-dir example/ \
--annotate fer.text

Thanks!

cnvkit • 557 views
ADD COMMENTlink modified 11 months ago by toralmanvar750 • written 11 months ago by andrew.lamb0

Hello andrew.lamb,

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.
code_formatting

Thank you!

ADD REPLYlink modified 11 months ago • written 11 months ago by h.mon24k
1
gravatar for toralmanvar
11 months ago by
toralmanvar750
toralmanvar750 wrote:

You do not require access or antitarget file, because batch command runs all necessary subcommands to generate access and antitarget file.

Here they have mentioned, which subcommands are executed while using batch command.

ADD COMMENTlink written 11 months ago by toralmanvar750
1

I think there's an error in the documentation because it shows that you have to pass an access file but first you need to run cnvkit.py access to create it. Do what Paul suggests in his answer.

ADD REPLYlink modified 4 months ago • written 4 months ago by nick90
0
gravatar for Paul
11 months ago by
Paul1.3k
European Union
Paul1.3k wrote:

Hi,

first please create your access file (need to download exclude regions file):

cnvkit.py access hg38.fa -x excludes.bed -o access-excludes.hg38.bed

Then you can use batch command:

cnvkit.py batch --target-avg-size 70 -p 12 TUMOR/*.bam --normal NORMAL/*.bam --targets capture_region.bed --fasta hg38/genome.fa --access access-excludes.hg38.bed --output-reference try_reference_bed_70.cnn --output-dir output_splitted_bed_70

You can set up target average size - default is I think 200.

Paul.

ADD COMMENTlink written 11 months ago by Paul1.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 929 users visited in the last hour