Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Using phylogenetic tools in v5.0 #304

Open
jeffreybm opened this issue Oct 5, 2023 · 0 comments
Open

Using phylogenetic tools in v5.0 #304

jeffreybm opened this issue Oct 5, 2023 · 0 comments

Comments

@jeffreybm
Copy link

Hi Jody,

I am interested in using the phylogenetic tools you had added to version 5. I quite dont understand how to use the functionality. If I run:
tb-profiler profile
--no_trim
--read1 ${reads_dir}/${sample}.trim.1.fastq.gz --read2 ${reads_dir}/${sample}.trim.2.fastq.gz
--prefix ${sample}
--call_whole_genome
--update_phylo
--dir tb_profiler_out

I get an error, see attached SRR7516395.errlog.txt
SRR7516395.errlog.txt

If I run without the --update_phylo flag run completes and there is an empty phylo.pb.lock file in results/

I assume we would need at least 3 samples to run on for a tree to be made, so I guess I dont understand how we would run the first sample and then update this phylo file with subsequent samples?

Thanks for the help!

Brendan

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant