Question: How do I process multiple pileup files at once with Varscan to get the vcf files?
0
gravatar for ppurkayastha2010
7 months ago by
ppurkayastha201030 wrote:

I have a list of 100 SRA IDs for which I need to get the VCF files using VARSCAN.

VARSCAN takes an input in pileup format.

I could get the pileup files for all 100 SRAs using NCBI SRA Tools. But not able to give the List of pileup files to Varscan, in order to generate VCF files.

Is there a way I could generate VCF files from all 100 pileup files at one go or a way to generate VCF files directly from NCBI SRA ?

any advice and suggestions will be appreciated..

sra snp varscan pileup • 329 views
ADD COMMENTlink modified 7 months ago by Jennifer Hillman Jackson25k • written 7 months ago by ppurkayastha201030
0
gravatar for Jennifer Hillman Jackson
7 months ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

To make calls for multiple samples, Varscan expects merged pileup output from Samtools Mpileup with read groups assigned. This is different from what is extracted from NCBI (individual samples). Instead, download the fastq data from NCBI, map with BWA, run Mpileup, then Varscan. Dataset collections can be used to process the data in batch mode.

Galaxy example tutorials that cover how to use collections and perform variant analysis: https://galaxyproject.org/learn/

Most Varscan command line functions/parameters are incorporated into the Galaxy wrapper. The same is true for the upstream tools also wrapped for Galaxy. This means that if a particular use-case is not covered by a tutorial, the original tool's manuals can be useful references. Varscan user manual: http://varscan.sourceforge.net/using-varscan.html#pileup-howto

Thanks! Jen, Galaxy team

ADD COMMENTlink written 7 months ago by Jennifer Hillman Jackson25k
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 16.09
Traffic: 136 users visited in the last hour