reads数越多，滑窗内GC含量降低 - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

634: 浏览

reads数越多，滑窗内GC含量降低

[size=15]代码：samtools mpileup -f ~/reference/genome/hg19/hg19.fa ../P_jmzeng.final.bam|[b]head -1000000[/b] |perl -alne '{$pos=int($F[1]/20000); $key="$F[0]\t$pos";$GC{$key}++ 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}'[/size] [size=15]当取前100万行时GC含量正常，但取的行数越多GC含量越低，请问是代码问题还是数据本身原因[/size]

好问题 0 评论收藏举报

1 回答

孟浩巍超级管理员用户来自于: 北京市
2018-08-29 09:58

perl代码…… pass下一题

赞同 0 0评论

关于作者

: 皆坐注册会员
这家伙很懒，还没有设置简介

1: 回答

0: 文章

2: 问题

问题动态

发布时间: 2018-08-29 09:02

更新时间: 2018-08-29 09:58

关注人数: 2 人关注

相关问题

如何提取比对上某一个基因或者某一段序列的所有reads: 558 浏览 2 关注 1 回答 0 评论

sgRNA的GC content 是什么意思，该如何计算: 943 浏览 2 关注 2 回答 0 评论

RNA-seq reads和表达量在转座子上的分布meta图: 593 浏览 2 关注 1 回答 0 评论

问题：如何筛选WGCNA分析得出的consensus module中的hub基因？: 872 浏览 2 关注 1 回答 1 评论

RNA Seq中什么样的reads是chimeric reads？是不是多数的chimeric reads都是circRNA的？: 685 浏览 2 关注 1 回答 0 评论

featureCounts assigned reads率低，正常吗？: 1184 浏览 1 关注 1 回答 0 评论

[求助]WGCNA: 527 浏览 2 关注 1 回答 0 评论

sudo apt-get install gcc++出现如下报错该如何解决？系统是：ubuntu14 server 版: 717 浏览 2 关注 1 回答 0 评论

WGCNA: 771 浏览 2 关注 2 回答 0 评论

WGCNA筛选到感兴趣的模块后，如何找这个模块的hub基因呢？: 629 浏览 2 关注 1 回答 0 评论

推荐内容

问 T to T测序中ONT测序起到的作用: 1478 浏览 2 关注 2 回答 0 评论

问可否用sanger测序测某一基因的上下游序列: 703 浏览 2 关注 1 回答 0 评论

问从bam文件还能提取出来原来比对前的fastaq文件吗？: 1207 浏览 2 关注 2 回答 0 评论

问第二代测序碱基出错偏向的原因是什么为什么第二代测序碱基出错主要以替换为主稀有碱基会影响测序吗: 865 浏览 2 关注 2 回答 0 评论

问高度杂合变异材料的遗传变异研究采取方法讨论: 707 浏览 2 关注 1 回答 0 评论

问 samtools view筛选cellranger比对结果: 677 浏览 2 关注 1 回答 0 评论

问 bowtie2比对结果报告best，以及aligned concordantly: 718 浏览 2 关注 2 回答 0 评论

问常用的对测序序列进行聚类分析的软件有什么？及其使用的了什么算法？: 1015 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024