Rsubread: Aligning multiple single and paired-end reads from multiple files (lanes)
Entering edit mode
4 months ago


I am new to bioinformatics and looking for some help. I have 27 files from an Illumina output. There are 4 paired end and 23 single read files. I am trying to align them using Rsubread in R directly.

I want to align all of the files at once using a code similar to this:

align(index = "reference_index", readfile1="/path to file/forward_reads_list",readfile2="/path to_file/reverse_reads_list", input_format="FASTQ", output_file="rsubread.bam", type = "dna)

Can I use a single function to align both the single read and paired-end reads? Or do I have to do this separately? Further, if I do have to do two separate functions, how do I merge them later on?

The ultimate goal is to get FASTQ files to align to BAM then pull out variants (using VariantAnnotation) and detect similar SNPs between this dataset and another dataset from DArT sequencing.

Thanks in advance!

pair-end multiple alignment Rsubread • 276 views
Entering edit mode

Did you mean you have 27 samples, 4 of them are paired end and 23 are single end?


Login before adding your answer.

Traffic: 586 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6