Dear RSubread developers and users:
We (my group) appreciate the Rsubread aligner! We have found it to be very useful and its great being able to do all of our genome analysis, including alignment, completely in R.
We are trying to use/optimize Rsubread to conduct some microbiome data analysis. Naturally, our microbial reference libraries are >>4GB. We are getting the following error:
"ERROR: The chromosome data contains too many bases. The size of the input FASTA files should be less than 4G Bytes"
So I am assuming that I need to break my reference genome library into chunks <4GB, do the alignments, and then merge the .bam files? We have been using Bowtie2, and have to do the same thing. Just asking if anyone knows an easier or more concise way to do it. If not, we will move forward with the separate libraries/merging .bam approach.