SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why so many unaligned reads inspite of "similarity" of any 2 human genomes? nadir Bioinformatics 0 04-06-2011 03:15 AM
Solexa mapping - extraction of "not in Ref" reads Oxygen81 Bioinformatics 6 11-24-2010 08:18 AM
Peak finding with MACS: Questions on # of reads and "strand bias" jjw14 Bioinformatics 0 07-29-2010 06:17 AM
Duplicate reads ("same start" reads) in 454 FLX/Titanium shotgun runs [c]oma 454 Pyrosequencing 20 08-28-2009 06:12 AM
SEQanswers second "publication": "How to map billions of short reads onto genomes" ECO Literature Watch 0 06-29-2009 11:49 PM

Reply
 
Thread Tools
Old 03-06-2013, 11:48 AM   #1
picabo
Junior Member
 
Location: San Diego, CA

Join Date: Jun 2009
Posts: 5
Default Bowtie: unaligned reads when mapping "max" reads

Hello.

I am mapping 2x100 Illumina reads against hg19 using Bowtie using a tiered approach. I want separate output files for reads that map to 1 location and reads that map to 2 locations.

In Tier 1, I set m=1 and use --max to save the reads that mapped to more than 1 location. (eg: m 1 --max max.fastq)

In Tier 2, I set m=2 and map the reads from the max.fastq file. Usually, the number of unaligned reads in this step is 0 (as would be expected). Sometimes however there are unaligned reads. Does anyone know why that is?

Thanks.
picabo 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 03:14 PM.


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