juliema / aTRAM

BSD 3-Clause "New" or "Revised" License
35 stars 14 forks source link

stitcher does not recognise contig file #312

Open Upasana1991 opened 2 years ago

Upasana1991 commented 2 years ago

I have list of contig file produced in this format S0057.S0057_OG0007378.fasta.all_contigs.fasta S0057.S0057_OG0007378.fasta.filtered_contigs.fasta the target file is named OG0007378.fasta which has a protein sequence ">ADIR005593.1"

I run this command python atram_stitcher.py \ --assemblies-dir=~/scratch/aTRAM-2.4.3/S0057.S0057_OG0007378.fasta.filtered_contigs.fasta\ --reference-genes=~/scratch/aTRAM-2.4.3/OG0007378.fasta \ --taxa=~/scratch/aTRAM-2.4.3/Taxa.txt \ --temp-dir=~/scratch/aTRAM-2.4.3/directory \ --output-prefix=Dirus

And I get this error

10902) 2022-05-27 20:37:08.436959 INFO : aTRAM version: v2.4.3 ( 10902) 2022-05-27 20:37:08.436989 INFO : Python version: 3.9.7 | packaged by conda-forge | (default, Sep 29 2021, 19:23:11) [GCC 9.4.0] ( 10902) 2022-05-27 20:37:08.437008 INFO : ~/scratch/aTRAM-2.4.3/atram_stitcher.py --assemblies-dir=~/scratch/aTRAM-2.4.3/S0057_db/ --reference-genes=~/scratch/aTRAM-2.4.3/ADIR/OG0007378.fasta --taxa=~/scratch/aTRAM-2.4.3/Taxa.txt --temp-dir=~/scratch/aTRAM-2.4.3/directory --output-prefix=Dirus ( 10902) 2022-05-27 20:37:08.575545 INFO : Preparing taxa ( 10902) 2022-05-27 20:37:08.577960 INFO : Preparing reference genes: ~/scratch/aTRAM-2.4.3/ADIR/OG0007378.fasta ( 10902) 2022-05-27 20:37:08.768478 INFO : Preparing reference gene files for exonerate ( 10902) 2022-05-27 20:37:08.772011 INFO : First contig insert: ~/scratch/aTRAM-2.4.3/S0057_db/ ( 10902) 2022-05-27 20:43:22.439124 INFO : Write contig files ( 10902) 2022-05-27 20:43:22.440153 INFO : Exonerate run for: ADIR005593.1 ( 10902) 2022-05-27 20:43:22.444345 FATAL: There were no hits. Please make sure your contig file names contain both a reference gene name and taxon name in them.

the taxon name was S0057

Could you help me figure out what's the error or give an example of the name of the files (assembly file, taxon name, protein file) I should have before I run stitcher.py