1.解释

RT,但是顺序一致

2.代码：

hadoop@Master:~/disk2/xubo/project/alignment/SparkSW/SparkSW20161114/alluxio-1.3.0$ vi sparkSWSSWQ10.sh dir=$(ls -l query/* |awk '!/^d/ {print $NF}' |grep -i '/[9-9]' )
for i in 1 2 3 4 5 6 7 8 9 10
do
echo 'SparkSW:'
#for j in 1 2 3 4 5 6 7 8 9 10
for j in 1 2 3
doecho 'time:'$jy='/xubo/project/SparkSW/input/D'$i'Line.fasta'for file in $dirdosh test.sh 8G /xubo/project/SparkSW/hash50.csv $file $y 128 1 5done
doneecho 'SSW'
#for j in 1 2 3 4 5 6 7 8 9 10
for j in 1 2 3
doecho 'time:'$jy2='/xubo/project/SparkSW/input/D'$i'Line.fasta'#y2='alluxio://Master:19998/xubo/project/SparkSW/input/D'$i'Line.fasta'for file in $dirdosh testSSW.sh 8G /xubo/project/SparkSW/hash50.csv $file $y2 128 1 5done
done
done

3.结果：

hadoop@Master:~/disk2/xubo/project/alignment/SparkSW/SparkSW20161114/alluxio-1.3.0$ tail -f sparkSWSSWQ10time201611271348.txt
SparkSW:
time:1
(9363,UniRef100_A1CLY8)
time:2
(9363,UniRef100_A1CLY8)
time:3
(268435464,UniRef100_Q5AWI5)
(9363,UniRef100_A1CLY8)
SSW
time:1
(9588,UniRef100_A1CLY8)
time:2
(9588,UniRef100_A1CLY8)
time:3
(9588,UniRef100_A1CLY8)
16/11/27 14:45:28 WARN QueuedThreadPool: 3 threads could not be stopped
SparkSW:
time:1
(9363,UniRef100_A1CLY8)
(7756,UniRef100_S0EEY3)
(7695,UniRef100_W7MLD7)
time:2
(9363,UniRef100_A1CLY8)
(7756,UniRef100_S0EEY3)
(7695,UniRef100_W7MLD7)
time:3
(9363,UniRef100_A1CLY8)
(7756,UniRef100_S0EEY3)
(7695,UniRef100_W7MLD7)
SSW
time:1
(9588,UniRef100_A1CLY8)
(8014,UniRef100_S0EEY3)
(7958,UniRef100_W7MLD7)
time:2
(9588,UniRef100_A1CLY8)
(8014,UniRef100_S0EEY3)
(7958,UniRef100_W7MLD7)
time:3
(9588,UniRef100_A1CLY8)
(8014,UniRef100_S0EEY3)
(7958,UniRef100_W7MLD7)
SparkSW:
time:1
(268435461,UniRef100_B1YQQ1)
(9363,UniRef100_A1CLY8)
(7756,UniRef100_S0EEY3)
(7695,UniRef100_W7MLD7)
(6660,UniRef100_Q4WAZ9)
time:2
(268435455,UniRef100_Q48ME3)
(9363,UniRef100_A1CLY8)
(7756,UniRef100_S0EEY3)
(7695,UniRef100_W7MLD7)
(6660,UniRef100_Q4WAZ9)
time:3
(9363,UniRef100_A1CLY8)
(7756,UniRef100_S0EEY3)
(7695,UniRef100_W7MLD7)
(6660,UniRef100_Q4WAZ9)
(5404,UniRef100_Q0C8M3)
SSW
time:1
(9588,UniRef100_A1CLY8)
(8014,UniRef100_S0EEY3)
(7958,UniRef100_W7MLD7)
(6905,UniRef100_Q4WAZ9)
(5620,UniRef100_Q0C8M3)
time:2
(9588,UniRef100_A1CLY8)
(8014,UniRef100_S0EEY3)
(7958,UniRef100_W7MLD7)
(6905,UniRef100_Q4WAZ9)
(5620,UniRef100_Q0C8M3)
time:3
(9588,UniRef100_A1CLY8)
(8014,UniRef100_S0EEY3)
(7958,UniRef100_W7MLD7)
(6905,UniRef100_Q4WAZ9)
(5620,UniRef100_Q0C8M3)

参考

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

基因数据处理122之SSW和SparkSW评分不一致，query为Q9相关推荐

基因数据处理121之SSW的score matrix调整，使得与SparkSW评分一致
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释 SSW的评分矩阵是128*128的,是按char的int值来进行计算的.而blosum50是蛋白质的,而 ...
基因数据处理123之SSW代码不正确，到时比SparkSW时间长
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释由于要生成新的score matrix:blosum50,第一次使用静态方法,直接传给align,到时每 ...
基因数据处理118之SSW运行
更多代码请见:https://github.com/xubo245 基因数据处理系列 1.解释 SSW是一个更快的SW算法,并且提供了c语言lib和java的调用代码: https://github ...
基因数据处理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 ...
基因数据处理12之samtool的tview来查看sam的匹配文件
基因数据处理12之samtool的tview来查看sam的匹配文件具体的之前有文章讲过:http://blog.csdn.net/xubo245/article/details/50836185 记 ...
基因数据处理8之BWA_MEM小数据集处理（成功）
基因数据处理8之BWA_MEM小数据集处理环境:ubuntu14.04 6G内存参考基因:GRCH38 来源请参考[1] 1.fastq数据:SRR003161.fastq 的头20行,即5条re ...
基因数据处理1之mapping_to_cram
基因数据处理1之mapping_to_cram 参考资料: A Worked Example Obtain some public data We will use the first 100,000 ...
基因数据处理44之cloud-scale-bwamem安装
cloud-scale-bwamem是在spark等云环境上实现bwa-mem算法,加速对全基因组匹配的处理 1.下载: git clone https://github.com/ytchen0323 ...

基因数据处理122之SSW和SparkSW评分不一致，query为Q9

1.解释

2.代码：

3.结果：

基因数据处理122之SSW和SparkSW评分不一致，query为Q9相关推荐

最新文章

热门文章