Skip to content

running sintax classification and utax  #7

@Liakumat

Description

@Liakumat

I recently emailed you, please see below query:

Issue No// 1
I am using linux on windows and have tried to run utax in ubuntu but I have a 64 bit see erro below:

./usearch8.0.1623_i86linux32/usearch8.0.1623_i86linux32 -utax PoJ100_S196_L001.q20.fq -db ../analyses/viridiplantae_all_2014.utax.udb -utax_rawscore -tt ../analyses/viridiplantae_all.utax.tax -utaxout filename1.utax
usearch v8.0.1623_i86linux32, 4.0Gb RAM (7.8Gb total), 16 cores
(C) Copyright 2013-15 Robert C. Edgar, all rights reserved.
http://drive5.com/usearch

License: personal use, non-transferrable

00:01 183Mb  100.0% Rows
00:01 184Mb Reading pointers...done.
00:01 218Mb Reading seqs...

./usearch8.0.1623_i86linux32/usearch8.0.1623_i86linux32 -utax PoJ100_S196_L001.q20.fq -db ../analyses/viridiplantae_all_2014.utax.udb -utax_rawscore -tt ../analyses/viridiplantae_all.utax.tax -utaxout filename1.utax

---Fatal error---
ReadStdioFile failed, attempted 55187803 bytes, read 51307735 bytes, errno=2

Issue No// 2:

for f in ../raw/filtered/*.fq; do out=$(basename $f .fq); usearch -sintax $f           -db precomputed/viridiplantae_all_2014.sintax.fa           -tabbedout $out.sintax.txt           -strand both           -sintax_cutoff 0.90; done

This produced output files the txt output files- *.sintax.txt, however to use the utax count.pl below i renamed the files to *.utax

for i in ../*.utax; do     BASE=$(basename $i .utax);     perl ../../code/count_taxa_utax.pl --in $i --cutoff 20 > $BASE.count; done

I am getteting the following error

x.pl line 139, <IN> line 18896.
Argument "k:Viridiplantae_33090(1.0000" isn't numeric in numeric ge (>=) at ../../code/count_taxa_utax.pl line 139, <IN> line 18897.
Argument "k:Viridiplantae_33090(1.0000" isn't numeric in numeric ge (>=) at ../../code/count_taxa_utax.pl line 139, <IN> line 18898.
Argument "k:Viridiplantae_33090(1.0000" isn't numeric in numeric ge (>=) at ../../code/count_taxa_utax.pl line 139, <IN> line 18899.
Argument "k:Viridiplantae_33090(1.0000" isn't numeric in numeric ge (>=) at ../../code/count_taxa_utax.pl line 139, <IN> line 18900.

Any code to use sintax.txt output file to count reads and follow the tutorial

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions