代码如下:
说明:sys.argv[1]传入的第一个参数是包含引物名称,染色体及区间的文件。sys.argv[2]传入的第二个参数是由软件samtools分析获得的测序深度文件。运行该脚本的方法:python绝对路径 python脚本 参数文件1 参数文件2。该脚本的运行结果存储在meanandposdepth.txt中。
运行结果:
说明:部分结果截图。第一列:引物名称;第二列:染色体;第三列和第四列:引物在染色体上的起始和终止位置;第五列:引物区间测序深度;其余列两两一对,表示引物区间某位点的测序深度。