该问题已被锁定!
2
关注
2958
浏览

Chip-seq bam文件的处理

为什么被折叠? 0 个回复被折叠
孟浩巍 超级管理员 用户来自于: 北京市
2023-05-30 09:56

现在的测序质量都很高,基本上不用做之前那种trim操作,就直接cut adapter序列之后比对即可。

chatGPT机器人 机器人 用户来自于: IANA
2023-05-29 15:09
对于Chip-seq bam文件的处理,其中 trim_galore 这部分的质控标准需要注意以下几点: 1. 保留高质量的reads:使用 -q 参数控制reads的最低质量阈值,一般设置为25或30,以保留高质量的reads。 2. 保留适当长度的reads:使用 --length 参数控制reads的最短长度,过短的reads可能会影响比对和后续分析,一般设置为25或30。 3. 去除低质量的reads:使用 -e 参数控制reads的错误率,一般设置为0.1或0.05,以去除低质量的reads。 4. 保留高保真性的reads:使用 --stringency 参数控制reads的保真性,一般设置为4或5,以保留高保真性的reads。 此外,根据fastqc结果可以做灵活调整,如根据reads长度分布情况调整 --length 参数,根据reads错误率分布情况调整 -e 参数等。 对于enhancer的鉴定,fastqc结果红色报错可能涉及到以下几个方面: 1. 过低的reads质量:如Per base sequence quality分布不均,Per sequence quality scores过低等,可能导致reads质量不高,影响后续分析结果。 2. 过低的reads数量:如Total Sequences数量过低,可能导致后续分析结果不够准确。 3. 过高的reads重复率:如Sequence Duplication Levels过高,可能导致后续分析结果不够准确。 因此,在进行enhancer鉴定时,需要对fastqc结果进行综合考虑,尽可能选择高质量、高覆盖度、低重复率的reads进行分析,以提高鉴定准确性。

关于作者

Zed999 注册会员

这家伙很懒,还没有设置简介

问题动态

发布时间
2023-05-29 15:02
更新时间
2023-05-31 22:18
关注人数
2 人关注

相关问题

公司双端测序的数据R1R2处理
从bed文件获取注释
两个samtools命令之间用管道连接,最后的输出文件被吞了。
人源细胞系ATAC-seq比对率较低,显示支原体污染,如何处理?
sambamba使用
如何按一个列表对基因型文件进行过滤,剔除不需要的样本?
是不是单细胞的chip-seq,cut run这些call peak是不需要对照的,我看代码里没有-c.所以Rna-seq的峰图也是不需要control来call peak的吗?
chip-seq数据下载有多个SRA
求助一个shell脚本问题,如何批量处理下面这种情况?
缺失值处理

推荐内容

ChIP-Seq deeptools热图信号比较为啥一般不做统计检验比较?
bowtie2 参考基因组注释 比对
bivalent domains 如何计算
chip-seq数据下载有多个SRA
WES数据下游分析和可视化结果展示
TWAS联合PWAS做 conditional analysis
chipseq中一共一般有多少峰,可以认为实验是比较靠谱的呢
染色体重叠区域问题
All Rights Reserved Powered BY WeCenter V4.1.0 © 2025