jtamames / SqueezeMeta

A complete pipeline for metagenomic analysis
GNU General Public License v3.0
365 stars 78 forks source link

stop in step10 #880

Open mingxues2024 opened 1 month ago

mingxues2024 commented 1 month ago

Dear sir, I met problem in step10 Mapping with Bowtie2 (Langmead and Salzberg 2012, Nat Methods 9(4), 357-9) Creating reference from contigs Working with sample 1: A7 Getting raw reads Aligning to reference with bowtie [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0011.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0007.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0009.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0010.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0008.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0005.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0002.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0003.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0006.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0004.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0001.bam" : File exists [E::hts_open_format] Failed to open file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0000.bam" : File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0000.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0001.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0002.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0003.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0004.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0005.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0006.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0007.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0008.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0009.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0010.bam": File exists samtools sort: failed to create temporary file "/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A7.bam.tmp.0011.bam": File exists Error running samtools at /lustre/BIF/nobackup/sun091/miniforge3/envs/SqueezeMeta/SqueezeMeta/scripts/10.mapsamples.pl line 272, line 0.000000. Stopping in STEP10 -> 10.mapsamples.pl. Program finished abnormally

If you don't know what went wrong or want further advice, please look for similar issues in https://github.com/jtamames/SqueezeMeta/issues Feel free to open a new issue if you don't find the answer there. Please add a brief description of the problem and upload the /lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/syslog file (zip it first) Died at /lustre/BIF/nobackup/sun091/miniconda3/envs/SqueezeMeta/bin/SqueezeMeta.pl line 941. Could you give me some suggestions?

syslog.zip

fpusan commented 1 month ago

It could be this bug in samtools (https://github.com/samtools/samtools/issues/1739) but it hasn't triggered for us so far so I am not sure. It could also be something related to writing permissions... did you perhaps run out of disk quota? Has it worked ok with the test dataset that we provide with the databases?

mingxues2024 commented 1 month ago

It could be this bug in samtools (samtools/samtools#1739) but it hasn't triggered for us so far so I am not sure. It could also be something related to writing permissions... did you perhaps run out of disk quota? Has it worked ok with the test dataset that we provide with the databases?

Thanks for your suggestions . I have run step10 successfully, but there is still some wrong with step14 [0 seconds]: STEP14 -> BINNING: 14.runbinning.pl [E::bgzf_read] Read block operation failed with error 4 after 0 of 4 bytes ERROR: fail to open index BAM file '/lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam/Q19.A8.bam' Traceback (most recent call last): File "/lustre/BIF/nobackup/sun091/miniforge3/envs/SqueezeMeta/SqueezeMeta/bin/CONCOCT-1.1.0/scripts/concoct_coverage_table.py", line 91, in generate_input_table(args.bedfile, args.bamfiles, samplenames=samplenames) File "/lustre/BIF/nobackup/sun091/miniforge3/envs/SqueezeMeta/SqueezeMeta/bin/CONCOCT-1.1.0/scripts/concoct_coverage_table.py", line 41, in generate_input_table sys.stderr.write(out) TypeError: write() argument must be str, not bytes Bus error sh: 1: cannot create /lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/intermediate/binners/concoct/concoct_int/clustering_merged.csv: Directory nonexistent wc: /lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/intermediate/binners/concoct/: Is a directory Could you give me some suggestions?

fpusan commented 1 month ago

What did you change in order to make step 10 work? Can you share the /lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/results/10.Q19.mappingstat file here?

mingxues2024 commented 3 weeks ago

Dear sir.

I think it was also the problem of server. And I have attached the file.

Best


发件人: Fernando Puente-Sánchez @.> 发送时间: 2024年8月15日 9:51 收件人: jtamames/SqueezeMeta @.> 抄送: Sun, Mingxue @.>; Author @.> 主题: Re: [jtamames/SqueezeMeta] stop in step10 (Issue #880)

What did you change in order to make step 10 work? Can you share the /lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/results/10.Q19.mappingstat file here?

― Reply to this email directly, view it on GitHubhttps://github.com/jtamames/SqueezeMeta/issues/880#issuecomment-2290838654, or unsubscribehttps://github.com/notifications/unsubscribe-auth/BFTDQ6KCHZZDESYCHKXHTVLZRRMZ7AVCNFSM6AAAAABMNSSHECVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEOJQHAZTQNRVGQ. You are receiving this because you authored the thread.

fpusan commented 3 weeks ago

I did not receive the file. Please attach it in GitHub.

mingxues2024 commented 3 weeks ago

image

fpusan commented 3 weeks ago

I see, mapping seems to have worked fine, but it didn't generate the sorted BAM files. Can you post here the output of ls /lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/data/bam?

Probably you will need to repeat step 10 to make sure the BAM files get generated properly with

10.mapsamples.pl /lustre/BIF/nobackup/sun091/metagenome/coassembly_per_accession/Q19/