This package can be used to compare indiviual genes of the two genomes. This package makes use of coding sequences format that can be found in Genbank. Right now, the package allows pairwise comparison of the genes of two genomes. You find detail manual in the github wiki section. You will need FASTA Nuceotide sequence files of the two genomes. These can be obtained from the Genbank. If you click on Send to button on Genbank, then you can select "Coding Sequences" and then FASTA Nucleotide or FASTA Protein sequences. Then you can run pairwise_genes_comparison() funtion in R. You have to provide path for two files and then you can run pairwise_genes_comaprison function on those two files. It will create a nice summary table output. It will also save the table as a csv file called "summary_table,csv" in the current folder. It uses DECIPHER package for pairwise alignment and to generate similarity matrix.
filepath1= path to the first genome file
filepath2= path to the second genome file
Dependency : DECIPHER, tidyverse,formattable, seqinr
Make sure you have all the depedencies installed.
Lok Raj Joshi
1 2 3 4 5 6 7
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.