unique mapped reads

就是指唯一比对的reads

现在人们已经开始避免使用unique mapped reads这个概念了，而转向使用mapq值来保留高质量的比对结果。因为mapq值反应了一组比对结果发生的可能性，MapQ = -10 log10(P), 比如结果为10，那就是1/10的概率会出现这个比对结果，如果我们认为0.05%是一个小概率的话，那个mapq值为15就可以用于筛选了, 如果认为0.01%是个小概率的话，mapq值为20就可以用于筛选了。但是人们往往从30这个值开始试起（1/1000的概率），如果它的筛选结果符合你的测序要求，就可以使用它了。如果不行，可以适当的调整这个筛选值。

samtools view -bhS -q 30 input.sam > output.bam

参考：关于map当中的unique mapped reads问题

unique mapped reads相关推荐

文献阅读---玉米干旱响应和耐受性基因表达的调控变异定位
文献阅读-玉米干旱响应和耐受性基因表达的调控变异定位 Mapping regulatory variants controlling gene expression in drought respon ...
samtools 的使用
1,对fasta文件建立index samtools faidx ref.fasta 注意, 如果你执行命令的地方和参考序列不在同一个目录,参考序列用全路径过相对路径,最后的index结果和参考序列 ...
RNA-seq流程学习笔记（10）-使用HTSeq-count软件对reads进行计数
参考文章: 转录组入门(6): reads计数 RNA-seq练习第二部分 RNA-seq(6): reads计数转录组学习六(reads计数与标准化) 1. reads计数的原理对我们测序得到 ...
RNA-seq第四期——HTSeq-count对reads进行计数
对测序得到的reads进行计数,即基因表达的定量过程.根据reads和基因位置的overlap,以此来判断reads到底属于哪一个基因,同时对该reads总数进行计数,生成counts矩阵. 今日内容 ...
生物信息学（高通量测序）名词
什么是高通量测序高通量测序技术(High-throughput sequencing,HTS)是对传统Sanger测序(称为一代测序技术)革命性的改变, 一次对几十万到几百万条核酸分子进行序列测定, ...
signature=1c1223d10236c18f125cb9f7a490588e,Microscaled proteogenomic methods for precision oncology
Patient-derived xenografts and drug treatment For PDX studies, all animal procedures were approved b ...
RNA 1. SCI 文章中读取 GEO 数据
依稀记得10年前一个样本的RNA-SEQ费用还是蛮高的,而现在就是洒洒水啦,所以这种转录组的数据也已经成为文章的主流,基本上就是WES结合RNA,但是个人感觉,真的能关联到突变基因和表达的水平上去的, ...
Coelho2021_GMGCv1 计算丰度的方法
文章链接复旦大学赵兴明领衔团队Nature发文构建全球微生物基因目录 Coelho, L.P., Alves, R., del Río, Á.R. et al. Towards the biogeo ...
【生信】R语言在RNA-seq中的应用
R语言在RNA-seq中的应用文章目录 R语言在RNA-seq中的应用生成工作流环境读取和处理数据由targets文件提供实验定义对实验数据进行质量过滤和修剪生成FASTQ质量报告比对 ...
RNA-seq中的基因表达量计算和表达差异分析
差异分析的步骤: 1)比对: 2) read count计算: 3) read count的归一化: 4)差异表达分析: 背景知识: 1)比对: 普通比对: BWA,SOAP 开大GAP比对:Toph ...

unique mapped reads

unique mapped reads相关推荐

最新文章

热门文章