Question: Using hg38 assembly with hg19 exome capture data
1
gravatar for jdm204
10 months ago by
jdm20410
jdm20410 wrote:

Hello all,

I have some exome sequencing data which used an exome capture kit (http://emea.support.illumina.com/downloads/nextera-rapid-capture-exome-v1-2-product-files.html) from Illumina which is based on hg19 coordinates.

I aligned this data to hg38 and called variants etc which worked out fine. However, I now want to look at copy number, which requires explicitly providing the exome regions that were targeted. Since this information is based on hg19 coordinates and my alignments (.bam) are based on hg38, I'm assuming I can't simply ignore the difference and hope to get meaningful results.

My question is, is it perfectly OK to simply lift over the targeted regions .BED file to hg38 (though the liftover tool reports that some regions are lost) or, is it a general rule that if my exome capture is based on hg19 coordinates I should be sticking to the hg19 assembly for analysis?

Many thanks,

Jamie

ADD COMMENTlink written 10 months ago by jdm20410
3

It's absolutly fine to liftover the bed file to hg38. Your DNA doesn't care about whether it is hg19 or hg38 ;)

But as you said, you have to live with the fact that some regions cannot be liftover.

ADD REPLYlink written 10 months ago by finswimmer13k
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: 873 users visited in the last hour