04 2019 档案
python学习——读取染色体长度（七：读取fasta文件）
摘要：读取fasta文件genome_test.fa，并计算染色体总长，同时输出最长染色体编号、序列以及长度 fasta文件genom_test.fa的内容如下： >chr1ATATATATAT>chr2ATATATATATCGCGCGCGCG>chr3ATATATATATCGCGCGCGCGATATAT
阅读全文
python学习——读取染色体长度（六：读取含有染色体长度的文件）
摘要：含有染色体长的文件chr_len.txt chr1 10chr2 20chr3 30chr4 40chr5 50 python脚本 #传递命令行参数 import sys # 导入模块 # 从命令行获取文件名称 f_chr_len = sys.argv[1] # 定义命令行参数，1表示变量1 # 打
阅读全文
python学习——读取染色体长度（五：从命令行输入染色体长度）
摘要：# 传递命令行参数 # 导入sys模块 import sys print(sys.argv) 命令行操作 python argv.py 10 20 30 40 50 回车输出 ['argv.py', '10', '20', '30', '40', '50'] # 方法一，删掉第0个元素 del(sy
阅读全文
python学习——读取染色体长度（四：获取最长染色体的编号）
摘要：# 读取fasta # 解析每条序列的长度 chr_len = {'chr1':10,'chr2':20,'chr3':30,'chr4':40,'chr5':15} # 求和 total_len = sum(chr_len.values()) # 获取键用keys for chr in chr_l
阅读全文
python学习——读取染色体长度（三、用循环或者函数求总长并获取最长染色体长度）
摘要：# 读取fasta # 解析每条序列的长度 chr_len = [10,20,30,40,50] # 求和 # 方法一：通过循环 total_len = 0 #定义total_len的初始长度 for len in chr_len: # 从列表chr_len中每次取一个值交给len total_le
阅读全文
python学习——读取染色体长度（二、向前一步：通过染色体序列获得长度信息）
摘要：# 读取fasta # 解析每条序列的长度 chr1_seq = 'ATATATATAT' chr2_seq = 'ATATATATATCGCGCGCGCG' chr3_seq = 'ATATATATATCGCGCGCGCGATATATATAT' chr4_seq = 'ATATATATATCGCG
阅读全文
python学习——读取染色体长度（一、简化问题）
摘要：# 读取fasta # 解析每条序列的长度 chr1_len = 10 chr2_len = 20 chr3_len = 30 chr4_len = 40 chr5_len = 50 # 求和 total_len = chr1_len + chr2_len + chr3_len + chr4_len
阅读全文
R语言学习——数据合并及绘制密度分布曲线图
摘要：setwd("E:/08_cooperation/07_X-lab/06-Crosstalk/Aadapter_primer")# 读取lane01.txt,并对其按列进行相加处理，然后对列进行进行命名d1=read.table("lane01.txt",header = FALSE,sep = "
阅读全文
