search

smehringer / SViper

Swipe your Structural Variants called on long (ONT/PacBio) reads with short exact (Illumina) reads.
BSD 3-Clause "New" or "Revised" License
32 stars 8 forks source link

Use SViper without short reads #3

Open tungufoss opened 6 years ago

tungufoss commented 6 years ago

I want to be able to run sviper without giving it a short read bam file, to get the consensus :)

smehringer commented 6 years ago

Hi :) So you want to give it a long read BAM file and an VCF file and have sviper output a multi FASTA file containing all the consensus sequences for the appropriate positions?

tungufoss commented 6 years ago

Yes, exactly. I'm working with an tandem repeat insertion that we are worried would get polished away if we use the short reads as well.

smehringer commented 6 years ago

@tungufoss Alright :)

Last question: Do you want a flanking reference region around the consensus (analogous to the polished consensus I use)? The consensus sequence will be around 800 bp long, if you plan to remap the consensus to the reference you might want the flanking region because it will create an artificial long read (18.000bp) and will ensure correct remapping to the corresponding reference position.

hannespetur commented 6 years ago

Why don't we have both? It could be an option in the project's graceful argument parser.

smehringer commented 6 years ago

Please checkout the current develop branch to test the new feature and tell me if it works.

Happy Holidays :)