call variants from cactus output
0
1
Entering edit mode
6 months ago
chelseaguoan ▴ 10

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 • 278 views
ADD COMMENT

Login before adding your answer.

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