timothyjamesbecker / SVE

Structural Variation Engine
GNU General Public License v3.0
7 stars 1 forks source link

Targets option for step 2 in advanced mode #10

Closed kyleLesack closed 4 years ago

kyleLesack commented 5 years ago

Hello,

I would like to use SVE/FusorSV for calling structural variants in some C. elegans data that I'm working with. Could you please clarify how to use to targets option or point me towards the correct documentation?

Also, regarding the FusorSV manual mentioned in the SVE documention (see below), are you just referring to the FusorSV repo's documentation (https://github.com/timothyjamesbecker/FusorSV)? I wasn't able to find the targets option instructions there.

-t or -targets option is used for target assembly for in silico SV callng validation or SV calling breakpoint refinement. This step is recommened after all calling is completed and the step (3) FusorSV data fusion and arbitration step has produce SV calls for each sample which will be passed to this -t argument. See the FusorSV manual for more information.

Thanks in advance!

timothyjamesbecker commented 4 years ago

I will have to revisit the documentation for the target option. The targets option is when you wish to use SVE and FusorSV with targeted assembly via TIGRASV. This will provide more accurate breakpoints and high confidence in those calls that TIGRASV can assemble.

kyleLesack commented 4 years ago

Thanks. I did manage to figure it out by looking through your source code. It worked after I pointed it to the Breakdancer vcf and calls files.

timothyjamesbecker commented 4 years ago

And of course, you can use first pass FusorSV calls as well.

On Tue, Feb 25, 2020 at 2:48 PM kyleLesack notifications@github.com wrote:

Thanks. I did manage to figure it out by looking through your source code. It worked after I pointed it to the Breakdancer vcf and calls files.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/timothyjamesbecker/SVE/issues/10?email_source=notifications&email_token=ABALAGOUHSEV47CZ6HSJSA3REVYZDA5CNFSM4JGBWHH2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEM5H25Q#issuecomment-591035766, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABALAGLQC4XCEQOEFIISU7LREVYZDANCNFSM4JGBWHHQ .

kyleLesack commented 4 years ago

Thanks for the tip. I was passing both a .calls and .vcf file to Tigra. Will it work with just a .vcf file (I'm not seeing a .calls file in the FusorSV output)?

timothyjamesbecker commented 4 years ago

It will work with FusorSV vcf calls the BD.call files and the BD vcf calls that the SVE converts from .call format. The FusorSV calls have less false positives so the runtime is greatly reduced compared to BD in human data. FOr humans this target assembly step can take up a day to run on a human for FusorSV vcf and up to a week on BD calls.

On Tue, Feb 25, 2020 at 3:39 PM kyleLesack notifications@github.com wrote:

Thanks for the tip. I was passing both a .calls and .vcf file to Tigra. Will it work with just a .vcf file (I'm not seeing a .calls file in the FusorSV output)?

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/timothyjamesbecker/SVE/issues/10?email_source=notifications&email_token=ABALAGIQOBVBDPQ5C7GIZS3REV6YVA5CNFSM4JGBWHH2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEM5NETQ#issuecomment-591057486, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABALAGKUYNAXK52MVYVR65DREV6YVANCNFSM4JGBWHHQ .