reducing memory consumption of circoletto
Hi,
- I splitted my query genome in 200 parts with pyfasta and I ran each of them separetly with blastall: blastall -p blastn -i ${i} -d ${db} -F F -e 1e-10 -E -1 -v 20000 -b 20000 -a 1 > output/${output}.blast
- I merged all blastall results and the output file is 8.8 Tb big: cat BAFB.contigs.*.fasta.clean.blast > BAFB.contigs.all.fasta.clean.blast
- Unfortunately, 50gb of RAM was not enough to run circoletto.pl perl /apps/Circoletto/circoletto.pl --blastout BAFB.contigs.all.fasta.clean.blast --cpus 8 --out_dir output
My blastall database is based on a reference genome with 10 chromosomes but my assembly has 16.000 contigs.
How can I reduce memory consumption of circoletto.pl ?
Thank you in advance.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Error Scaffolding with RagTag
written by Bipin •
-
BWA alignment failing after adding nextflow.config file?
written by Eliveri •
-
Question about blastx output
written by yinbinqiu •
-
Setting up local blast in webserver
written by sharmatina189059 •
-
OrthoMCL issue with BLAST and mySQL
written by merwright •
-
extract chromosome names from blast output
written by Ric
-
How to reduce the size of blastall output
written by Ric
-
Convert Blastall cmd to Blast+
written by cowmoomoo •
-
To to set the number of HSP in blastn tabular output
written by Mingkun •
-
What might explain the "[blastall] ERROR: Arguments must start with '-'" error in my code?
written by joreamayarom
