生信技能43 - 基于BAM文件性染色体的样本性别判断

已于 2024-03-25 20:05:04 修改
阅读量218 收藏
点赞数
分类专栏： 生信分析项目实战技能合集 文章标签： 数据分析 python linux
于 2024-03-25 14:09:26 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/LittleComputerRobot/article/details/137003317
版权
本文介绍了如何利用BAM文件中的chrX和chrY reads比例来判断样本性别，通过分析高通量测序数据，当Y/X reads比例超过特定阈值时，判定为女性，反之为男性。详细步骤包括检测reads数量和预测性别。
摘要由CSDN通过智能技术生成

更多生信开发内容可查看本人专栏：

生信分析项目实战技能集合

1. 基于BAM文件性别判断原理

对于高通量测序数据，包括WGS和WES，根据bam文件里面的chrX和chrY的reas比例可以判断性别。如果chrX和chrY的reads比例超过一定阈值（例如100倍），则判定为女性，否则为男性。

2. 检测BAM文件chrX和chrY的reads数量

根据chrX / chrY染色体的reads比例，女性的Y/X染色体的reads比会显著大于男性。

# detect_gender.sh shell脚本
# 输入bam文件路径
bam_path=$1
# 输出包
了解本专栏 订阅专栏 解锁全文 超级会员免费看
生信与基因组学
关注
专栏目录
订阅专栏
生信技能51 - 基于BAM文件的微生物污染分析流程
LittleComputerRobot的博客
07-09 92
提取未比对到人类参考基因组的bam文件，将其比对到微生物参考基因组，根据比对情况查看是否存在污染
生信技能27 - BAM文件的批量合并
LittleComputerRobot的博客
09-27 809
当有多个样本的bam文件时，可以使用samtools的merge命令将这些bam文件进行合并为一个bam文件。
生信软件31 - Bcftools操作VCF/BCF文件高级用法合集
LittleComputerRobot的博客
08-24 1757
Bcftools操作VCF/BCF文件高级用法合集
生信技能26 - BWA比对、去除PCR重复和统计最佳实践
LittleComputerRobot的博客
08-14 365
BWA MEM 比对、去除PCR重复和统计最佳实践
linux怎么查看一个bam文件,生信分析过程中这些常见文件的格式以及查看方式你都知道吗？...
weixin_32052253的博客
04-28 6978
原标题：生信分析过程中这些常见文件的格式以及查看方式你都知道吗？生信分析过程中，会与很多不同格式的文件打交道，除了原始测序数据 fastq 之外，还需要准备基因组文件 fasta 格式和基因注释文件 gtf 格式。在分析的过程中还会有众多中间文件的生成，如 bed 、 bed12 、 sam 、 bam 、 wig 、 bigwig 、 bedgraph 等，生成后我们一般会查看下内容了解文件每一...
生信技能-高通量测序工具bam、samtools、bedtools及conda的下载和安装
weixin_43664814的博客
03-15 2839
用于建立 index；基于 BWT 算法，将 reads 比对到参考基因组；最新版本 bwa-mem2，Intel实验室对计算效率进行了优化。
生信搬运工-03-测序文件格式介绍
weixin_42149446的博客
09-13 4297
生信文件格式
生信软件14 - bcftools提取和注释VCF文件关键信息
LittleComputerRobot的博客
03-29 2552
VCFtools可用于变异信息的描述性统计，计算，过滤和格式转换。
生信分析过程中，这些常见文件的格式以及查看方式你都知道吗？
悟道西方
07-15 1082
生信分析过程中，会与很多不同格式的文件打交道，除了原始测序数据fastq之外，还需要准备基因组文件fasta格式和基因注释文件gtf格式。在分析的过程中还会有众多中间文件的生成，如bedbed12sambamwigbigwigbedgraph等，生成后我们一般会查看下内容了解文件每一列的含义，以此来决定需要提取哪些有用信息列来进行下一步分析。
医学数据分析实训 项目二 数据预处理预备知识(数据标准化处理,数据离差标准化处理,数据二值化处理,独热编码处理,数据PCA降维处理)
m0_73678713的博客
09-14 1386
使用StandardScaler进行数据预处理。
前端设计之 主页面、书架页面、数据分析页面
weixin_47505105的博客
09-14 1119
下面的代码我把导入第三方包相关的代码都省略了。
【Python 数据分析学习】Pandas的基础和应用（2）
2302_79527141的博客
09-14 1136
数据透视表是常用的数据分析工具之一，它可以根据一个或多个指定的维度对数据进行聚合。对于Pandas数据，由于其类似表格的特性，还支持许多累计统计的分析，为数据分析提供了很大的方便，也成Pandas的一大优势。的元素，而且不进行补齐操作，返回一个布尔型的对象。常见的比较运算有：‘>’，‘=’，‘
[Python数据可视化]Plotly Express: 地图数据可视化的魅力
最新发布
William数据分析的博客
09-17 889
在数据分析和可视化的世界中，地图数据可视化是一个强大而直观的工具，它可以帮助我们更好地理解和解释地理数据。Python 的 Plotly Express 库提供了一个简单而强大的方式来创建各种地图。本文将通过一个简单的示例，展示如何使用 Plotly Express 来创建一个交互式的地图，并探讨其在地图数据可视化方面的应用。
Tcl lnit error: Can’t find a usable init.tcl in the following directories 问题解决
梦想闹钟
09-15 544
实际研究后发现，其实py2exe已经把打包需要的lib放在dist文件夹下了，但是打包后的程序运行后却没有去lib下找，而是去找系统自带的环境变量里找，所以找不到。这个问题出现在我用py2exe打包了一个包含tkinter的图形化界面，在当前电脑上运行无问题，在移动到新电脑上后提示报错、getcwd用于获取当前工作目录绝对路径，在设置环境的变量的时候它用的是绝对路径-所以也导致了在当前电脑上能用而移动后不能用。解决方法是在你的程序里重新设置下环境变量，而且是用相对路径的形式。
Java企业面试题3
m0_74972727的博客
09-15 836
1. break和continue的作用(智*图) break：用于完全退出一个循环（如 for, while）或一个 switch 语句。当在循环体内遇到break语句时，程序会立即跳出当前循环体，继续执行循环之后的代码。 continue：用于跳过当前循环体中剩余的部分，并开始下一次循环。如果是在 for 循环中使用continue，则会直接进行条件判断以决定是否执行下一轮循环。 2. if分支语句和switch分支语句的异同之处(智*图) 相同点：都是用来根据不同的条件执行不同的代码块。
Python——俄罗斯方块
2302_81225694的博客
09-14 2400
这段代码使用了Pygame库来实现游戏的图形界面，通过键盘控制方块的移动和旋转。游戏循环不断更新方块的位置和网格状态，并绘制在屏幕上。在方块达到底部或无法继续移动时，判断是否有满行，并清除满行的方块。游戏会根据方块的状态和移动情况不断更新，直到无法继续下落为止，游戏结束。俄罗斯方块游戏是一款经典的益智游戏，通常使用编程语言Python来实现。请注意，这只是一个简单的示例，可能还有一些功能和优化方面的改进。您可以根据自己的需求进行修改和扩展。
文档内容识别系统源码分享
xuehaishijue的博客
09-17 1656
数据集信息展示在当今信息爆炸的时代，文档内容识别系统的有效性和准确性愈发重要。为了提升YOLOv8在文档内容识别任务中的表现，我们采用了名为“FULL 2”的数据集进行训练和评估。该数据集专为文档图像的多样性和复杂性而设计，涵盖了多种类别的文本元素，使其成为优化深度学习模型的理想选择。
2024.9.16 day 1 pytorch安装及环境配置
m0_58285219的博客
09-16 636
pytorch安装及环境配置
【生信技能树】GEO数据库挖掘 P5
07-03
【生信技能树】中的GEO（Gene Expression Omnibus）数据库挖掘是指从大规模基因表达谱数据集中获取和分析信息的过程。P5可能指的是第5级或某个具体的步骤。GEO是NCBI（美国国家生物技术信息中心）提供的一个公共...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

生信与基因组学

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值