samtools得到mapping中各个位置覆盖度情况程序，samtools tview的使用

最新推荐文章于 2024-07-29 11:14:19 发布

morningCV

最新推荐文章于 2024-07-29 11:14:19 发布

阅读量7k

点赞数 2

分类专栏：生物编程jin

本文链接：https://blog.csdn.net/g_r_c/article/details/8740416

版权

编程jin 同时被 2 个专栏收录

52 篇文章 0 订阅

订阅专栏

生物

27 篇文章 4 订阅

订阅专栏

use strict;
use warnings;


system 'samtools tview /home/blackstar/lastz/GRC/Yeast/share/02.assembly/00.novoalign/yeast_set1/H1_1/H1_1_Clean.S288C.novo.pe.dedup.bam /home/blackstar/lastz/GRC/Yeast/share/ref/S288C_R64.fasta';

这里面，前面的参数我们基本不用看，你可以自己选择，

-f 后面跟ref序列

-l 后面跟你的问询序列，序列的格式就是前面文章中的bed文件格式

然后后面要跟bam文件

最后是输出文件。

use strict;
use warnings;

system 'samtools mpileup -ABQ0 -d100000 -f /home/blackstar/lastz/GRC/Yeast/share/ref/S288C_R64.fasta -l a.bed /home/blackstar/lastz/GRC/Yeast/share/02.assembly/00.novoalign/yeast_set1/H1_1/H1_1_Clean.S288C.novo.pe.dedup.bam >PREFIX.gatk.flt.spec.snp.around10.mpileup';

use strict;
use warnings;

system 'samtools tview -s H1_1_Clean /opt/share/yeast/02.assembly/02.processed/01.stampy/yeast_set1/S288C/yeast_set1.S288C.stampy.pe.realn.bam /opt/share/yeast/ref/S288C_R64.fasta';