get_seqs: Gets sequences from bam file

View source: R/metascope_blast.R

get_seqsR Documentation

Gets sequences from bam file

Description

Returns fasta sequences from a bam file with a given taxonomy ID

Usage

get_seqs(id, bam_file, n = 10, bam_seqs)

Arguments

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

Value

Biostrings format sequences


compbiomed/MetaScope documentation built on Sept. 23, 2024, 5:25 a.m.