fastquniq去除pcr重复效果

最新推荐文章于 2024-03-26 07:03:15 发布
VIP文章 最新推荐文章于 2024-03-26 07:03:15 发布
阅读量483 收藏
点赞数
文章标签： python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43504665/article/details/128447142
版权

这次组装的数据clean data差不多有40G的数据量，可惜使用megahit组装却只得到几十兆大小的文件，显然这里面有很多重复序列，刚好测试下fastuniq

默认参数

fastuniq -i 1.txt -t q -o output_1.fastq -p output_2.fastq -c 1

没什么改变

file            format  type     num_seqs         sum_len  min_len  avg_len  max_len
clean_R1.fastq  FASTQ   DNA   129,257,421  19,209,488,541       50
最低0.47元/天 解锁文章
ZZZpiinocchio
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
FastUniq:fastuniq-开源
05-30
用于配对短 DNA 序列的超快从头重复去除工具
MATLAB模拟PCR过程
07-25
MATLAB模拟PCR过程是指使用MATLAB软件来建立数学模型和进行数值计算，以模拟聚合酶链反应（PCR）的扩增过程。PCR是一种在生物学和分子生物学中常用的技术，用于扩增特定的DNA片段。 在MATLAB中，可以通过定义PCR模型方程和设置PCR反应参数，来模拟PCR过程中DNA的扩增情况。PCR模型方程通常基于指数函数或其他适当的数学表达式，用来描述DNA浓度随着PCR循环次数的变化。 通过设置初始DNA浓度、PCR循环次数和DNA延伸时间等参数，MATLAB可以计算每个PCR循环后DNA浓度的变化，并记录下所有的扩增结果。这样，研究人员可以通过模拟PCR过程来了解不同条件下PCR反应的动力学，例如扩增效率、最终产物浓度等。 MATLAB还提供了绘图工具，可以将PCR模拟结果可视化为曲线，便于观察和分析扩增效果。通过绘制PCR结果曲线，研究人员可以更直观地了解PCR过程中DNA的扩增情况，并根据需要对PCR实验条件进行优化和调整。 总之，MATLAB模拟PCR过程提供了一种灵活、可定制的方法，用于研究和优化PCR反应。它可以帮助研究人员更好地理解PCR过程中DNA扩增的动
基于reads的ID对fastq文件进行去重。
weixin_41869644的博客
01-22 2007
#!/usr/bin/perl use warnings; use Getopt::Long; my $usage = qq{ This script takes an input FASTQ file (formatted with four lines per sequence) and removes duplicate entries based on ID. Usage: perl...
MPB：微生物所蔡磊组-​​基于二代测序的真菌基因组组装和注释
刘永鑫的博客——宏基因组公众号
01-08 7354
为进一步提高《微生物组实验手册》稿件质量，本项目新增大众评审环节。文章在通过同行评审后，采用公众号推送方式分享全文，任何人均可在线提交修改意见。公众号格式显示略有问题，建议电脑端点击文末...
纯二代测序从头组装基因组
weixin_33737774的博客
03-09 3480
基因组组装 基因组组装一般分为三个层次，contig, scaffold和chromosomes. contig表示从大规模测序得到的短读(reads)中找到的一致性序列。组装的第一步就是从短片段(pair-end)文库中组装出contig。进一步基于不同长度的大片段(mate-pair)文库，将原本孤立的c...
使用二代数据进行基因survey-01cleandata（持续更新中）
weixin_44616693的博客
11-24 1421
使用二代数据进行基因survey（持续更新中） 这是本人自己测序的二代数据，公司提供的是rawdata。文库350bp。测序类型PE150。二代数据特点就是：短读长，低质量。因此第一步就是过滤数据，得到cleandata。 本人的坑：不管是处理什么数据，第一步都是观察数据，了解数据的来源，构造，是否有菌污染等。通常菌污染是进行blast，与nt库进行比对（可查看本人的另一篇文章）。 第一步：过滤数据（filter） 指标：①过滤N碱基占比>10%的reads；②过滤低质量碱基（质量低于<5）占比
数据质量什么叫做好，以及Trimmomatic sickle等过滤软件的比较
sunwanying123的博客
07-04 2638
在二代测序合成的过程中随着合成链的增长，DNA聚合酶的效率会不断下降，特异性也开始变差，这就会带来一个问题——越到后面碱基合成的错误率就会越高，以下几个指标用于观察数据的质量情况read各个位置的碱基质量值分布 ：大于30且波动较小碱基的总体质量值分布：大部分高于20read各个位置上碱基分布比例，目的是为了分析碱基的分离程度，AT CG在1%以内GC含量分布read各位置的N含量，一般不该出现r...
从fastq文件中批量提取/过滤序列【python】
每天都要学Python的博客
03-09 9557
博主也是刚刚接触生信，会将自己平时练习用到的python脚本发布到博客上，用来记录自己的学习之路。 介绍 测序回来的fastq文件通常在分析之前，需要进行过滤，该脚本利用python实现从压缩的fastq文件中提取指定ID的序列，并保存为新的压缩格式的fastq文件。 说明 输入文件为fq.gz文件，压缩的ID list文件。 必须是压缩格式的文件才可以，如果非压缩格式，可以压缩成gz格式后...
linux fastQC 操作命令,Linux shell合并fastq测序数据/批量fastqc小脚本|merge|multiqc
weixin_42470362的博客
04-29 3387
合并fastq测序数据不同泳道的同一个样品测序数据经过质量检查QC后是可以合并的。本例中文件命名情况如下：示例文件名：83b_S156_L004_R1_001.fastq.gz，其中83b_S156是样品名，L004是泳道，R1是正向序列。这时候只需要用这个小脚本即可：ls *R1* | cut -d _ -f 1,2 | sort | uniq \| while read id; do \cat...
生信技能26 - BWA比对、去除PCR重复和统计最佳实践
LittleComputerRobot的博客
08-14 234
BWA MEM 比对、去除PCR重复和统计最佳实践
PCR应用运行分析软件
01-15
PCR原始数据应用软件
PCR详细讲解PPT课件.ppt
04-26
PCR详细讲解，PCR详细讲解课件，PCR详细讲解PPT
PCR.zip_pcr_远程控制
09-22
这是一个高级远程控制程序的源代码，可以通过本机操纵对方鼠标，所有传输均通过压缩，确保速度
荧光定量PCR数据分析.doc
11-30
荧光定量PCR数据分析.doc
基因组重测序全流程(简易版)
weixin_73362123的博客
11-23 2649
基因组重测序全流程(简易版)
iMeta | 刘金鑫/徐珂琳/张安云-揭示生命早期微生物碳水化合物代谢的改变
最新发布
刘永鑫的博客——宏基因组公众号
03-26 97
点击蓝字 关注我们婴儿年龄与肠道携带抗性基因呈现负相关，揭示生命早期微生物碳水化合物代谢的改变iMeta主页：http://www.imeta.science研究论文●原文链接DOI: https://doi.org/10.1002/imt2.169●2024年1月31日，南京农业大学刘金鑫、复旦大学公共卫生学院徐珂琳和四川大学生命科学学院张安云团队在iMeta在线联合发表了题为 “Infan...
bam文件去除pcr重复和去除线粒体信息
m0_70611557的博客
04-26 554
sambamba markdup -r input.bam output.bamsamtools view -h input.bam |grep -v "chrM" |samtools view -bS -o output.bam
2021.12.04【宏基因组】|宏基因组流程搭建进展梳理
12-06 1204
分析框架 质控 过滤 trimmomatic PE {input.R1} {input.R2} {output.R1_PE} {output.R1_UNPE} {output.R2_PE} {output.R2_UNPE} LEADING:3 TRAILING:3 SLIDINGWINDOW:5:20 MINLEN:50 -phred33 去重 echo \"{input.R1}\" > {params.sample_ID}_merge.fastuniq ;ech..
linux uniq命令_如何在Linux上使用uniq命令
cum43546的博客
10-09 757
linux uniq命令Fatmawati Achmad Zaenuri/ShutterstockFatmawati Achmad Zaenuri / ShutterstockThe Linux uniq command whips through your text files looking for unique or duplicate lines. In this guide, we co...
如何去除bam 文件中的pcr 重复
03-31
其中，MarkDuplicates 工具可以检测和去除 PCR 重复。该工具可以标记 PCR 重复的序列，也可以将它们从 BAM 文件中删除。 安装 Picard 工具： ``` conda install -c bioconda picard ``` 使用 MarkDuplicates 工具...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值