SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA-MEM: output mapped reads larger than input reads Kennels Bioinformatics 9 09-05-2013 06:27 PM
BWA interpretation Inma Bioinformatics 2 05-09-2013 07:34 PM
how do I output the CS tag for BWA align of SOLID reads? KevinLam Bioinformatics 16 07-23-2011 11:06 PM
Does output-file of bwa-sw contain unmapping reads? louis7781x 454 Pyrosequencing 1 04-06-2011 08:25 AM
sam output from bwa for SOLiD reads in colorspace? nisha SOLiD 19 01-07-2010 05:05 AM

Reply
 
Thread Tools
Old 10-16-2014, 09:38 AM   #1
gwilymh
Member
 
Location: Milwaukee

Join Date: Dec 2011
Posts: 72
Default Interpretation of bwa output - has it used all the reads

I am analyzing an RNASeq data set. I have screened reads in Trimmomatic and aligned them using bwa mem. The reads are mostly 100bp long, but a minority are as short as 30bp. Reads are submitted in two fastq files and are paired end. The input files each have a total of 61,383,869 reads.

I get the following message in the bwa standard output:
[M::main_mem] read 109456 sequences (10000139 bp)...
[M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (45, 19766, 45, 52)

Can anyone tell me precisely what the above means? Why did bwa read only 109,456 sequences from the millions in the fastq files?

As only 19,908 "unique pairs" in total were reported, does this mean that the rest of the read data were all PCR duplicates?
gwilymh is offline   Reply With Quote
Reply

Tags
bwa, fastq, pcr duplicates, reads

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 09:01 PM.


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