yhlhhhhh -
每日与生物工程斗智斗勇到谢顶
综合讨论组
将微基因给出的全序CRAM文件转为bam
软件:samtools, cmd(windows), 终端(Mac), terminal(Linux)
温馨提示:
1. 此转换过程时间很长,至少我的数据就花了4.5h,所以如果是晚上转换,建议大家洗洗睡吧
2. 一定要准备一个比较大的U盘,因为转换出来的文件大概90G
3. 电脑容易发烫的同学一定要找个方法降温
4. 记得电脑接电源
安装:
Windows:
官网直接下载链接:http://www.htslib.org/download/
利用wget下载命令:
wget https://github.com/samtools/samtools/releases/download/1.11/samtools-1.11.tar.bz2
Mac:
利用brew下载命令:brew install samtools
Linux:
没试过,不知道23333
转换:
1. 下载参考序列:human_g1k_v37
2. 打开终端,输入命令:cd 你存储cram文件的文件夹路径
3. 之后输入命令:samtools view -T human_g1k_v37.fasta -b input.cram > output.bam
4. 等着吧
温馨提示:
1. 此转换过程时间很长,至少我的数据就花了4.5h,所以如果是晚上转换,建议大家洗洗睡吧
2. 一定要准备一个比较大的U盘,因为转换出来的文件大概90G
3. 电脑容易发烫的同学一定要找个方法降温
4. 记得电脑接电源
安装:
Windows:
官网直接下载链接:http://www.htslib.org/download/
利用wget下载命令:
wget https://github.com/samtools/samtools/releases/download/1.11/samtools-1.11.tar.bz2
Mac:
利用brew下载命令:brew install samtools
Linux:
没试过,不知道23333
转换:
1. 下载参考序列:human_g1k_v37
2. 打开终端,输入命令:cd 你存储cram文件的文件夹路径
3. 之后输入命令:samtools view -T human_g1k_v37.fasta -b input.cram > output.bam
4. 等着吧
7 个回复
想到之前有人来通过文件大小来质疑数据量,太不科学了
赞同来自: yhlhhhhh 、[已注销] 、大灰狼
赞同来自: yhlhhhhh 、蚀月
赞同来自: yhlhhhhh
若能找到且找对fasta/fa/fna(或者用bgzip压缩的对应gz格式)参考坐标文件则更好:
在Windows环境,建议首选Cygwin或安装Linux虚拟机,缺少前述条件的话再考虑在Powershell(加./)或cmd运行已经被编译好的samtools.exe
赞同来自: 大灰狼
赞同来自:
赞同来自:
赞同来自:
要回复问题请先登录或注册