Alignment of single/paired end reads with rsubread align
Hi biostars,
align function by Rsubread package allows to set parameters according to reads layout (single or paired), since i download fastq files from SRA, if the layout is paired, i obtain a single fastq files containing both pairs.
The question is: if i don't specify if fastq file is paired or not, align would detect "automatically" this and consequently align reads correctly? or align the reads as if it were a single layout?
any reply will be appreciated, thanks.
• 1,509 views
•
link
written
by
qstefano •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Adapting nf-core/sarek for Ion Torrent Single-End Data Analysis
written by Nour-eddine •
-
How many FASTQ files does one SRA Run correspond to?
written by hamarillo •
-
retrieving paired end sequencing data with fasts-dump
written by v.berriosfarias •
-
Rsubread: Aligning multiple single and paired-end reads from multiple files (lanes)
written by schaeferchristina15 •
-
getting paired end datasets from SRA
written by from the mountains
-
Aligning Paired-end file in R/Rstudio using align function
written by student •
-
How get read group information from fastq files in SRA
written by zamani_2012 •
-
fastq-dump: failure to get two fastq files
written by msrk04011 •
-
How many paired FASTQ reads to input for alignment if 4 different SRA IDs belong to one sample-whol…
written by flight22121 •
-
BWA MEM alignment output of splited fastq files differ from the original(unsplit) fastq file
written by manojkumar_bhosale
You can easily split the files using
reformat.shfrom BBMap suite.You could have also used
--split-filesoption withfastq-dumpto get pair of files.@genomax just as an aside, would
reformat.shproduce an emptyR2.fqif the input is unpaired reads?Yes. If inputs can't be recognized by
/1 or /2(old Illumina data) or1:N: or 2:N:(current format) at the end of read headers.That's clear, thanks!
That is for interleaved files. Be sure that fastq-dump did not simply concatenate the mates into a single long read which it sometimes does. Confirm that you have one entry per mate otherwise alignment is going to be a mess.