Biopython for NGS
Can I use Biopython to automate the NGS pipeline. If yes, which Biopython library can be used to run GATK, SAMTOOLS, bwa, PICARD. ? I hope my question is clear.
sequencing
next-gen
snp
• 2,668 views
•
link
written
by
shuksi1984 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
BWA alt-aware pipeline troubleshooting
written by Duy •Hi everyone, I have a question regarding the BWA alt-aware pipeline since I haven't seen much guideline online other than [https://gatk.broadinstitute.org/hc/en-us/articles/360037498992--How-to-Map-reads-to-a-reference-with-alternate-contigs-like-GRCH38#2][1] (though it might a …
-
Strange speed up in GATK LeftAlignIndels
written by jonas.andersson •Hi! I noticed a strange thing, I have been running a DNA-seq pipeline like this: reads -> bwa-mem2 -> picard SortSam -> picard MergeSamFiles -> …
-
docker image building error
written by leukippus0116 •Hello, I'm working on NGS data analysis and now I'm trying to build docker image of Variant Calling Pipeline using GATK4(https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/). I have encountered an …
-
Chromosome position in VCF file
written by tothepointI am working on wgs dataset and after running bwa, samtools, picard, gatk pipeline and filtering snp, indels got vcf file for further analysis. Now …
-
AddOrReplaceReadGroups (picard) vs reheader (samtools)
written by shuksi1984 •What is the difference between output generated by AddOrReplaceReadGroups (picard) and reheader (samtools)? Can AddOrReplaceReadGroups (picard) be used instead of reheader (samtools) or vice versa …
-
Unit-testing for NGS data analysis pipeline
written by NBHello, I was wondering if anyone does unit-testing for NGS pipeline that follows GATK best practice and is written in bash ? From what I …
-
Help in writing NGS pipeline (or: the road from fasta to vcf)
written by Hervé •Hi, for various reasons I decided to try to understand better the variant calling process. I am interested in human data, however I decided to …
-
Bwa Mem -M Option
written by Varun Gupta<p>Hello</p> <p>I want to align my reads to the reference genome and then use GATK for downstream analysis. I am using <a href='http://bio-bwa.sourceforge.net/'>bwa</a> mem protocol.</p> …
-
Whole Genome Sequencing Pipeline: Fixmate
written by mathiasf52 •<p>Hello, </p> <p>I have a question regarding WGS pipeline. should we use Fixmate for WGS? if yes, which software is the fastest? PICARD FixMateInfo.jar seems …
-
Bwa And Gatk Runtimes
written by aniketd86I am trying to optimize a NGS pipeline for Human Genome sequence analysis. I would like to know the estimated runtimes for bwa and gatk. …
For automated pipelines, people use what is known as snakemake
Creating workflows with snakemake and conda
Although, I frequently use shell scripting
While browsing, I also came across Ruffus. How about Ruffus? Which one would be better sankemake or Ruffus? Although, both are python library but there must be some difference
I'd suggest snakemake. As far as I know it's more widely used.
Other options you could explore are NextFlow/FlowCraft