1.解释

SSW是一个更快的SW算法，并且提供了c语言lib和java的调用

代码：

https://github.com/mengyao/Complete-Striped-Smith-Waterman-Library

论文：

Zhao, M., et al. (2013). "SSW Library: an SIMD Smith-Waterman C/C++ library for use in genomic applications." PloS one 8(12): e82138.

2.代码：

git clone https://github.com/mengyao/Complete-Striped-Smith-Waterman-Library.git
cd Complete-Striped-Smith-Waterman-Library
cd src
make java
java  -Djava.library.path=. -jar ssw.jar ssw.Example

3.结果：

xubo@xubo:~/xubo/tools/Complete-Striped-Smith-Waterman-Library/src$ java  -Djava.library.path=. -jar ssw.jar ssw.Example
Aligning nucleotides
score1=21 score2=8 ref_begin1=8 ref_end1=21 read_begin1=0 read_end1=14 ref_end2=4 cigar=9M1I5M

遇到的问题：

xubo@xubo:~/xubo/tools/Complete-Striped-Smith-Waterman-Library/src$ java -jar ssw.jar -Djava.library.path=.  ssw.Example
Cannot find libsswjni.so. Has the library been built and LD_LIBRARY_PATH or -Djava.library.path set appropriately?
java.lang.UnsatisfiedLinkError: no sswjni in java.library.path
Exception in thread "main" java.lang.UnsatisfiedLinkError: no sswjni in java.library.pathat java.lang.ClassLoader.loadLibrary(ClassLoader.java:1886)at java.lang.Runtime.loadLibrary0(Runtime.java:849)at java.lang.System.loadLibrary(System.java:1088)at ssw.Example.main(Example.java:12)

原因：

-Djava.library.path 应该放在-jar前面

参考

【1】https://github.com/xubo245
【4】http://blog.csdn.net/xubo245/

基因数据处理118之SSW运行相关推荐

基因数据处理123之SSW代码不正确，到时比SparkSW时间长
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释由于要生成新的score matrix:blosum50,第一次使用静态方法,直接传给align,到时每 ...
基因数据处理121之SSW的score matrix调整，使得与SparkSW评分一致
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释 SSW的评分矩阵是128*128的,是按char的int值来进行计算的.而blosum50是蛋白质的,而 ...
基因数据处理122之SSW和SparkSW评分不一致，query为Q9
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释 RT,但是顺序一致 2.代码: hadoop@Master:~/disk2/xubo/project/a ...
基因数据处理54之bwa-mem运行paird-end（1千万条100bp的reads）
指令: ```hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub/cs-bwamem$ bwa mem GRCH38BWAindex/GRCH38chr1L3 ...
基因数据处理56之bwa运行paird-end（1千万条100bp的reads）.md
(1)pair1.fq>sai bwa aln GRCH38BWAindex/GRCH38chr1L3556522.fasta g38L100c10000000Nhs20Paired1.fq & ...
基因数据处理106之bwa-mem运行paird-end（1千万条100bp的reads g38L100c10000000Nhs20Paired12）
脚本: hadoop@Master:~/xubo/project/alignment/sparkBWA$ cat g38L100c10000000Nhs20Paired12Bwamem.sh echo ...
基因数据处理59之snap运行single-end（1千万条100bp的reads）
记录: hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub/cs-bwamem$ snap-aligner single snap/snapindex g38 ...
基因数据处理120之scala调用SSW在linux下运行
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释先有java提供转换,使用jni调用c 然后scala调用java 2.代码: 2.1 java: pa ...
基因数据处理119之java调用SSW在linux下运行
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释测试自带Example: xubo@xubo:~/xubo/tools/Complete-Striped ...

基因数据处理118之SSW运行

1.解释

2.代码：

3.结果：

基因数据处理118之SSW运行相关推荐

最新文章

热门文章