SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SOLiD paired end data and BWA - cleaning/mapping Kath Bioinformatics 4 09-20-2013 04:24 AM
Error found while mapping Illumina data using BWA nirav99 Bioinformatics 15 01-14-2013 06:30 PM
BWA align SOLiD PE data gives poor mapping & 0.5% properly paired alig Bioinformatics 3 07-08-2011 08:44 AM
BWA behaviour with Mate Pair data + Multi read mapping apratap Bioinformatics 5 06-23-2011 04:38 PM
BWA mapping fastq files with Illumina quality maricu Bioinformatics 3 11-19-2010 11:18 AM

Reply
 
Thread Tools
Old 01-02-2014, 07:33 AM   #1
Genomics101
Member
 
Location: Maryland, USA

Join Date: May 2012
Posts: 60
Default Mapping multiple FASTQ data types using BWA?

Greetings!

I plan to do mapping using multiple types of Illumina FASTQ files for the same sample.

I have one set of FASTQs made from a library with 400bp inserts, 100bp reads (trimmed for quality), and processed with the 1.5 pipeline; the second set is 600bp inserts, 250bp reads (trimmed for quality), and processed with the 1.9 pipeline. Because of the trimming process, I have paired-end and single-end FASTQs for each sample as well.

As a result I have 8 FASTQs for mapping, e.g. -

Insert size 400:
Sample1_Illumina_1.5_100bp_R1.fastq.pe
Sample1_Illumina_1.5_100bp_R2.fastq.pe
Sample1_Illumina_1.5_100bp_R1.fastq.se
Sample1_Illumina_1.5_100bp_R2.fastq.se

Insert size 600:
Sample1_Illumina_1.9_250bp_R1.fastq.pe
Sample1_Illumina_1.9_250bp_R2.fastq.pe
Sample1_Illumina_1.9_250bp_R1.fastq.se
Sample1_Illumina_1.9_250bp_R2.fastq.se

I have only been mapping with the 100bp/1.5 so far, using BWA, but reading it looks like I should use BWA-MEM if I am using the longer reads. Is it possible to map all of these together?

Thanks!
Genomics101 is offline   Reply With Quote
Old 01-02-2014, 11:54 AM   #2
adamdeluca
Member
 
Location: Iowa City, IA

Join Date: Jul 2010
Posts: 95
Default

I would map separately for each insert size, then merge the resulting bams for further analysis.
adamdeluca is offline   Reply With Quote
Old 01-02-2014, 11:55 AM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

No, don't map them together. Also, bwa mem doesn't support illumina 1.5.
dpryan is offline   Reply With Quote
Reply

Tags
bwa, bwa-mem, fastq, paired end mapping, single end 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:43 AM.


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