Question: Using hg38 assembly with hg19 exome capture data
gravatar for jdm204
15 days ago by
jdm2040 wrote:

Hello all,

I have some exome sequencing data which used an exome capture kit ( 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,


ADD COMMENTlink written 15 days ago by jdm2040

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 15 days ago by finswimmer10k
Please log in to add an answer.


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