This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error with gzipped FASTA file in mash.

Hi i am trying to run mashtree to calculate WGS distance and i encountered this error by using following commands. If anyone know anything about this error, please let me know.

## mashtree_bootstrap.pl --reps 100 --numcpus 30 *.fna.gz -- --min-depth 0 > mashtree.bootstrap.dnd##

Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. Use of uninitialized value $dist in concatenation (.) or string at /home/bvs/anaconda3/envs/myenv/lib/site_perl/5.26.2/Mashtree/Db.pm line 263. mashtree_bootstrap.pl: main: Running mashtree single-threaded among 100 replicates using 30 total threads. mashtree_bootstrap.pl: mashtreeRepWorker(TID1): Rep 0 mashtree_bootstrap.pl: mashtreeRepWorker(TID2): Rep 1 mashtree_bootstrap.pl: mashtreeRepWorker(TID3): Rep 2 mashtree_bootstrap.pl: mashtreeRepWorker(TID4): Rep 3 mashtree_bootstrap.pl: mashtreeRepWorker(TID5): Rep 4 mashtree_bootstrap.pl: mashtreeRepWorker(TID6): Rep 5 mashtree_bootstrap.pl: mashtreeRepWorker(TID7): Rep 6 mashtree_bootstrap.pl: mashtreeRepWorker(TID8): Rep 7 mashtree_bootstrap.pl: mashtreeRepWorker(TID9): Rep 8 mashtree_bootstrap.pl: mashtreeRepWorker(TID10): Rep 9 mashtree_bootstrap.pl: mashtreeRepWorker(TID11): Rep 10 Thread 1 terminated abnormally: mashtree_bootstrap.pl: main::mashtreeRepWorker: ERROR on rep 0 mashtree: main: Found mash version 2 - /home/bvs/neelam/mash/mash mashtree: main: Temporary directory will be /tmp/MASHTREE_BOOTSTRAP.il93zT/0 mashtree: main: mashtree on 359 files mashtree: mashSketch(TID1): This thread will work on 359 sketches mashtree: mashSketch(TID1): Working on file 1 out of 359

gzip: P_aeruginosa_12-4-4_59.fna.gz: not in gzip format Thread 1 terminated abnormally: mashtree: main::mashSketch: ERROR with gzip -cd P_aeruginosa_12-4-4_59.fna.gz Stopped at /home/bvs/anaconda3/envs/myenv/bin/mashtree line 236. mashtree: mashDistance: Waiting to join thread (1/1, TID2) mashtree: mashDistance: Databasing distances (1/1, TID2) mashtree: mashDistance: Converting to phylip format into /tmp/MASHTREE_BOOTSTRAP.il93zT/0/distances.phylip mashtree: Mashtree::createTreeFromPhylip: Can't call method "as_text" on an undefined value Stopped

Thank you!

mash
• 1,551 views • link
written by Neel • Neel

Please type file P_aeruginosa_12-4-4_59.fna.gz in gnu-linux OS and post the output here.

• 0 views • link
written by cpad0112 cpad0112

Thank you so much for your time...Now it solve.

• 0 views • link
written by Neel • Neel

1 answer

I have highlighted the significant portion of the error message. It is likely that you simply have to download the fna.gz file. It is likely damaged or not in gzip format. Try to download it again using wget or curl. If you are unsure of its format try to gzip -t it manually.

• 0 views • link
written by Michael Dondrup Michael Dondrup

Thank you it was very helpful. Now it solve.

• 0 views • link
written by Neel • Neel

Log in to answer this question.

More posts like this
  • ValueError: not enough values to unpack (expected 2, got 1)
    written by Neel • Neel
    Hi i am using Clinker and it gives following error- clinker P_aeruginosa_AR_0110.gbk P_aeruginosa_AR_0111.gbk P_aeruginosa_AR_455.gbk P_aeruginosa_ATCC_15692.gbk -r 4_dataforcli.txt -o out_4 -i 70 -p 4_plot.html [15:35:03] INFO …
  • BRAKER - Genome annotation, software error ?
    written by sunnykevin97 sunnykevin97
    Hi, Try to run BRAKER after installing all the dependencies, unable to run couldn't really understand the error. In-addition I installed the `Scalar::Util::Numeric` perl module. …
  • NCBI annotated sequence use as input of Roary
    written by Manoj • Manoj
    NCBI sequence (.fna) and gene annotation file (.gff) downloads using curl --remote-name --remote-time ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/192/045/GCF_000192045.2_ASM19204v3/GCF_000192045.2_ASM19204v3_genomic.fna.gz These two files was merged to create a gff3 file using …
  • Issues installing vep
    written by pmw • pmw
    I am having issues with installing ensembl's vep on the HPC. It requires perl dependencies which I have successfully installed most of via anaconda. However …
  • miRDeep2.pl error in make_html2.pl
    written by tikshyadav19 • tikshyadav19
    I have installed cpan install PDF::API2 And still got error Can't locate PDF/API2.pm in @INC (@INC contains: /home/nbri/Tikshana/Ath/small_RNA/mirDeep/mirdeep2/lib/perl5 /home/nbri/vaishali/TOMATO_SRNA/new_analysis/mirdeep2-master/lib/perl5 /usr/local/lib64/perl5 /usr/local/share/perl5 /usr/lib64 /perl5/vendor_perl /usr/share/perl5/vendor_perl /usr/lib64/perl5 …
  • VEP installation problem
    written by jimmy0958073736 • jimmy0958073736
    Hi,everyone I tried to install one of package, but something in my environment. Could any one give me a help I download the VEP package …
  • How to install Bioperl on Ubuntu 18.x
    written by nadarajan_v • nadarajan_v
    I am trying to install bioperl on my Ubuntu 18 machine. I ran the following lines of code as suggested by https://zoomadmin.com/HowToInstall/UbuntuPackage/bioperl sudo apt-get update …
  • Error while running vcf2maf tool
    written by Jordan Jordan
    Hi, I'm running vcf2maf.pl on about 100 VCF files. But some files (about 10) do not report any MAF, and they report the below error: …
  • Apollo installtion for maker
    written by ashaneev07 • ashaneev07
    Hiii.. As i'm installing maker genome annotation pipeline, when it comes with building dependencies it shows the following errors while the installation of Apollo. have …
  • Help in NGS QC Toolkit
    written by caelyn5201314 • caelyn5201314
    Hi All! I was use NGS QC Toolkit IlluQC.pl),and it keep throw the error: ``` Use of uninitialized value $avgQRangeF2 in split at /leofs/zhangz_group/xial/NGSQCToolkit_v2.3.3/QC/lib/html.pl line …