Advantages of split by chromosome?
0
0
Entering edit mode
6 weeks ago
kwanghoon ▴ 20

Hi, I'm a baby at this field, and trying to do some tutorials like GATK best practice, CNV analysis, HLA-typing etc.

I have a question about files that are split by chromosome.

What advantages do they have?

And If I use them, is my command correct like below?

ex) bwa mem REF_Chr1.fa A_1.fastq A_2.fastq > A.sam

After many process, do I have to merge them at the end?

Thank you!

cnv split gatk • 269 views
ADD COMMENT
1
Entering edit mode

What advantages do they have?

parallelisation.

ADD REPLY
0
Entering edit mode

is my command correct like below?

ex) bwa mem REF_Chr1.fa A_1.fastq A_2.fastq > A.sam

No never map a bam on a contig only. you'll get false positives. You must use the whole reference genome.

ADD REPLY
0
Entering edit mode

Then do I need to split after bam processing using whole reference?

Sorry I have no knowledge.

ADD REPLY
1
Entering edit mode

you don't have to _split_ the bam. For example, for calling variants, most tools are able to call a specific region/chromosome if the bam is indexed (=bai file)

ADD REPLY

Login before adding your answer.

Traffic: 850 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6