How to calculate GC/AT content in R language
What this article shares with you is about how to calculate GC/AT content in R language. The editor thinks it is very practical, so I share it with you. I hope you can get something after reading this article.
Today I would like to introduce to you a package that can handle FASTA files-Biostrings. This package mainly deals with some sequence information of the genome, including sequence translation, DNA/RNA interconversion, statistics of each base content, triplet content. These can be solved by one-line command. Today, let's first teach you how to calculate the GC/AT content.
The first is the installation, with the following code:
Source ("http://bioconductor.org/biocLite.R")biocLite("Biostrings")"
You need to wait patiently for a few minutes after entering the code.
After installation, you only need to type a few lines of code to achieve GC/AT content visualization.
# sequence file storage path filepath