每日犯蠢例行记录

君子知夫不全不粹之不足以为美也

Posted by Chevy on January 8, 2021

background: 所里的老师另谋高就了,他负责的服务器也不给用了,那两天我帮大家备份了一波数据,没想到忘了把我的数据备份起来(事实上我实在的觉得我自己备份了)

今天整理的数据的时候,发现我之前测的一批RNA-seq数据好像没有备份在硬盘里,惊出一身冷汗。

平台老师走了,服务器也搬离了,原来这个世界上没有什么是可以永远存在的,当我意识到这件事的时候我的数据已经被我rm -rf ./*了。

好在之前导出到电脑里有两个bam文件:cry:(只占到四个数据里面的两个),所以亡羊补牢的把bam转化为fastq再备份一下吧,简单的做个记录:

# 首先使用samtools将bam文件按照name排序
samtools sort Ka-12-1.bam -n -O BAM -o Ka-12-1.sorted.bam -@ 4

# 随后再使用bedtools的bam2fastq将bam转化为fastq
bedtools bamtofastq -i Ka-12-1.sorted.bam -fq Ka-12-1_1.fastq -fq2 Ka-12-1_2.fastq

最后还得重新比对一下,反正有啥问题(到时候还得上传GEO,一切顺利吧希望)