amplab实验室的SNAP算法运行和操作
2016-01-08 17:07
1071 查看
amplab实验室的SNAP算法运行和操作:
环境:ubuntu14.04
1.安装:
参考【1】和【2】,需要:
g++ version 4.6
zlib 1.2.8 from http://zlib.net/
安装:
sudo apt-get install g++
sudo apt-get install zlib1g-dev
zlib没有直接的安装包,
2.编译:在snap目录下运行make
3.操作:
(1)构建参考基因的索引:
./snap-aligner index tests/datatest/datatest.fa index1运行结果:
hadoop@Mcnode1:~/cloud/adam/xubo/snap$ ./snap-aligner index tests/datatest/datatest.fa index1
Welcome to SNAP version 1.0beta.22.
Hash table slack 0.300000
Loading FASTA file 'tests/datatest/datatest.fa' into memory...0s
Saving genome...0s
Computing bias table.
Bias computation: 0 / 1202
Computed bias table in 0s
Allocating memory for hash tables...0s
Building hash tables.
0(0%) seeds occur more than once, total of 0(0%) genome locations are not unique, 498(41%) bad seeds, 0 both complements used 500 no string
Hash table build took 0s
Building overflow table.
Overflow table build and hash table save took 0s
Saving overflow table...0s
Index build and save took 0s (1202 bases/s)
(2)比对:
单端比对:
./snap-aligner single index1 tests/datatest/datatest.fq -o output.sam
运行结果:
hadoop@Mcnode1:~/cloud/adam/xubo/snap$ ./snap-aligner single index1 tests/datatest/datatest.fq -o output.sam
Welcome to SNAP version 1.0beta.22.
Loading index from directory... 0s. 1202 bases, seed size 20
Aligning.
Total Reads Aligned, MAPQ >= 10 Aligned, MAPQ < 10 Unaligned Too Short/Too Many Ns %Pairs Reads/s Time in Aligner (s)
2 2 (100.00%) 0 (0.00%)
双端比对:
./snap-aligner paired index1 tests/datatest/datatest.fq tests/datatest/datatest.fq -o output2.sam
运行结果:
hadoop@Mcnode1:~/cloud/adam/xubo/snap$ ./snap-aligner paired index1 tests/datatest/datatest.fq tests/datatest/datatest.fq -o output2.sam
Welcome to SNAP version 1.0beta.22.
Loading index from directory... 0s. 1202 bases, seed size 20
Aligning.
Total Reads Aligned, MAPQ >= 10 Aligned, MAPQ < 10 Unaligned Too Short/Too Many Ns %Pairs Reads/s Time in Aligner (s)
4 4 (100.00%) 0 (0.00%) 0 (0.00%) 0 (0.00%) 0.00%% 1,000 0
4.其他:snap文件夹下的docs有教程
参考文献:
【1】 https://github.com/amplab/snap
【2】 http://snap.cs.berkeley.edu/
环境:ubuntu14.04
1.安装:
参考【1】和【2】,需要:
g++ version 4.6
zlib 1.2.8 from http://zlib.net/
安装:
sudo apt-get install g++
sudo apt-get install zlib1g-dev
zlib没有直接的安装包,
2.编译:在snap目录下运行make
3.操作:
(1)构建参考基因的索引:
./snap-aligner index tests/datatest/datatest.fa index1运行结果:
hadoop@Mcnode1:~/cloud/adam/xubo/snap$ ./snap-aligner index tests/datatest/datatest.fa index1
Welcome to SNAP version 1.0beta.22.
Hash table slack 0.300000
Loading FASTA file 'tests/datatest/datatest.fa' into memory...0s
Saving genome...0s
Computing bias table.
Bias computation: 0 / 1202
Computed bias table in 0s
Allocating memory for hash tables...0s
Building hash tables.
0(0%) seeds occur more than once, total of 0(0%) genome locations are not unique, 498(41%) bad seeds, 0 both complements used 500 no string
Hash table build took 0s
Building overflow table.
Overflow table build and hash table save took 0s
Saving overflow table...0s
Index build and save took 0s (1202 bases/s)
(2)比对:
单端比对:
./snap-aligner single index1 tests/datatest/datatest.fq -o output.sam
运行结果:
hadoop@Mcnode1:~/cloud/adam/xubo/snap$ ./snap-aligner single index1 tests/datatest/datatest.fq -o output.sam
Welcome to SNAP version 1.0beta.22.
Loading index from directory... 0s. 1202 bases, seed size 20
Aligning.
Total Reads Aligned, MAPQ >= 10 Aligned, MAPQ < 10 Unaligned Too Short/Too Many Ns %Pairs Reads/s Time in Aligner (s)
2 2 (100.00%) 0 (0.00%)
双端比对:
./snap-aligner paired index1 tests/datatest/datatest.fq tests/datatest/datatest.fq -o output2.sam
运行结果:
hadoop@Mcnode1:~/cloud/adam/xubo/snap$ ./snap-aligner paired index1 tests/datatest/datatest.fq tests/datatest/datatest.fq -o output2.sam
Welcome to SNAP version 1.0beta.22.
Loading index from directory... 0s. 1202 bases, seed size 20
Aligning.
Total Reads Aligned, MAPQ >= 10 Aligned, MAPQ < 10 Unaligned Too Short/Too Many Ns %Pairs Reads/s Time in Aligner (s)
4 4 (100.00%) 0 (0.00%) 0 (0.00%) 0 (0.00%) 0.00%% 1,000 0
4.其他:snap文件夹下的docs有教程
参考文献:
【1】 https://github.com/amplab/snap
【2】 http://snap.cs.berkeley.edu/
相关文章推荐
- Basis小知识02
- sqlite中datetime 判断相等时的异常情况
- 精品软件 推荐 Ashampoo Snap 7.0.10 - 精美、实用的屏幕截图软件
- 遍历交换机里面的mac地址,与公司登记的合法mac地址进行对比
- LVM 快照对性能的影响
- Linux之LVM
- 希望博客新增加的功能——“网页预览”&SNAP
- 比对布线标准,推进行业发展
- 快速浏览Silverlight3 beta:鸡肋一样的WritableBitmap 推荐
- 对lvm快照的理解
- 利用LVM快照备份mysql数据库 推荐
- oVirt的VM在线快照实现浅析
- 基于mysqldump及lvm snapshot备份mysql数据库
- 11GR2新功能:DATAGURAD STANDBY的SNAPHOST测试 推荐
- 支持二次开发的Zigbee模块(SNAP技术)
- svg使用小计(二)
- 自动校验文件夹及其内容shell程序(md5 校验文件夹中所有数据)
- LVM 逻辑卷管理器
- SNAP格式
- SNAP建立索引的时候出现问题