测序质量linux命令?
使用deeptools将bam文件转换为bw文件
1、比对后得到的bam文件可以转化为bigWig(bw)格式,通过可视化软件进行展示。deeptools软件可以实现bw格式转化和可视化展示。
2、https://deeptools.readthedocs.io/en/develop/content/installation.html 此工具将读取BAM文件,并生成bigWig或bedGraph。覆盖率计算为每个bin的读取次数,bin为定义的连续读取的窗口大小。
3、该工具可用于评估给定样品的测序深度。它采样100万bp,计算重叠读数的数量,并可以报告一个直方图,告诉你有多少碱基被覆盖多少次。接受多个BAM文件,但它们都应该对应于相同的基因组组装。
samtools常用命令详解
samtools还有个非常重要的命令mpileup,以前为pileup。该命令用于生成bcf文件,再使用bcftools进行SNP和Indel的分析。bcftools是samtool中附带的软件,在samtools的安装文件夹中可以找到。
默认情况,samtools view输出序列到屏幕,可以重导向到别的文件。-b :声明输出为bam格式的文件。
建立索引 建立索引只需在Linux下输入命令:samtools faidx input.fa 这里序列文件为 input.fa,生成的索引文件以 .fai 结尾。
RSeQC使用
第一阶段是基础知识学习,找一本覆盖面广但是又不是很难啃的教材先对生物信息所涉及各个方面有所了解,比如人卫版李霞主编那本《生物信息学》。第二阶段是一个逐步深入的过程,这个过程中要学会工具的使用。
③ SAM 文件处理 使用 samtools 对 SAM 文件排序并转化为 BAM 文件。samtools是一个用于操作sam和bam文件的工具合集,包含有许多命令。④比对结果可视化 比对结果使用 IGV 、Genome Maps 和Sacant 等可视化查看。

相关推荐HOT
更多>>
命令查看linux日志?
linux查看日志的命令是?具体怎么用less 命令:可以打开大文件,支持上下翻页、搜索等功能,也可以用于查看日志文件。如下图所示,先cd到我们需...详情>>
2023-10-24 20:48:53
linux下编号的命令?
在Linux系统中使用fuser命令来列出进程编号在Linux中,可以使用pgrep命令来查询Firefox的进程号。pgrep命令可以根据进程名或者其他属性来查找进...详情>>
2023-10-24 19:43:03
linux追寻ip命令?
linux查看本机ip命令1、首先打开手机【设置】应用,然后找到【WLAN】,点击该项进入详细页面。接着在WLAN管理页面,则可以看到已经连接的WiFi网...详情>>
2023-10-24 18:07:23
linux命令大全中cd?
linux下cd是什么命令linux cd是linux光盘的意思 cd 命令设置某一进程的当前工作目录。用户必须具有指定目录中的执行(搜索)许可权。Linux cd ...详情>>
2023-10-24 16:56:05