QIIME pipeline 采用 fastq-join 和 SeqPrep 合并Illumina paired-end reads。其安装方法如下: Currently, there are two methods that can be selected by the user to join paired-end data: fastq-join - Erik Aronesty, 2011. ea-utils : “Command-line tools for processing biological sequencing data” ( ...
本人是科学网新用户,今天在写博文时想对不同段落设置不同行距。在工具栏找了好久还是没能找到有效调整不同行距的按钮,冒昧在此向各位博友求助!谢谢! 例如,如何减小下面英文部分行间距让它看起来更紧凑? Spend all your time waiting for that second chance For a break that would make it OK There’s ...
# Written and tested in Python 2.7.2 and bipython #Extract_annotation_using_BlastSubjectIds(fast).py # 本脚本用来从 blast 结果(m6 fomat)中,根据第二列subject id 从原始database中抽取相关注释及参考序列。 import os from Bio import SeqIO while True: Parameters=raw_input( ...
最近从基于高通量测序的16S分析转换到 环境样品的 Metagenome 分析时,需要统计组装后contig的GC含量和N50。于是花了点时间写了个小脚本,方便今后数据处理工作。 #GC_N50.py print 'Python and Biopython needed for running this script!' print Script for calculating N50 of assembly ...