Question regarding Liftover error fixing (GRCh38 conversion to GRCh37)

50 views
Skip to first unread message

Olivia Cao

unread,
Jun 22, 2022, 12:06:14 PM6/22/22
to gen...@soe.ucsc.edu
Hello,

Hope you are having a lovely day, I encountered several error messages when attempting to transfer a set of SNPs from GRCh38 to GRCh37. Here is a screenshot of the reported error, there are mainly three types of errors reported:
1. split in new;
2. deleted in new;
3. partially deleted in new.
 
I found a post with a solution to a similar problem (https://groups.google.com/a/soe.ucsc.edu/g/genome/c/egfMAIRfQBQ). However, our team is trying to keep the integrity of original data and convert the assembly without splitting one SNP into multiple new ones. If you would be so kind to help, the parameters are attached in a screenshot and reported error in the text file.

Hope you have a wonderful rest of today.

Best wishes,
Olivia
Screenshot 2022-06-22 at 12.37.20.png
error_694.rtf

Jairo Navarro Gonzalez

unread,
Jun 23, 2022, 7:44:32 PM6/23/22
to Olivia Cao, UCSC Genome Browser Discussion List

Hello,

Thank you for using the UCSC Genome Browser and sending your inquiry.

We do not recommend using the LiftOver tool to lift the SNPs from hg38 to hg19. Lifting a single base between assemblies is not always a trivial task as it may not get a high enough score to be considered a successful conversion.

If you have the rsIDs for all of your SNPs, we recommend pulling the hg19 coordinates from the dbSNP 153 track by searching for the rsID on hg19 or using the Table Browser. If you are unfamiliar with using the Table Browser to retrieve SNPs from the dbSNP 153 track, you can get more information from the following help page:

https://genome.ucsc.edu/goldenPath/help/hgTablesHelp.html#BatchQuery

If you wish to use the LiftOver tool, then you will have to add some bases upstream and downstream of the SNP location. However, if you choose this approach, we recommend using the LiftOver & ReMap track to look at the chains in those regions to see where they map between assemblies. You could also try to BLAT the surrounding sequence to see if there are multiple matches in hg19.

I hope this is helpful. If you have any further questions, please reply to gen...@soe.ucsc.edu.
All messages sent to that address are archived on a publicly accessible Google Groups forum.
If your question includes sensitive data, you may send it instead to genom...@soe.ucsc.edu.

Jairo Navarro

UCSC Genome Browser

Want to share the Browser with colleagues?

Host a workshop: https://bit.ly/ucscTraining


--

---
You received this message because you are subscribed to the Google Groups "UCSC Genome Browser Public Support" group.
To unsubscribe from this group and stop receiving emails from it, send an email to genome+un...@soe.ucsc.edu.
To view this discussion on the web visit https://groups.google.com/a/soe.ucsc.edu/d/msgid/genome/CAATQm6q4UCXaObVf3ZH0kFb8WkFQm83e3i04kFF6L3EfOsEzKQ%40mail.gmail.com.
Reply all
Reply to author
Forward
0 new messages