SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
too small output .sai file created by bwa aln? ..leads to bwa sampe hanging? Stina Bioinformatics 12 12-05-2012 08:32 AM
Problems with BWA XA output on kbhit General 1 02-29-2012 02:16 PM
the output of bwa milkingroom Bioinformatics 1 09-21-2011 02:48 PM
Question on BWA output rboettcher Bioinformatics 3 04-18-2011 04:21 AM
bwa output scami Bioinformatics 1 12-22-2010 03:10 AM

Reply
 
Thread Tools
Old 04-30-2013, 01:48 PM   #1
DrWorm
Member
 
Location: St. Louis, MO

Join Date: Apr 2013
Posts: 17
Default BWA output

Hi all,

I'm trying to screen paired Illumina RNAseq reads against a database of putative contaminants prior to assembly. BWA (version 7.0.3a, in case it matters) ran beautifully. Now I'm trying to extract the names of all the reads that mapped to something in my contaminant database. Basically, if a read OR its mate matched a putative contaminant, I want to toss it out. A colleague told me that the output is easy to parse, given that all the unmapped reads will have a "*" in the third column of the output line. I wrote up a quick perl script to parse the output based on the "*".

Later, I learned that its possible to convert BWA's sam output to bam (via samtools) and then count the mapped reads with the following command:
samtools view -c -F4 <bam file>
(I found that here: http://seqanswers.com/forums/showthread.php?t=16743)

The issue is that I get two different answers, and I can't figure out why. What's wrong with my approach?
DrWorm is offline   Reply With Quote
Old 04-30-2013, 01:51 PM   #2
DrWorm
Member
 
Location: St. Louis, MO

Join Date: Apr 2013
Posts: 17
Default

I re-read this and wonder if the description of my perl script is unclear. I basically counted anything that did NOT have an "*" in the 3rd column. I just didn't want anyone to think I was getting inconsistent results because I had this backwards.

Thanks in advance for your feedback!
DrWorm 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 04:23 AM.


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