Question: call variants from cactus output
1
gravatar for chelseaguoan
29 days ago by
chelseaguoan10
chelseaguoan10 wrote:

I am new on vg tools. I have 100+ fungi assemblies and raw reads. The assemblies were aligned by cactus and the output hal file was converted to vg graph by hal2vg. I the vg github page, it seems there was 2 ways to call variants.

  1. mapping reads to the vg graph and use vg call to get variants. Do I need to map the reads of each isolates to graph individually? and get vcfs of each isolates, then merge them together?
  2. use vg deconstruct directly to get the vcf. This seems suitable for my cactus output?

I do not know if my understanding is right? What is the difference of the two methods? Especially for the first method, is it right to map raw reads to graph individually?

Thank you very much!

vg • 107 views
ADD COMMENTlink written 29 days ago by chelseaguoan10
Please log in to add an answer.

Help
Access

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