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