Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA Soft Clipping Bio.X2Y Bioinformatics 11 03-09-2015 10:08 AM
No soft-clipping in BWA 0.7.10 anymore? svos Bioinformatics 12 02-10-2015 09:10 AM
High Percentage of Soft Clipping, BWA-MEM, MiSeq logicthief Bioinformatics 9 09-11-2014 06:19 PM
Mapping quality and clipping in BWA WhatsOEver Bioinformatics 2 08-11-2014 11:57 AM
Soft-clipping in Stampy msalm Bioinformatics 0 01-23-2014 04:41 AM

Thread Tools
Old 01-12-2018, 07:43 AM   #1
Junior Member
Location: Brazil

Join Date: Jun 2017
Posts: 5
Question BWA -q soft clipping: Which alignment quality threshold to choose?

I wish to align to human genome samples sequenced in a Illumina HiSeq 4000. Most reads (single end) are between 20-40bp long, clean from adaptors. The base calls have 38 quality on average.

But some reads have some overrepresented repetitive kmers (like AAAAAA, CCCCCC, AGCGGGG) at the end, mostly after the 35th base, detected by FastQC. I assume these are from sequencing errors at the end of the reaction.

I'm new to sequencing and wish to align with BWA on Galaxy main. The program offers a convenient option:

"-q quality threshold for read trimming down to 35bp"

Which alignment quality threshold should I chose to soft clip the unwanted kmers at the end of some reads? Any other hint about soft clipping and other program parameters I could fine tune for this alignment? I'm currently using the defaults for these.

Last edited by Eurioste; 01-12-2018 at 07:56 AM.
Eurioste is offline   Reply With Quote

bwa, quality read alignment, repeated sequences, soft clipping

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 10:03 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO