Hi,--I am having a problem installing ABySS on CentOS. When I execute the 'make' command I got the following error:Making all in DataBase
make[2]: Entering directory `/Shared/002-Software/programs/abyss-1.9.0/DataBase'
g++ -Wall -Wextra -g -O2 -o abyss-db-csv abyss_db_csv-DB.o abyss_db_csv-db-csv.o -lsqlite3 -ldl -lm
/bin/ld: cannot find -lsqlite3
collect2: error: ld returned 1 exit status
make[2]: *** [abyss-db-csv] Error 1
make[2]: Leaving directory `/Shared/002-Software/programs/abyss-1.9.0/DataBase'
make[1]: *** [all-recursive] Error 1
make[1]: Leaving directory `/Shared/002-Software/programs/abyss-1.9.0'
make: *** [all] Error 2
I have sqlite3 installed and in my path. Any help will be very appreciated.Thank youRiccardo
You received this message because you are subscribed to the Google Groups "ABySS" group.
To unsubscribe from this group and stop receiving emails from it, send an email to abyss-users...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
./configure --with-sqlite=/home/software/lib/
configure: WARNING: boost/functional/hash.hpp: present but cannot be compiledconfigure: WARNING: boost/functional/hash.hpp: check for missing prerequisite headers?
configure: WARNING: boost/functional/hash.hpp: see the Autoconf documentation
configure: WARNING: boost/functional/hash.hpp: section "Present But Cannot Be Compiled"
configure: WARNING: boost/functional/hash.hpp: proceeding with the compiler's result
configure: WARNING: ## ----------------------------------- ##
configure: WARNING: ## Report this to abyss...@bcgsc.ca ##
configure: WARNING: ## ----------------------------------- ##
Uncompress.cpp: In function 'bool uncompress_init()':
Uncompress.cpp:220:9: error: 'HAVE_LIBDL' was not declared in this scope
return HAVE_LIBDL;
^
Uncompress.cpp:221:1: warning: control reaches end of non-void function [-Wreturn-type]
}
^
make[2]: *** [libcommon_a-Uncompress.o] Error 1
make[2]: Leaving directory `/media/sequentia/NAS/software/001-Software/abyss-1.9.0/Common'
make[1]: *** [all-recursive] Error 1
make[1]: Leaving directory `/media/sequentia/NAS/software/001-Software/abyss-1.9.0'
make: *** [all] Error 2
Could you please help me?
Thank you very much
Riccardo
Hi Ben,finally I managed to install Abyss and then TransAbyss. Unfortunately now I have another problem, I am running some assembly tests and when I execute the command I get this message and then an error:make: Entering directory `/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/K21_prova'ABYSS -k21 -q3 -e2 -E0 -c2 --SS --coverage-hist=coverage.hist -s AbyssK21-bubbles.fa -o AbyssK21-1.fa /Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R1.fastq /Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R2.fastqABySS 1.9.0ABYSS -k21 -q3 -e2 -E0 -c2 --SS --coverage-hist=coverage.hist -s AbyssK21-bubbles.fa -o AbyssK21-1.fa /Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R1.fastq /Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R2.fastqReading `/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R1.fastq'...`/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R1.fastq': reversed 1000000 reads`/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R1.fastq': discarded 2 reads containing non-ACGT charactersReading `/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R2.fastq'...`/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/prova_R2.fastq': discarded 68 reads containing non-ACGT charactersLoaded 47800009 k-merMinimum k-mer coverage is 44Using a coverage threshold of 2...The median k-mer coverage is 3The reconstruction is 30231282The k-mer coverage threshold is 1.73205Generating adjacencyAdded 95526490 edges.Eroding tipsEroded 13008220 tips.Eroded 0 tips.Pruning tips shorter than 1 bp...Pruned 1811 k-mer in 1811 tips.Pruning tips shorter than 2 bp...Pruned 3280 k-mer in 2118 tips.Pruning tips shorter than 4 bp...Pruned 11595 k-mer in 4283 tips.Pruning tips shorter than 8 bp...Pruned 41623 k-mer in 8550 tips.Pruning tips shorter than 16 bp...Pruned 138839 k-mer in 15622 tips.Pruning tips shorter than 21 bp...Pruned 107598 k-mer in 8591 tips.Pruning tips shorter than 21 bp...Pruned 5 k-mer in 2 tips.Pruning tips shorter than 21 bp...Pruned 40977 tips in 7 rounds.Marked 714517 edges of 348850 ambiguous vertices.Removing low-coverage contigs (mean k-mer coverage < 2)Found 34478236 k-mer in 596193 contigs before removing low-coverage contigs.Removed 4643987 k-mer in 143655 low-coverage contigs.Split 258356 ambigiuous branches.Eroding tipsEroded 20080 tips.Eroded 0 tips.Pruning tips shorter than 1 bp...Pruned 1470 k-mer in 1470 tips.Pruning tips shorter than 2 bp...Pruned 2065 k-mer in 1648 tips.Pruning tips shorter than 4 bp...Pruned 5889 k-mer in 2827 tips.Pruning tips shorter than 8 bp...Pruned 16319 k-mer in 4224 tips.Pruning tips shorter than 16 bp...Pruned 47417 k-mer in 6326 tips.Pruning tips shorter than 21 bp...Pruned 42081 k-mer in 3326 tips.Pruning tips shorter than 21 bp...Pruned 116 k-mer in 17 tips.Pruning tips shorter than 21 bp...Pruned 19838 tips in 7 rounds.Popping bubblesRemoved 11695 bubbles.Removed 11695 bubblesMarked 207407 edges of 99661 ambiguous vertices.Left 12898 unassembled k-mer in circular contigs.Assembled 29412231 k-mer in 252740 contigs.Removed 18092395 k-mer.The signal-to-noise ratio (SNR) is 2.15372 dB.make: *** No rule to make target `AbyssK21-1.adj'. Stop.make: Leaving directory `/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/K21_prova'Elapsed time: 0 h 5 m 42 sERROR: CMD ended with status code 2I understand that the file .adj should be generated by the command AdjList. I tried to to use it separately and it works.AdjList -k 21 AbyssK21-1.fa > AbyssK21-1.adjThis try has been done with a small dataset with 1M of pe reads. I started with a full dataset with 35M pe reads and I had the same error.Could you please help me to solve this problem?Thank youRiccardo--
Dr. Riccardo Aiese Cigliano
Scientific Director and Co-funderSequentia Biotech SL
Parc Cientific de Barcelona
Torre R 02A6/A7
c/ Baldiri Reixac 4-8
08028 Barcelona
Cel. +34 633042534
raiesecigliano@sequentiabiotech.com
www.sequentiabiotech.com
This message is intended exclusively for its addressee and may contain information considered CONFIDENTIAL and protected by professional privilege. If you are not the intended recipient you are hereby notified that any dissemination, copy or disclosure of this communication is strictly prohibited by law. If this message has been mistakenly received, please immediately notify us via e-mail and delete it.
n n:200 L50 min N80 N50 N20 max sum
252740 50500 16292 200 262 389 625 2358 18.76e6 K21_prova/AbyssK21-1.fa
Then I checked the fasta file with an head:
>0 143 322
CTCTGTTGTTGGCCATATATTGAGTACTCTGTTATAGCCTTAGAGGAACTGCACTGCCTTCTCTCACTCTGTTAGCCAACAATTTAGACCATCTGTTTGTTCAAACAAGGCTAAACTTACAAAATGTGCTATAGCATGTGCTA
>1 274 1050
GTGGAGATCATTGGGAAGAATGCAGTGGTAATAGGGAGAAGCAAGATTGCTGGGTTATCCACTTCTTTGCTCTTGCAGTCCTGGACAGAGGCACCACGCGACGGTTAGCACTCTCCATTCGTTCACCAACAACCCAGAACAAATTACTCGTCGAGCTGATATTGTTGTATCAGATGTTGGCATTCCAAATAGAGTCGGATGCAATTGGCTGATGCCAGGGGCAGTTGTAGTCGACATGGGGACAAATTCAGTTAAGGATCCCAATAGTCCCCGA
And then with a tail:
>252738 44 48
ATGATGGGCGCCTAATAATTTGTGACGAAGCAGAATCGTACAAA
>252739 75 199
TCCAACCGCTACTACGACTTAAAGACCTTGCTCAGCGGGGCAATGACCAACCAGTACACGTGTCTTGATGGGTTT
I am not an expert but I don't see anything strange. Finally I repeated only the abyss-pe command as follows (with the verbose mode):
abyss-pe k=21 name=K21_prova lib=pe in='prova_R1.fastq prova_R2.fastq' e=2 E=0 c=2 SS=--SS j=16 v=-v
make: Entering directory `/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/K21_prova'
make: `AbyssK21-1.fa' is up to date.
make: Nothing to be done for `AbyssK21-1.adj'.
make: Leaving directory `/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/K21_prova'
CHECKPOINT: De Bruijn graph assembly completed.
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
Assembly stopped at stage 'dbg'
Total wallclock run time: 0 h 0 m 0 s
Traceback (most recent call last):
File "/Shared/002-Software/programs/transabyss-1.5.3/transabyss", line 748, in <module>
__main__()
File "/Shared/002-Software/programs/transabyss-1.5.3/transabyss", line 553, in __main__
if perform_pe_assembly and not has_edges(adj3):
File "/Shared/002-Software/programs/transabyss-1.5.3/utilities/adj_utils.py", line 1017, in has_edges
with open(adj_file, 'r') as fh:
IOError: [Errno 2] No such file or directory: '/Synology/002-Projects/084-Cherry/Trimmed_Reads/Transabyss/K21_prova/AbyssK21-3.adj'
Hi Riccardo,
www.sequentiabiotech.com
This message is intended exclusively for its addressee and may contain information considered CONFIDENTIAL and protected by professional privilege. If you are not the intended recipient you are hereby notified that any dissemination, copy or disclosure of this communication is strictly prohibited by law. If this message has been mistakenly received, please immediately notify us via e-mail and delete it.
Hi Ben,I made a new installation of Abyss and Transabyss and I started over using the test dataset that is in the installation folder of Transabyss. So first I ran the the assembly.sh script that produced again the error:make: *** No rule to make target `test-1.adj'. Stop.make: Leaving directory `/Shared/002-Software/programs/transabyss-1.5.3/sample_dataset/test/assembly'
ERROR: CMD ended with status code 2
I executed the abyss-pe command: abyss-pe k=32 in='reads/rnaseq_1.fq.gz reads/rnaseq_2.fq.gz' name=test v=-v and I got this log:ABYSS -k32 -q3 -v --coverage-hist=coverage.hist -s test-bubbles.fa -o test-1.fa reads/rnaseq_1.fq.gz reads/rnaseq_2.fq.gzABySS 1.9.0ABYSS -k32 -q3 -v --coverage-hist=coverage.hist -s test-bubbles.fa -o test-1.fa reads/rnaseq_1.fq.gz reads/rnaseq_2.fq.gzReading `reads/rnaseq_1.fq.gz'...Read 10000 reads. Hash load: 5006 / 1073741824 = 4.66e-06 using 359 MBReading `reads/rnaseq_2.fq.gz'...Read 10000 reads. Hash load: 6361 / 1073741824 = 5.92e-06 using 359 MBLoaded 6361 k-merHash load: 6361 / 16384 = 0.388 using 852 kBMinimum k-mer coverage is 5Coverage: 5 Reconstruction: 3563Coverage: 14.8 Reconstruction: 3499Coverage: 14.8 Reconstruction: 3499Using a coverage threshold of 15...The median k-mer coverage is 219The reconstruction is 3499The k-mer coverage threshold is 14.8Setting parameter e (erode) to 15Setting parameter E (erodeStrand) to 1Setting parameter c (coverage) to 14.8Generating adjacencyAdded 12766 edges.Eroding tipsEroded 4295 tips.Eroded 0 tips.Hash load: 2066 / 8192 = 0.252 using 918 kB
Pruning tips shorter than 1 bp...
Pruning tips shorter than 2 bp...
Pruning tips shorter than 4 bp...
Pruning tips shorter than 8 bp...
Pruning tips shorter than 16 bp...
Removed 11 marked k-mer.Pruned 11 k-mer in 2 tips.Pruning tips shorter than 32 bp...Pruned 2 tips in 5 rounds.Hash load: 2055 / 8192 = 0.251 using 918 kBMarked 60 edges of 30 ambiguous vertices.Removing low-coverage contigs (mean k-mer coverage < 14.8)Found 2055 k-mer in 47 contigs before removing low-coverage contigs.Removed 416 k-mer in 13 low-coverage contigs.Split 26 ambigiuous branches.Hash load: 1639 / 8192 = 0.2 using 918 kBEroding tipsEroded 0 tips.Eroded 0 tips.Hash load: 1639 / 8192 = 0.2 using 918 kB
Pruning tips shorter than 1 bp...
Pruning tips shorter than 2 bp...
Pruning tips shorter than 4 bp...
Pruning tips shorter than 8 bp...
Pruning tips shorter than 16 bp...
Pruning tips shorter than 32 bp...Pruned 0 tips in 5 rounds.Hash load: 1639 / 8192 = 0.2 using 918 kBPopping bubblesRemoved 2 bubbles.Removed 2 bubblesMarked 0 edges of 0 ambiguous vertices.Assembled 1575 k-mer in 2 contigs.Removed 4722 k-mer.The signal-to-noise ratio (SNR) is -4.6 dB.AdjList -v -k32 -m50 --dot test-1.fa >test-1.dotReading `test-1.fa'...Finding overlaps of exactly k-1 bp...V=4 E=0 E/V=0Degree: ▇012340: 1e+02% 1: 0% 2-4: 0% 5+: 0% max: 0abyss-filtergraph -v --dot -k32 -g test-2.dot1 test-1.dot test-1.fa >test-1.pathLoading graph from file: test-1.dotGraph stats before:V=4 E=0 E/V=0Degree: ▇012340: 1e+02% 1: 0% 2-4: 0% 5+: 0% max: 0Removing shim contigs from the graph...Shim removal stats:Removed: 0 Too Complex: 0 Tails: 2 Too Long: 0 Self Adjacent: 0 Parallel Edges: 0Reading `test-1.fa'...Edge removal stats:Removed: 0Graph stats after:V=4 E=0 E/V=0Degree: ▇012340: 1e+02% 1: 0% 2-4: 0% 5+: 0% max: 0MergeContigs -v -k32 -g test-2.dot -o test-2.fa test-1.fa test-2.dot1 test-1.pathReading `test-2.dot1'...Read 4 vertices. Using 446 kB of memory.Reading `test-1.fa'...Read 2 sequences. Using 446 kB of memory.Reading `test-1.path'...Read 0 paths. Using 446 kB of memory.Writing `test-2.dot'...V=4 E=0 E/V=0Degree: ▇012340: 1e+02% 1: 0% 2-4: 0% 5+: 0% max: 0PopBubbles -v --dot -j2 -k32 -p0.9 -g test-3.dot test-2.fa test-2.dot >test-2.pathReading `test-2.dot'...V=4 E=0 E/V=0Degree: ▇012340: 1e+02% 1: 0% 2-4: 0% 5+: 0% max: 0Reading `test-2.fa'...Bubbles: 0 Popped: 0 Scaffolds: 0 Complex: 0 Too long: 0 Too many: 0 Dissimilar: 0V=4 E=0 E/V=0Degree: ▇012340: 1e+02% 1: 0% 2-4: 0% 5+: 0% max: 0MergeContigs -v -k32 -o test-3.fa test-2.fa test-2.dot test-2.pathReading `test-2.dot'...Read 4 vertices. Using 446 kB of memory.Reading `test-2.fa'...Read 2 sequences. Using 446 kB of memory.Reading `test-2.path'...Read 0 paths. Using 446 kB of memory.awk '!/^>/ {x[">" $1]=1; next} {getline s} $1 in x {print $0 "\n" s}' \test-2.path test-1.fa >test-indel.faln -sf test-3.fa test-unitigs.faabyss-map -v -j2 -l32 reads/rnaseq_1.fq.gz reads/rnaseq_2.fq.gz test-3.fa \|abyss-fixmate -v -l32 -h test-3.hist \|sort -snk3 -k4 \|DistanceEst -v -j2 -k32 -l32 -s200 -n10 -o test-3.dist test-3.histReading `test-3.fa'...Using 442 kB of memory and 2.21e+05 B/sequence.Reading `test-3.fa'...Building the suffix array...Building the Burrows-Wheeler transform...Building the character occurrence table...Read 1.67 kB in 2 contigs.Using 442 kB of memory and 265 B/bp.Reading from standard input...Mapped 11384 of 20000 reads (56.9%)Mapped 11374 of 20000 reads uniquely (56.9%)Read 20000 alignmentsMateless 0Unaligned 1965 19.6%Singleton 4686 46.9%FR 3349 33.5%RF 0FF 0Different 0Total 10000FR Stats mean: 495.6 median: 495 sd: 48.7 n: 3325 min: 353 max: 642 ignored: 24_▁_▁▁▁▂▂▂▃▄▃▄▄▆▇▅▆▇█▆█████▇████▇▇▆▅▅▆▅▄▃▄▂▃▂▂▂▁▁_ ▁ __Mate orientation FR: 3349 (100%) RF: 0 (0%)The library test-3.hist is oriented forward-reverse (FR).Stats mean: 495.6 median: 495 sd: 48.7 n: 3325 min: 353 max: 642_▁_▁▁▁▂▂▂▃▄▃▄▄▆▇▅▆▇█▆█████▇████▇▇▆▅▅▆▅▄▃▄▂▃▂▂▂▁▁_ ▁ __Minimum and maximum distance are set to -31 and 642 bp.DistanceEst: DistanceEst.cpp:619: int main(int, char**): Assertion `in' failed./bin/bash: line 3: 100465 Done abyss-map -v -j2 -l32 reads/rnaseq_1.fq.gz reads/rnaseq_2.fq.gz test-3.fa100466 | abyss-fixmate -v -l32 -h test-3.hist100467 | sort -snk3 -k4100468 Aborted | DistanceEst -v -j2 -k32 -l32 -s200 -n10 -o test-3.dist test-3.histmake: *** [test-3.dist] Error 134make: *** Deleting file `test-3.dist'The abyss-fix problem has disappeared, so I suppose that now the installation was ok and the error seems to be related to the DistanceEst command.Thank you again
Riccardo--
Dr. Riccardo Aiese Cigliano
Scientific Director and Co-funderSequentia Biotech SL
Parc Cientific de Barcelona
Torre R 02A6/A7
c/ Baldiri Reixac 4-8
08028 Barcelona
Cel. +34 633042534
raiesec...@sequentiabiotech.com
www.sequentiabiotech.com
This message is intended exclusively for its addressee and may contain information considered CONFIDENTIAL and protected by professional privilege. If you are not the intended recipient you are hereby notified that any dissemination, copy or disclosure of this communication is strictly prohibited by law. If this message has been mistakenly received, please immediately notify us via e-mail and delete it.