Closed Michaelijesse closed 10 months ago
Could you please sort this out. why its happening. Even this is sample dataset
[0 seconds]: STEP10 -> MAPPING READS: 10.mapsamples.pl Reading samples from /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/00.Samples.samples Metagenomes found: 2 Mapping with Bowtie2 (Langmead and Salzberg 2012, Nat Methods 9(4), 357-9) Creating reference from contigs Working with sample 1: SRR1927149 Getting raw reads Aligning to reference with bowtie BAM file already found in /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.bam, skipping Calculating contig coverage Reading contig length from /home/ec2-user/han/20230828_WMS_test_2/output/Samples/intermediate/01.Samples.lon Illegal division by zero at /home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/scripts/10.mapsamples.pl line 445. 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 /home/ec2-user/han/20230828_WMS_test_2/output/Samples/syslog file (zip it first)
What's your OS and version?
Can you share the /home/ec2-user/han/20230828_WMS_test_2/output/Samples/syslog
file with us here?
I am working in ec2 server. Instance type : M5.2xlarge syslog.zip
So what OS and version are you running there?
Can you run manually the following command and let me know what error is raised?
/home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/bowtie2/bowtie2 -x /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/Samples.bowtie -1 /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_1.gz -2 /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_2.gz --quiet -p 12 -S /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.sam --very-sensitive-local
Sure I will run and post the output here. Thank you.
I tried it successfully. No Error. It took approx 30mins. SAM file created in the bam folder.
$ cd /home/ec2-user/han/20230828_WMS_test_2/ $ /home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/bowtie2/bowtie2 -x /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/Samples.bowtie -1 /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_1.gz -2 /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_2.gz --quiet -p 12 -S /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.sam --very-sensitive-local $
Ok, then can you run
/home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/samtools sort /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.sam -o /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.bam -@ 12
/home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/samtools index /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.bam -@ 12
It should fail in one of these two
OS : amazon-linux (centos)
I have chosen this tool for metagenomic assembly for multiple good reasons. but due to errors at some steps, I cant able to come to conclusion. Could you please sort this out.?
Lines 445,446,447 my $mapperc=($mappedreads/$totalreadcount)*100; if($mapperc<50) { $warnmes=1; } printf outfile1 "$thissample\t$totalreadcount\t$mappedreads\t%.2f\t$totalreadlength\n",$mapperc; #-- Mapping statistics
Hi again
Did you run these two specific commands as I suggested above?
/home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/samtools sort /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.sam -o /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.bam -@ 12
/home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/samtools index /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.bam -@ 12
If so, which one of those failed?
Also was this installed through conda? In the screenshot above I can see that you are in the (base)
environment, which suggests that you installed it manually.
SqueezeMeta has not been tested in amazon linux so certain binaries may not work properly there. I pretty much suspect that samtools is failing in your case. If we confirm that (see my post above to test it) we can start testing for solutions.
Helllo, Output from Bowtie2.
Ok, I have uploaded a dev version that may fix this issue
Try
mamba create -n SqueezeMeta_dev -c conda-forge -c bioconda -c anaconda -c fpusan squeezemeta-dev=1.6.3.beta1 --no-channel-priority
And then
conda activate SqueezeMeta_dev
and run the test data again using the SqueezeMeta.pl
executable present in your PATH
Note that you will need to start a new project, restarting won't work in this case since the run was bugged
And of course after making the new install you'll need to link it to the old databases using configure_nodb.pl
(check the ReadMe for details)
configure_nodb.pl
will also trigger a series of tests, including some tests to check that certain binaries can be executed in your system. Can you share with me the output of that after running it on the new SQM install?
Hello,
Running configure_nodb.pl, failed mothur configuration. Even I installed Mothur through conda but Test failed.
Installed mothur.
That's encouraging. mothur is not needed for the main SqueezeMeta pipeline (though I see you made It work regardless) Then can you try running the pipeline again using the conda version?
Hello, Thank you for your guidance and support. The whole process completed successfully.
So closing the issue
Great! We will release a new version incorporating these changes shortly
This is SYSLOG out
[17 hours, 18 minutes, 44 seconds]: STEP10 -> 10.mapsamples.pl Creating Bowtie reference: python3 /home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/bowtie2/bowtie2-build --quiet /home/ec2-user/han/20230828_WMS_test_2/output/Samples/results/01.Samples.fasta /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/Samples.bowtie Getting raw reads for SRR1927149: cp /home/ec2-user/han/20230828_WMS_test_2/input/raw/SRR1927149_1.fastq.gz /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_1.gz; cp /home/ec2-user/han/20230828_WMS_test_2/input/raw/SRR1927149_2.fastq.gz /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_2.gz; Aligning with bowtie: /home/ec2-user/han/20230828_WMS_test_2/SqueezeMeta/bin/bowtie2/bowtie2 -x /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/Samples.bowtie -1 /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_1.gz -2 /home/ec2-user/han/20230828_WMS_test_2/output/Samples/temp/Samples.SRR1927149.current_2.gz --quiet -p 12 -S /home/ec2-user/han/20230828_WMS_test_2/output/Samples/data/bam/Samples.SRR1927149.sam --very-sensitive-local Stopping in STEP10 -> 10.mapsamples.pl. Program finished abnormally