
將微基因給出的全序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 個回復
?
贊同來自:
贊同來自:
? ?
若能找到且找對fasta/fa/fna(或者用bgzip壓縮的對應gz格式)參考坐標文件則更好:
? 在Windows環境,建議首選Cygwin或安裝Linux虛擬機,缺少前述條件的話再考慮在Powershell(加./)或cmd運行已經被編譯好的samtools.exe
贊同來自: 大灰狼
贊同來自:
贊同來自: yhlhhhhh 、蝕月
贊同來自: yhlhhhhh
想到之前有人來通過文件大小來質疑數據量,太不科學了
贊同來自: yhlhhhhh 、[已注銷] 、大灰狼 、微信用戶_740549
要回復問題請先登錄或注冊