how to compare mapping of WES samples to human pangenome?
0
1
Entering edit mode
16 months ago
bitpir ▴ 250

Hi,

I'm still trying to wrap my head around the new human pangenome reference and would like some advice on how to go about analyzing some of the WES (hg38/hg19 baits) that I currently have.

  1. How should I calculate the coverage of a specific gene, say, BRAC1, based on the human pangenome annotation? Would it make sense to compare the coverage at all or should I be focusing on the coverage of the best ancestry-matched pangenome reference?

  2. Is there a visualization tool to show how reads are mapped to the pangenome at certain gene? I am envisioning the 47 pangenomes being displayed in each row and reads are piled up under each reference.

Thanks so much for the guidance and advice!

WES pangenome vg • 1.1k views
ADD COMMENT
1
Entering edit mode

I don't know of any tool to visualize all of the haplotypes simultaneously like you are describing. The most capable visualization tool for mapped reads is probably sequence tube map, but it's a graph-based layout. You can use it query reference regions if you want to look at a particular region. If you want to calculate coverage, there are tools to do that in vg pack, but first you'll need to convert the reference region into node IDs.

ADD REPLY
0
Entering edit mode

great response - in your first sentence you mention not knowing of such a tool - in what ways does the vg sequence tube map you link to not do this?

ty

ADD REPLY
0
Entering edit mode

You will not see all of the haplotypes separately as independent, linear tracks, each with its own read alignments. Instead, the (distinct) haplotypes are shown as walks through a graph. Also, the alignments are shown one time (rather than once for each haplotype), also threaded through the graph.

ADD REPLY
0
Entering edit mode

Thank you for the response! What is the best way to convert the reference into node ID's? Is there a liftover tool to convert hg38's coordinates to the new reference?

ADD REPLY
0
Entering edit mode

I also aligned my reads to the human pangenome reference, and I am interested at the HLA area and will also be happy to hear how to evaluate the coverage, visualize the mapping and create a vcf from my alignments

ADD REPLY

Login before adding your answer.

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