Question: VirusClip script Virus_clip.sh
0
gravatar for KVC_bioinfo
18 months ago by
KVC_bioinfo380
Boston
KVC_bioinfo380 wrote:

Hello,

I am trying to use VirusClip for one of my projects. However, I am not sure what exactly we are supposed to provide in the following part of the virus_clip.sh script:

"resources"

refFasta=/path/to/hbv_ref/hbv.fa

blastdb=/pathool01/disk1/software/ref/hg19

annovardb=/pathool01/disk1/data/danielho/ireneng/ANNOVAR/annovar/humandb

If anyone has used virusclip before:

  1. can someone help me with this issue.
  2. also, how long did it take if you have used Virusclip before.

    any help will be greatly appreciated. Thank you!

#rna-seq #virusclip • 459 views
ADD COMMENTlink modified 18 months ago • written 18 months ago by KVC_bioinfo380
0
gravatar for genomax
18 months ago by
genomax65k
United States
genomax65k wrote:

What part of that do you find confusing?

refFasta = Path (use real paths from your environment) to the fasta file containing sequence of virus of interest (you should have this)
blastdb = Path to the blast database for human genome (in this case hg19, but you can use any other) (you will need to make the index)
annovardb = Path to the annovar human database that you can download following the instructions on this page.

BTW: Don't use # signs in your post tags.

ADD COMMENTlink modified 18 months ago • written 18 months ago by genomax65k

Thank you. I did not understand annovardb. what exactly that database is?

ADD REPLYlink written 18 months ago by KVC_bioinfo380
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: 1133 users visited in the last hour