- 综合排序
- 最热优先
- 最新优先
- 来自专栏生信技能树
去rRNA可以解决GC含量双峰的右峰
前些天我在生信技能树提出来了一个转录组数据分析的疑难杂症：RNA-seq的fastq文件里面为什么有gc含量的双峰，就是fastq测序数据质量控制的时候发现了GC含量的双峰，然后我简单分析了那些高重复的 可以看到，GC含量最后一个峰是由rRNA导致，因为8百万条reads被去除后，该峰就消失了。 最后剩下的问题，就是GC含量的另外一个峰。我们后续再谈它！
1.9K40发布于 2020-06-11 - 来自专栏小汪Waud
测序数据GC含量异常该如何处理？
我们在对测序数据进行fastqc质控时，会比较关注样本的GC含量，较好的数据如下图所示 Normal data 事实上，我们的测序结果会受到建库方式、样本质量等很多因素的影响，其GC含量也会千差万别 下图是真实测序结果的fastqc报告，很明显其GC含量有多个峰值，相当的“红”。我们在进行后续分析之前，首先要对其进行一定处理，去掉杂峰的影响。 rawdata_qc 下图是双端测序结果的fastqc质控报告中的GC含量部分，我们可以看到有几个异常峰值。 read id do echo prefetch ${id} -O ./ done 参考资料： RNA-seq数据分析完全指北-03：去除奇怪的RNA RNA-seq的fastq文件里面为什么有gc 含量的双峰
4.2K20编辑于 2023-02-16 - 来自专栏生信技能树
生信基础50问-GC含量是否应该成正态分布呢
对NGS测序得到的成千上万条reads的GC含量进行统计，并且与正态分布形式对比已经是常规操作了，一般用fastqc软件即可分析，如下： ? 但现实情况往往不是这样。 Statistics Per base sequence quality Per sequence quality scores Per base sequence content Per base GC content Per sequence GC content Per base N content Sequence Length Distribution Sequence Duplication Levels Overrepresented sequences Kmer Content 那么问题来了，不同物种的ngs组学的GC含量理想范围是什么？
2K20发布于 2018-09-21 - 来自专栏生信技能树
RNA-seq的fastq文件里面为什么有gc含量的双峰
fastq测序数据质控的时候 首先fastq测序数据质量控制表格就发现质量差的可怜，而且居然有GC含量的双峰，如下： ? 遇到这样的情况，就必须单独看具体的每个样本，上面的GC含量图表是项目里面全部的样本的multiqc汇总图表。 我随机抽一个样本的fastqc报告看了看，如下： ? 这个59%到67%的GC含量reads有点多啊，就是第一个GC峰值。 然后80%附近GC含量reads也有点多，就是第二个GC含量峰值。 这两个GC含量峰值就是需要解决的问题，正常的RNA-seq数据不会出现这样的情况。 我们有必要仔细看看具体那些动辄一两百万的同一个GC含量的reads是什么，看了看，如下所示： Sequence Count Percentage Possible Source CCGGCCCGGAGCGGACGAGCCGCCCCGGCGGTGAACGGGGAGGAGGCGGG
2K50发布于 2020-06-09 - 来自专栏生信技能树
【直播】我的基因组47:测序深度和GC含量的关系
含量片段偏好性。 if $F[2]=~/[GC]/;$counts_sum{$key}+=$F[3];$number{$key}++;}END{print "$_\t$number{$_}\t$GC{$_}\t$counts_sum {$_}" foreach sort{$a<=>$b} keys %number}' 上面的代码写的不好，跑10万行需要 4s，跑一百万行需要36s，我估计把这8.9亿行的bam运行完，这样推算是10小时即可 我觉得我这次画的图还不错，很明显能看到这个趋势，GC含量比较高的窗口，有着相应比较高的测序深度！ ? 至此，完美的证明了文章开头的结论！ /tmp.txt') a$GC = a[,4]/a[,3] a$depth = a[,5]/a[,3] a = a[a$depth<100,] plot(a$GC,a$depth) library(ggplot2
1.8K110发布于 2018-03-08 - 来自专栏生信技能树
【直播】我的基因组47:测序深度和GC含量的关系
含量片段偏好性。 if $F[2]=~/[GC]/;$counts_sum{$key}+=$F[3];$number{$key}++;}END{print "$_\t$number{$_}\t$GC{$_}\t$counts_sum {$_}" foreach sort{$a<=>$b} keys %number}' 上面的代码写的不好，跑10万行需要 4s，跑一百万行需要36s，我估计把这8.9亿行的bam运行完，这样推算是10小时即可 我觉得我这次画的图还不错，很明显能看到这个趋势，GC含量比较高的窗口，有着相应比较高的测序深度！ ? 至此，完美的证明了文章开头的结论！ /tmp.txt')a$GC = a[,4]/a[,3]a$depth = a[,5]/a[,3]a = a[a$depth<100,]plot(a$GC,a$depth)library(ggplot2
3.3K90发布于 2018-03-08 - 来自专栏生信技能树
对参考基因组按照200k分区间统计测序深度及GC含量
以前是自己写脚本： 【直播】我的基因组47:测序深度和GC含量的关系 可能是太复杂，大多数读者表示看不懂，所以我重新使用已有的轮子来做这件事。 下载hg38参考基因组 直接谷歌搜索即可： ? 命令，代码如下 pip install pyfaidx faidx hg38.fasta -i chromsizes > sizes.genome 结果如下： chr1 248956422 chr10 含量 因为使用的是bedtools这样成熟的轮子， 所以就是一行代码而已： bedtools nuc -fi hg38.fa -bed 200k.bed | cut -f 1-3,5 > 200k_gc.bed # 4_pct_at 5_pct_gc 6_num_A 7_num_C 8_num_G 9_num_T 10_num_N 文件如下： $head 200k_gc.bed #1_usercol 2_usercol 3_usercol 5_pct_gc chr1 0 200000 0.420110 chr1 200000 400000 0.220065 chr1
4.4K20发布于 2019-06-19 - 来自专栏InCerry
为 .NET 10 GC（DATAS）做准备
字数 8728，阅读大约需 44 分钟 原作者:maoni 原文链接：https://maoni0.medium.com/preparing-for-the-net-10-gc-88718b261ef2 启用 DATAS 后，他们观察到吞吐量下降了 6.8%，同时工作集减少了 10%。 假设我们希望在触发 GC 前允许分配 20MB，这意味着乘数是 10。但如果 LDS 是 20GB，我们也不希望分配 200GB 才 GC，这意味着乘数要小得多。 以下是不同 LDS 值的一些示例 — LDS (MB) m m 限制后 BCD (MB) 1 15.00 10.00 10 5 6.71 6.71 34 10 4.74 4.74 47 50 2.12 如果指定了 GCHeapCount[10]，DATAS 会被禁用，因为它告诉 GC 不要修改堆数量。而修改堆数量是 DATAS 调整性能的关键机制之一，所以这是关闭 DATAS 的信号。
33900编辑于 2025-11-12 - 来自专栏简说基因
生物信息学算法之Python实现|Rosalind刷题笔记：005 GC含量计算
DNA 序列的 GC 含量是指序列中'G'和'C'所占的百分比。 一条 DNA 序列很容易表示，但是如果有多条 DNA 序列放在一起，则每条序列必须被标记，通常的做法是保存为 FASTA 格式文件。 需得：GC 含量最高的序列名称及其 GC 含量（各占一行行输出）。 ("cout_gc_content:Failed") sys.exit(1) item = max_gc_content('rosalind_gc.txt') print ，特别是为其构造一个 key 函数并传入，这是解本题的关键，GC 含量本身是很容易理解的。 Given: At most 10 DNA strings in FASTA format (of length at most 1 kbp each).
1.5K20发布于 2020-12-14 测试没技术含量？
哈哈，这个话题带点江湖气息：开发有时觉得“自己写代码是造房子”，而测试“就是来找茬的”，于是容易说出“测试没技术含量”这种欠扁的话。 面对“测试没技术含量”这类轻视言论，不需要情绪化“怼”，而是用理性、专业、有深度的回应让对方闭嘴甚至肃然起敬。 他们没技术含量？那你去应聘试试？”“微软的测试工程师要会写编译器插件做静态分析，Amazon的QA要会建分布式压测平台 —— 你觉得这是‘点点点’？” ✍️ 总结：你可以这样优雅回击“测试当然有技术含量 —— 它是在不确定性中寻找确定性，在混沌中建立秩序，在沉默的系统中听见即将爆炸的声音。undefined我们不生产代码，我们生产信心。 谁敢说显微镜没技术含量，那是眼界问题。
15410编辑于 2025-09-18- 来自专栏生物信息学
R语言计算GCAT含量
这个包主要是处理基因组的一些序列信息，包括：序列翻译、DNA/RNA互转、统计各个碱基的含量、三连字母的含量.....这些都是一行命令可以解决的。今天就先来教大家怎样计算GC/AT含量。 ? 安装完毕，只需敲几行代码，就可以实现GC/AT含量可视化。 /dell/Desktop/sequence.fasta"#读取文件（FASTA格式）x<-readDNAStringSet(filepath)chrom<-x[[1]]#每100个碱基为窗口计算AT含量
2K10发布于 2020-04-13 - 来自专栏小明的数据分析笔记本
使用Django框架开发的第一个简易web程序：计算fasta序列长度和GC含量
Django是python写成的开放源代码的web应用框架 本次使用的基础配置 python版本：3.6.3 Django:2.2 Windows10系统 如何查看Django版本 import django
{% for line in seqLength 'fasta'): d = {"SeqName":rec.id,"SeqLength":str(len(rec.seq)),\ "GCcontent":str(GC序列名称 序列长度 GC 含量
83010发布于 2020-03-02来自专栏小明的数据分析笔记本杂记：ggpairs更改配色；ggplot2极坐标添加直线；seqkit计算fasta序列的长度和gc含量
计算gc含量 seqkit.exe fx2tab --name --only-id --gc output.fasta -o gc.txt ?
2.1K30发布于 2021-07-12来自专栏老男孩成长之路Minor GC、Major GC、Full GC的区别
今天主要谈谈JVM GC的类型和策略，特别是大家经常混淆的Minor GC、Major GC、Full GC，年轻代GC、老年代GC，之间有什么区别和联系。 Major GC 老年代的垃圾收集叫做Major GC，Major GC通常是跟full GC是等价的，收集整个GC堆。 Minor GC和Major GC其实就是年轻代GC和年老年GC的俗称。 而在Hotspot VM具体实现的收集器：Serial GC, Parallel GC, CMS, G1 GC中，大致可以对应到某个Young GC和Old GC算法组合。 分代GC并不收集整个GC堆的模式，而是只专注分代收集 Young GC：只收集年轻代的GC Old GC：只收集年老代的GC(只有CMS的concurrent collection是这个模式) Mixed 以上就是Minor GC、Major GC、Full GC的一个介绍，后续我们再重点介绍JVM GC相关的垃圾回收算法以及底层的实现。
8.1K73发布于 2021-07-21来自专栏生信开发者胚胎DNA含量低的原因
Low fetal fraction has been linked to very early gestations, high maternal BMI, maternal medications, smoking and factors which lead to a smaller placenta, such as trisomies 13 and 18 (Kuhlmann-Capek et al. 2019)
69310发布于 2021-01-05来自专栏程序猿DDMinor GC、Major GC和Full GC之间的区别
如果正好相反，Eden 区大部分新生对象不符合 GC 条件，Minor GC 执行时暂停的时间将会长很多。 所以 Minor GC 的情况就相当清楚了——每次 Minor GC 会清理年轻代的内存。 首先，许多 Major GC 是由 Minor GC 触发的，所以很多情况下将这两种 GC 分离是不太可能的。 这使得我们不用去关心到底是叫 Major GC 还是 Full GC，大家应该关注当前的 GC 是否停止了所有应用程序的线程，还是能够并发的处理而不用停掉应用程序的线程。 34048.0 34048.0 34048.0 0.0 272640.0 129480.8 1756416.0 745100.2 20608.0 19704.5 2944.0 2678.4 10
1.1K30发布于 2018-12-24来自专栏爱撸猫的杰Minor GC、Major GC和Full GC之间的区别
Minor GC 从年轻代空间（包括 Eden 和 Survivor 区域）回收内存被称为 Minor GC。这一定义既清晰又易于理解。 如果正好相反，Eden 区大部分新生对象不符合 GC 条件，Minor GC 执行时暂停的时间将会长很多。 所以 Minor GC 的情况就相当清楚了——每次 Minor GC 会清理年轻代的内存。 首先，许多 Major GC 是由 Minor GC 触发的，所以很多情况下将这两种 GC 分离是不太可能的。 这使得我们不用去关心到底是叫 Major GC 还是 Full GC，大家应该关注当前的 GC 是否停止了所有应用程序的线程，还是能够并发的处理而不用停掉应用程序的线程。 34048.0 34048.0 34048.0 0.0 272640.0 129480.8 1756416.0 745100.2 20608.0 19704.5 2944.0 2678.4 10
2K31发布于 2019-03-28来自专栏福大大架构师每日一题2020-10-04：java中GC Root 对象有哪些？
福哥答案2020-10-04：#福大大架构师每日一题# 简单回答： 1.局部变量。 2.静态引用的对象。 3.常量引用的对象。 4.JNI引用的对象。
1.8K31发布于 2020-10-04来自专栏余生大大Minor GC、Major GC和Full GC有何不同
三种GC GC类型 GC区域 触发条件 Stop The World时间 Minor GC Eden 和 Survivor 区域 Eden区域 > 设定内存阈值 对于大部分应用程序，Minor GC停顿导致的延迟都是可以忽略不计的 Major GC Old区域 根据不同的垃圾收集器配置由Minor GC触发 MajorGC 的速度一般会比 Minor GC 慢 10倍以上。 Full GC 整个Heap空间包括年轻代和永久代 调用System.gc时Old老年代空间不足方法区空间不足通过Minor GC后进入老年代的平均大小大于老年代的可用内存 Full GC作用于整个堆空间的 首先，许多 Major GC 是由 Minor GC 触发的，所以很多情况下将这两种 GC 分离是不太可能的。 34048.0 34048.0 34048.0 0.0 272640.0 129480.8 1756416.0 745100.2 20608.0 19704.5 2944.0 2678.4 10
80970编辑于 2022-11-02来自专栏芋道源码1024Minor GC、Major GC 和 Full GC 之间的区别
Minor GC 从年轻代空间（包括 Eden 和 Survivor 区域）回收内存被称为 Minor GC。这一定义既清晰又易于理解。 如果正好相反，Eden 区大部分新生对象不符合 GC 条件，Minor GC 执行时暂停的时间将会长很多。 所以 Minor GC 的情况就相当清楚了——每次 Minor GC 会清理年轻代的内存。 首先，许多 Major GC 是由 Minor GC 触发的，所以很多情况下将这两种 GC 分离是不太可能的。 这使得我们不用去关心到底是叫 Major GC 还是 Full GC，大家应该关注当前的 GC 是否停止了所有应用程序的线程，还是能够并发的处理而不用停掉应用程序的线程。 34048.0 34048.0 34048.0 0.0 272640.0 129480.8 1756416.0 745100.2 20608.0 19704.5 2944.0 2678.4 10
70910发布于 2019-10-17腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059
粤公网安备44030502008569号腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号