After I use the "-q 20" option in bwa aln I start getting errors down the road eg when i run ValidateSamFile in picard I get lots of "Mate negative strand flag does not match read negative strand flag of mate" and then "Sequence and quality are inconsistent" errors as well. I have run all the same commands without the -q option without these problems. I am very new to all this and not sure where I'm going wrong!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
This is a bug in BWA. If you use -q it puts the trimmed nucleotide read in the output SAM file but the quality sequence is not trimmed at all and thus there is a difference in the length of the read and the length of the quality string. Either use some third party tool like trimgalore to trim your fastq reads before aligning and don't use BWA -q option at all OR use some script that takes SAM file and trims off the base quality string accordingly.
PS:This post may help- http://www.biostars.org/post/show/13...or-space-data/
Comment
Latest Articles
Collapse
-
by seqadmin
The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...-
Channel: Articles
02-24-2025, 06:31 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 03-03-2025, 01:15 PM
|
0 responses
173 views
0 likes
|
Last Post
by seqadmin
03-03-2025, 01:15 PM
|
||
Started by seqadmin, 02-28-2025, 12:58 PM
|
0 responses
265 views
0 likes
|
Last Post
by seqadmin
02-28-2025, 12:58 PM
|
||
Started by seqadmin, 02-24-2025, 02:48 PM
|
0 responses
651 views
0 likes
|
Last Post
by seqadmin
02-24-2025, 02:48 PM
|
||
Started by seqadmin, 02-21-2025, 02:46 PM
|
0 responses
266 views
0 likes
|
Last Post
by seqadmin
02-21-2025, 02:46 PM
|
Comment