首页 > 其他分享 >bam文件去重复

bam文件去重复

时间:2023-09-22 23:13:16浏览次数:32  
标签:文件 utm 重复 sample01 gatk sorted bam

 

建库过程PCA扩增过程中引入重复序列,会对变异检测结果产生影响,重复的DNA片段会比对到参考基因组的相同位置,根据这一特点来进行去重复。

 

001、gatk(picard标记重复)

gatk MarkDuplicates -I sample01.sorted.bam -O sample01.sorted.markdup.bam -M sample01.sorted.markdup_metrics.txt

 

 

002、samtools

 

 

参考:

01、https://www.jianshu.com/p/8cdbb89530c6?utm_campaign=maleskine&utm_content=note&utm_medium=seo_notes

02、

 

 

 

 

 

标签:文件,utm,重复,sample01,gatk,sorted,bam
From: https://www.cnblogs.com/liujiaxin2018/p/17723575.html

相关文章

  • Hadoop是什么? Hadoop是一个由Apache开发的开源分布式计算框架,它能够处理大规模数据并
    Hadoop是什么?Hadoop是一个由Apache开发的开源分布式计算框架,它能够处理大规模数据并行处理任务,支持大规模数据存储和处理。Hadoop的核心组件包括分布式文件系统HDFS和分布式计算框架MapReduce,它们使得Hadoop可以在廉价的硬件上并行地处理大量数据。Hadoop还包括很多相关的项目和子......
  • k8s yaml文件总结
    k8s支持yaml和JSON格式创建资源对象,json用于接口之间消息传递,适用于开发;yaml格式用于配置和管理,适用于云平台管理,yaml简洁非标记性语言1.yaml相关基础概念  yaml语法规则:  大小写敏感;缩进表示层级关系;缩进不允许使用tab键,只允许使用空格;#表示注释---为可选分隔符,当需要......
  • 基于transbigdata包获取线路和站点的邻接表及shp文件
    importpandasaspdimportnumpyasnpimportgeopandasasgpdimporttransbigdataastbdimportjsonimportcsvimportnetworkxasnx#输入需要的城市place='深圳'#输入需要的线路line_range=[14,15,16,20]#遍历线路列表,挨个处理forxinline_range:try:......
  • 大文件断点续传、快传秒传实现方案
    前言为什么视频、音频、大型文档等大文件不能也直接上传吗,简单又方便?遇到手动暂停、网络中断、网络不稳定或者服务端响应超时,当你终于半天到99%,网络突然断开喜提从0%再来一次再者一次服务接受如此大的数据传输,不说服务器肯同意接收,即使配置同意接受这常常会使服务器出现响应超时......
  • C# 获取文件的类型(不是后缀)
    最近在做一个文件映射功能,发现使用FileInfo只能获取到文件的后缀,并没有具体的类型描述 可以以下方式获取usingSystem;usingSystem.Runtime.InteropServices;usingSystem.Windows;namespaceGetChineseExtension{///<summary>///MainWindow.xaml的交互逻......
  • Linux系统文件管理读书笔记
        在我初学Linux系统的过程中,我学到了一些基础知识,包括了解Linux系统的目录结构、掌握处理文件的基本命令,以及掌握vim编辑器的使用。这些知识对于理解和使用Linux系统非常重要,下面是我学到的一些关键点:Linux系统目录结构:Linux系统的目录结构是层级化的,它包括了许多不同......
  • WTM(基于Blazor)访问共享文件夹
    有个需求,需要访问内网网络共享文件夹中的文件,有域控限制。一开始直接在本地映射一个网络驱动器,然后像本地磁盘一样访问共享文件夹里的文件,比如:Y:\。然后直接在程序中访问共享文件夹中的文件,如下代码:DirectoryInfodirectoryInfo=newDirectoryInfo(@"\\10.10.65.14\文件夹"......
  • 简单的DVWA文件包含
    1.LOW难度,没有难度,直接就是文件包含了。主要是理清思路,为什么会感觉这里有文件包含,这里有一个page参数,当我们切换页面时变得是他的参数,而且是xxx.php那么大概率就可以认为是用一个函数把后面这个xxx.php文件包含读了进去。在另一个文件夹下构造一个phpinfo文件<?phpphpinfo?>......
  • linux命令创建文件
    Linux命令创建文件 在Linux系统中,有多种命令可以用来创建文件。下面将介绍几个常用的方法。1.使用touch命令创建文件:touch文件名该命令会创建一个空文件,如果文件已存在,则会更新文件的访问和修改时间。2.使用echo命令创建文件:echo"内容">文件名该命令会将指定的内......
  • 如何实现大文件断点续传、秒传
    大家先来了解一下几个概念:「文件分块」:将大文件拆分成小文件,将小文件上传\下载,最后再将小文件组装成大文件;「断点续传」:在文件分块的基础上,将每个小文件采用单独的线程进行上传\下载,如果碰到网络故障,可以从已经上传\下载的部分开始继续上传\下载未完成的部分,而没有必要从头开始......