SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
cufflink doesn't like sam from bwasw lukas1848 Bioinformatics 1 01-16-2012 05:11 AM
Newbler Trim Status blindtiger454 De novo discovery 2 05-18-2011 04:46 AM
GATK - doesn't see files Kath Bioinformatics 3 12-08-2010 08:38 AM
Do I need to trim the sequences like this? days369 Bioinformatics 4 08-16-2010 08:19 PM
Soap2: the -w command doesn't work Yushan Hsiao Bioinformatics 0 08-10-2010 09:54 PM

Reply
 
Thread Tools
Old 11-27-2011, 04:52 PM   #1
CNVboy
Member
 
Location: boston

Join Date: Jun 2011
Posts: 27
Default BWA trim doesn’t work.

I tried to run BWA with trimming those low-quality bases around the end of read. And command is like:
Code:
/share/bin/bwa aln -t 7 -I -q 20 ****.txt.gz > ****.sai
Here -I is to convert Illumina format to sanger
However when I check the generated sam file, I can still see lines like:
Code:
HWI-ST150_0131:3:1:5473:1943#0  133 X   113713314   0   35M5S   =   113713314   0   TNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN    ########################################    XC:i:35
Shouldn't the crappy sequences trimmed? Or my command has some problems?thx
CNVboy is offline   Reply With Quote
Old 11-28-2011, 02:57 AM   #2
cjp
Member
 
Location: Cambridge, United Kingdom

Join Date: Jun 2011
Posts: 58
Default

In here for 133 - the samtools flags: http://picard.sourceforge.net/explain-flags.html

you get:

Summary:
read paired
read unmapped
second in pair

So the read is unmapped:

Also see this thread from today about bwa-mapped reads with the unmapped flag:

http://seqanswers.com/forums/showthread.php?t=15772

Chris

Last edited by cjp; 11-28-2011 at 02:58 AM. Reason: explain what 133 is
cjp is offline   Reply With Quote
Reply

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:54 AM.


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