soft-clipping during alignment and variant calling query · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

soft-clipping during alignment and variant calling query

It is known that BWA-MEM soft-clips bases to achieve a higher alignment rate. If a variant caller is used after that, say, GATK Haplotyecaller, is it necessary to keep the parameter --dont-use-soft-clipped-bases to TRUE? Or should we keep it in default (false) as it may affect the realignment step during Haplotycaller run?

bwa-mem gatk variant calling alignment

• 2,390 views • link

updated 6.7 years ago by Biostar

• written 6.8 years ago by prasundutta87

Please stick to the GATK best practices. If they don't recommend to use that parameter there is no need to overthink.

• 0 views • link

written 6.8 years ago by WouterDeCoster

0 answers

No answers yet.

Log in to answer this question.

More posts like this

BWA MEM paramenter
written 6.5 years ago by deniselavezzari •

Hi everyone! I have a question about BWA-mem parameters. I know that is usually better to use default ones, but I want to know if …
Soft-clipping of reads in Amplicon-sequenced data
written 9.0 years ago by jsneaththompson

I have a variant calling pipeline which I use to process amplicon-sequenced fastq files; it uses cutadapt to remove the adapter sequences on either the …
Tool: GRIDSS: the Genomic Rearrangement IDentification Software Suite
written 9.4 years ago by d-cameron

GRIDSS is typically used for detecting structural variation breakpoints from short read sequencing data but is a modular software suite containing a number of tools …
BWA and Soft clipping
written 9.8 years ago by christopher.patterson3 •

Hello, I'm new to the forum here. I am working with a set of Amplicon sequences produced from Illumina (paired-end, read length 251 bp) and …
What are soft clipped reads??
written 10.3 years ago by le.marcorin •

Hi! I'm working with Illumina Miseq reads and i'm having some trouble with variant calling. I used cutadapt for trimming adapters, bwa for alignment and …
No soft-clipping in BWA 0.7.10 anymore?
written 11.5 years ago by Sebastian •

Hi all! After moving from bwa version 0.5.8c to 0.7.10 I discovered differences that look like errors in the newer version: it seems as BWA …
Soft-Clipping And Variant Calling
written 12.4 years ago by DoubleDecker

<p>Having aligned my reads with <a href='http://bio-bwa.sourceforge.net/'>BWA</a>-mem, I noticed that a lot of alignments are dodgy, to put it mildly, e.g. with both ends massively …
Why Is Bwa Soft-Clipping Exactly Matching Reads?
written 12.7 years ago by ktpetrie •

I am using bwasw to align reads from a MiSeq run to a very short reference sequence. For some reason, [bwa](http://bio-bwa.sourceforge.net/) is soft-clipping (assigning S …
Disable Soft-Clipping Bwa Aln
written 13.5 years ago by Irsan

<p>Is there a way to tell <a href='http://bio-bwa.sourceforge.net/'>bwa</a> aln to do not soft-clip reads and try to align them? </p> <p>We see that many soft-clipped …
Bwa Soft-Clipping Pattern In Read Length Distribution Using "-Q 20" Parameter
written 13.8 years ago by toni

<p>Hi all,</p> <p>I just had to analyze a really bad lane. It is an old Illumina GAIIx lane with reads of 152 cycles. The last …