Python通过调用windows命令行处理sam文件

时间:2021-11-27 13:34:10

Python通过调用windows命令行处理sam文件

以samtools软件为例

一、下载或者索取得到windows版本的samtools软件,解压后如下:

Python通过调用windows命令行处理sam文件

进入文件内部,有如下几个文件:

Python通过调用windows命令行处理sam文件

二、将samtools设置环境变量:

Python通过调用windows命令行处理sam文件

上图是设置环境变量的窗口,建议在 系统变量里更改,这样即使变更用户名,环境变量依然有效。

在cmd中输入samtools,出现samtools的命令说明,表明环境变量成功

Python通过调用windows命令行处理sam文件

三、 调试 通过cmd的命令行调用samtools命令:

A.将sam转化为bam文件  Samtools view -bS path<in>\**.sam -o path<out>\**.bam

处理完毕,查看path<out>,应该有sam文件如下:

Python通过调用windows命令行处理sam文件

B.将bam文件sort     samtools sort path<in>\**.bam path<out>\**.sort

处理完毕,查看path<out>,应该有.sort.bam文件如下:

Python通过调用windows命令行处理sam文件

C.将sort.bam文件转化成sort.sam文件samtools view -h path<in>\**bam > path<out>\**.sam

处理完毕,查看path<out>,应该有.sort.sam文件如下:

Python通过调用windows命令行处理sam文件

四、用python调用cmd命令行实现以上操作

编写代码如下:

Python通过调用windows命令行处理sam文件

运行结果应该产生三个文件:

Python通过调用windows命令行处理sam文件