摘要： 转载自：https://blog.csdn.net/soyabean555999/article/details/62235577 一、转录组还是基因组？ map常用的工具有bowtie/bowtie2, BWA,SOAP1/SOAP2等。这个问题又会被分成两个问题，是基因组测序（DNA-seq）还 阅读全文

摘要： 想实现对多层文件中的某些指定数据的处理 如下： 一级目录 二级目录 三级目录 对所有lane文件夹下的所有cyc文件夹中的所有R001C002图像进行处理 可以使用父bat与子bat文件进行处理 子bat文件（test5.bat）代码如下： @echo off setlocal enabledela 阅读全文

摘要： 转自：https://www.cnblogs.com/followyourheart/articles/3349899.html 统计指标是数据分析的基本元素，变量之间的对比分析和综合分析是最基本、最常用的统计分析方法。当统计指标的量纲不同或性质不同时，如果直接用原始数据进行数据分析，往往会得到不合 阅读全文

摘要： 数据示例 代码实例 阅读全文

摘要： 提取fasta文件genome_test.fa中第14号染色体的序列，其内容如下： 用python以及命令行参数实现 新建.py文件“”GetSeqFromChrID.py”， python脚本如下： 命令行参数输入如下：红色字体是输入部分 结果如下： 阅读全文

摘要： 将读取文件的代码封装成函数，并使其作为模块可在其他程序运行 创建fasta_def.py文件，并输入如下代码： 新建一个test.py文件 输入如下代码调用上述模块fasta_def 阅读全文

摘要： 把代码封装成函数的好处是可以重复使用该段代码，并且会使代码结构清晰 例如要计算chr1以及chr2染色体的GC含量，代码如下： 阅读全文

摘要： argparse是python中用于传递和解析命令行参数的模块，例如： 从cmd输入命令行 E:\15_python\DEBUG>python fasta_argparse.py -avg -min_len 50 100 10 200 40 回车得到 total length:300.00avera 阅读全文

摘要： 以人类染色体个数为例： 通过range()函数快速生成一系列整数，将其一列表的形式存储；对其进行扩展；然后修改列表中的元素（染色体id）,生成Chr_形式 cmd终端python代码如下： 用VSCode编写python 阅读全文

摘要： 对fasta文件genome_test.fa中的染色体序列进行反向互补，并输出到文件genome_test_RC.fa genome_test.fa >chr1ATATATATAT>chr2ATATATATATCGCGCGCGCG>chr3ATATATATATCGCGCGCGCGATATATATAT 阅读全文