From 30af494bb190f044c2a522e1d72525b6824eb9ac Mon Sep 17 00:00:00 2001 From: Arthur Date: Tue, 3 Dec 2024 12:55:41 +0000 Subject: [PATCH] read targets as LongDNA --- src/utils.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/utils.jl b/src/utils.jl index a2e702f..412e6cb 100644 --- a/src/utils.jl +++ b/src/utils.jl @@ -3,7 +3,7 @@ #### Read in fasta files """ get_reference(ref_fasta) - + Takes a path to a reference fasta. Returns a `Vector{LongDNA{4}}` of reference sequences and a matrix taxonomic classifications. The reference fasta must be a DADA2-formatted reference database. See [here](https://benjjneb.github.io/dada2/training.html) for examples. @@ -34,7 +34,7 @@ function get_targets(seq_fasta) seqs = [] open(FASTA.Reader, seq_fasta) do reader for record in reader - push!(seqs,sequence(record)) + push!(seqs,LongDNA{4}sequence(record)) push!(ids,identifier(record)) end end