View source: R/metascope_blast.R
get_seqs | R Documentation |
Returns fasta sequences from a bam file with a given taxonomy ID
get_seqs(id, bam_file, n = 10, bam_seqs)
id |
Taxonomy ID of genome to get sequences from |
bam_file |
A sorted bam file and index file, loaded with Rsamtools::bamFile |
n |
Number of sequences to retrieve |
bam_seqs |
A list of the sequence IDs from the bam file |
Biostrings format sequences
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.