possible to convert hg19 coordinated position to hg38 coordinated position??
1
0
Entering edit mode
3.7 years ago
hs960201 ▴ 10

hi i'm just beginned graduate school student i have bunch of panel of normal files that made by lab predecessor(senior) who already graduated and those are hg19 coordinated. the looks like this :

#CHROM   POS   REF   ALT   DPsum;DPn;VARsum;VARn;VAFpct
1       10003   A    AC       4166;36;6;6;0.14
......

but nowadays quite many tools support hg38 and tools like Mutect2 they highly recommend using hg38 reference so i wanna change my hg19 coordinated PON file to hg38 coordinated PON file if it's possible but i can't find appropriate way i don't have any BAM or fastq file that used for PON but i really need it for VCF post filtering thanks for in advance for any reply

genome snp gene sequencing • 1.5k views
ADD COMMENT
2
Entering edit mode
3.7 years ago
Sam ★ 4.8k

What you are looking for is LiftOver

ADD COMMENT
0
Entering edit mode

Thanks for reply. i tried it yesterday with Crossmap.py but there were some unmaped region about ~2% of total VCF called region how can i handle this unmapped region? could be manually liftover it but it's gonna be real labor-intensive work. is there any solution for it?

thanks you for reply again.

ADD REPLY
0
Entering edit mode

I've never used Crossmap and have alway used LiftOver. However, it is normal to have regions that you cannot lift. You can try running LiftOver and see if you have similar level of unmapped region. Otherwise, I am not aware of any other method that can perform a 100% liftover.

ADD REPLY
0
Entering edit mode

hi sam thank you for rapid reply OK i'm gonna try LiftOver and nevertheless if they can not map that unmapped in Crossover.py results it's worth to try manually liftover for some variants

thanks have a good day

ADD REPLY

Login before adding your answer.

Traffic: 1736 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