alignMitogenomes: alignMitogenomes

In chutter/MitoCap: Package For extracting, assembling, annotating and aligning mitochondrial genomes from high-throughput sequencing data

alignMitogenomes

Description

Function for batch trimming a folder of alignments, with the various trimming functions available to select from

Usage

alignMitogenomes(
  alignment.folder = NULL,
  genbank.file = NULL,
  draft.contigs = "draftContigs",
  output.dir = "Genomes",
  dataset.name = "untrimmed",
  overwrite = FALSE
)

Arguments

alignment.folder	path to a folder of sequence alignments in phylip format.
genbank.file	available input alignment formats: fasta or phylip
draft.contigs	contigs are added into existing alignment if algorithm is "add"
output.dir	available output formats: phylip
dataset.name	remove samples too divergent from consensus, values 0-1 for proportion similar sites
overwrite	TRUE to supress mafft screen output

Value

an alignment of provided sequences in DNAStringSet format. Also can save alignment as a file with save.name

Examples

your.tree = ape::read.tree(file = "file-path-to-tree.tre")
astral.data = astralPlane(astral.tree = your.tree,
                          outgroups = c("species_one", "species_two"),
                          tip.length = 1)

chutter/MitoCap documentation built on July 17, 2025, 12:04 a.m.