vg surject for long reads
Can vg surject be used for long reads like ONT or pacbio ?
• 1,894 views
•
link
updated
by
Ram
•
written
by
ved_vyas •
•
written
by
ved_vyas •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
minimap2 for polyA site sequenced with long read data (both cDNA ONT, PacBio)
written by shinyjj •
-
How many sequences per treatment is required for long reads in hybrid approach of microbial metagen…
written by ranatekan91 •
-
BAM creation - vg surject vs vg mpmap output
written by AshleeThomson •
-
GraphAligner/vg surject - alignment output nonsense?
written by Andrew •
-
Order of reads in Long read FASTQ file
written by shinyjj •
-
Long read RNA # of reads interpretation help
written by shinyjj •
-
vg add introduces cycles to a directed acyclic graph
written by ved_vyas •
-
Polishing PACBIO corrected assembly using ONT and ILLUMINA data (.fastq)
written by K •
-
Segmentation fault (core dumped) during vg surject
written by LilithElina •
-
Tools for metagenomics classification using long reads
written by Medhat
I can't say that we've extensively characterized the performance with long reads, but in principle yes, it should work. As of right now, VG itself doesn't have a mapping tool that works well for long reads (although we probably will in the future). However, you could perhaps use alignments from GraphAligner. One speed bump we've discovered in the past is that GraphAligner's GAF output isn't compatible with vg surject, so if you go that direction you should use GAM as the alignment format.
Thank you for the prompt reply, I did exactly what you have described, but got the following error:
warning[vg::Surjector]: Refusing to perform very large alignment against 119596 bp strand split subgraph for read Read_name ; suppressing further warnings. terminate called after throwing an instance of 'std::out_of_range' what(): basic_string::substr: __pos (which is 22) > this->size() (which is 1)
Originally I had the gfa file, which I converted to vg and xg index using "vg convert" module: In GraphAligner vg index was used for mapping and in vg surject xg index was used. Command used for gam to bam conversion: vg surject -x ${xg_file} -p chr21 -t 20 -b ${gam_file} > ${bam_file}