I am pretty new to R. So I apologize for asking maybe a very basic question. Let's say I have a fasta file with sequence below:
>sequence_1 ACCTGC--A >sequence_2 ACC-GCTTA >sequence_3 ACCTGCTTA
Is there a function or method to calculate entropy for each column in R? I did some research and found a function that calculate entropy for Amino Acids but not for DNA nucleotides.
I would like to have the output be stored in a vector or in a list.
Thanks for the help in advance!