SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Number of reads in paired end fastq files bong28 Bioinformatics 2 04-08-2018 10:53 PM
Paired end fastq files contamination HSmith General 2 12-02-2016 10:05 AM
Merging Paired-End FastQ Files jmpi Bioinformatics 3 05-22-2013 07:21 AM
splitting big paired fastq files JahnDavik Bioinformatics 2 04-14-2013 11:16 PM
Split Large FASTQ file in small FASTQ files with user defined number of reads Windows deepbiomed Bioinformatics 3 04-04-2013 07:14 AM

Reply
 
Thread Tools
Old 06-21-2019, 11:39 AM   #1
lorenzo.C
Junior Member
 
Location: berkeley

Join Date: Jan 2018
Posts: 4
Default isolate reads from a paired-ended fastq files

hi folks,
I have pair ended fastq files (let's call them 01_R1 and 01_R2). I need to make a new pair of fastq files that contains just the reads that fulfill this requirement: the reads of R2 should be within a certain range (let's say no more than 100bp far) from the respective read in R1.
any suggestion on how to approach this?
lorenzo.C is offline   Reply With Quote
Old 06-26-2019, 04:20 AM   #2
bluesmartini
Junior Member
 
Location: Montpellier

Join Date: Jun 2014
Posts: 4
Default

Hello,

I would say that most of the aligners have this option. Which aligner do you want to use?
bluesmartini is offline   Reply With Quote
Old 06-26-2019, 09:56 AM   #3
lorenzo.C
Junior Member
 
Location: berkeley

Join Date: Jan 2018
Posts: 4
Default

BOWTIE2 would be great! thanks for the help!
lorenzo.C is offline   Reply With Quote
Old 06-27-2019, 12:10 AM   #4
bluesmartini
Junior Member
 
Location: Montpellier

Join Date: Jun 2014
Posts: 4
Default

with Bowtie2, the (very extensive) documentation mentions the following options :

--al-conc <path>
Write paired-end reads that align concordantly at least once to file(s) at <path>. These reads correspond to the SAM records with the FLAGS 0x4 bit unset and either the 0x40 or 0x80 bit set (depending on whether it’s mate #1 or #2). .1 and .2 strings are added to the filename to distinguish which file contains mate #1 and mate #2


So as output you will have 2 fastq files that contain the reads that align concordantly to your options.

For the range between pairs, take a look at --maxins option.

hope it will help

aubin
bluesmartini is offline   Reply With Quote
Old 06-27-2019, 11:22 AM   #5
lorenzo.C
Junior Member
 
Location: berkeley

Join Date: Jan 2018
Posts: 4
Default

thanks for the help good sir, if I understand correctly --al-conc <path> will produce a new couple of fastq that fulfil the --maxins option. the maxins command will take just reads that are spaced to a maximum number of base pairs in between. the problem I have is that the second read should have the opposite orientation of the first read to be selected by maxins (if I understand correctly), is there a way to not have this option selected? in other words I need to isolate reads that are within a certain range of the other direction, no matter the orientation (both should be considered valid).
lorenzo.C 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 05:38 PM.


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