Biostar Test Site

This is site is used for testing only. Visit: https://www.biostars.org to ask a question.

How to do damage calibration after vg map
0
0
Entering edit mode
3 months ago

Hi all,

I am currently using vgtools on ancient DNA. My objective is to map reads to the variation graph and call variants using vg. In the past, I used bwa to do the alignments, and remove duplicates (Markduplicates) and calibrate PMD (MapDamage) on the BAM file. However, in the vg workflow, I put in two paired FASTQ file to do the mapping, and then I got a GAM file. I haven't figure out a way to put GAM into further modification.

Any suggestions on dealing with duplicates and PMD using vg?

Many thanks!

Monica

vgtool vg PMD aDNA • 103 views
ADD COMMENT
0
Entering edit mode

It seems like you want to go from gam to bam using something like this: vg surject -x x.xg -b aln.gam > aln.bam

ADD REPLY
0
Entering edit mode

Yes exactly, thank you!

ADD REPLY

Login before adding your answer.

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