This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Qscript Example For Pipelining Gatk With Queue

I am totally a newbie to Scala/Qscripts. I want to use Queue to automate the multiple steps of running GATK (from Bam to VCF files). Could someone share their QScript examples for achieving this? I want to learn from actual examples. I've found a few examples from https://github.com/broadgsa/gatk/tree/master/public/scala/qscript/org/broadinstitute/sting/queue/qscripts/examples, but I don't think it's really helpful for building my scripts.

gatk
• 4,800 views • link
updated by Johan Johan • written by newDNASeqer newDNASeqer

1 answer

I've been working some with QScripts, in our pipeline project called Piper. The code is available here on github: https://github.com/johandahlberg/piper

If you want to start off with some simple example I'd recommend looking at MergeBamsBySample. However, there are some more advanced examples in there as well, including running the GATK best practice data processing pipeline and doing variant calling.

• 554 views • link
written by Johan Johan

Log in to answer this question.

More posts like this
  • Are the reads in the GTEx BAM/FASTQ files trimmed?
    written by Andrew • Andrew
    newbie I am want to do some RNA seq differential expression experiments using our lab-specific reference index. I need FASTQ files. I know the raw …
  • GATK for gVCF compilation issue
    written by win win
    I want to generate a gVCF file using Google tools and in order to use that a modified version of GATK has to be used. …
  • Search A Dna Sequence Motif In Bam File?
    written by newDNASeqer newDNASeqer
    <p>I need to find a short DNA sequence motif from a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file, because the wet biologist found a potential mutation from Sanger sequencing. …
  • Read File Format .Bam
    written by hudarembang • hudarembang
    <p>My research is related to the manufacture of snp database. constraints that occurs is how to read the file. <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a>. in this case, I …
  • Detecting Structural Variant From Chip Seq Data
    written by sara • sara
    <p>Hi,</p> <p>I am newbie to NGS data analysis. </p> <p>I used <a href='http://bio-bwa.sourceforge.net/'>BWA</a> for aligning to reference genome and then samtool. Now am having <a …
  • Get Random Alignments From Sam File
    written by komal.rathi komal.rathi
    <p>Hi all,</p> <p>I have a .sam file and I want to extract random alignments from it. </p> <p>I googled a bit and there are quite …
  • Paired-End Reads Alignment For Variant Calling ?
    written by newDNASeqer newDNASeqer
    <p>I'm trying to do variant calling (SNPs, Indels) from exome-sequencing data, and the sequencing was done with paired end reads. I would like to use …
  • How To Interprete/Display .Vcf Data File?
    written by newDNASeqer newDNASeqer
    <p>I am learning to do Variant Calling by using pipelines such as <a href='http://bio-bwa.sourceforge.net/'>BWA</a> + PiCard + GATK or <a href='http://bio-bwa.sourceforge.net/'>BWA</a> + <a href='http://samtools.sourceforge.net/'>Samtools</a>/Varscan. Now, …
  • What Are The (Non-Standard) Bd And Bi Tags In Samformat?
    written by Johan Johan
    <p>I'm running the <a href='https://github.com/broadgsa/gatk/blob/master/public/scala/qscript/org/broadinstitute/sting/queue/qscripts/DataProcessingPipeline.scala'>dataprocessing pipeline qscript</a> provided by the GATK people to do dataprocessing according to their <a href='http://gatkforums.broadinstitute.org/discussion/1186/best-practice-variant-detection-with-the-gatk-v4-for-release-2-0'>best practices</a>. Looking at the files …
  • Softwares For Calculation Of Genome-Wide Alignments Coverage
    written by Bioscientist Bioscientist
    <p>Any recommendation for packages to calculate alignment coverage from <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a>/sam files? 1. I need whole genome-wide calculation instead of certain genomic regions. 2. I …