Read Biological Sequences File Formats with R.
Several (standardized) file formats have been published and are actively used
to store, exchange, and analyze biological sequences. The seqreadr
package
aims to provide read/write functions to work with this variety of sequence file formats.
To avoid re-inventing the wheel the seqreadr
package relies on existing interface functions implemented in a variety of R packages. However, additional functionality for
sequence quality control is provided by the functions implemented in seqreadr
.
Furthermore, seqreadr
aims to provide a complete collection of sequence file formats and is based on a unifying terminology and nomenclature to avoid learning different package styles etc.
Use Cases for each Sequence File Format:
- FastA/FastQ
- SAM/BAM
- GFF/GTF
- BED/Big Bed
- PSL
- Pairwise Interaction File Format
- Variation File Format
- WIG File Format
The current status of the package as well as a detailed history of the functionality of each version of seqreadr
can be found in the NEWS section.
# install developer version of seqreadr
library(devtools)
install_github("HajkD/seqreadr", build_vignettes = TRUE, dependencies = TRUE)
I would be very happy to learn more about potential improvements of the concepts and functions provided in this package.
Furthermore, in case you find some bugs or need additional (more flexible) functionality of parts of this package, please let me know:
https://github.com/HajkD/seqreadr/issues
read.genome()
: Read the genome of a given organism (DNA)read.proteome()
: Read the proteome of a given organism (Protein)read.cds()
: Read the CDS of a given organism (DNA)
read.gtf()
: Read GTF (General Transfer Format) Fileswrite.gtf()
: Write GTF files to hard drive
read.bed()
: Read BED fileswrite.bed()
: Write BED files to hard drive