生信脚本练习（9）合并文件 ②

最新推荐文章于 2024-03-10 21:27:23 发布
最新推荐文章于 2024-03-10 21:27:23 发布
阅读量586 收藏
点赞数
分类专栏： python linux bioinfo 文章标签： 脚本
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hxoxh/article/details/76856103
版权
这篇博客介绍了如何合并两个文件，通过将样本信息转化为字典，利用匹配的键替换值为'-'的方法，展示了巧妙的脚本处理技巧。
摘要由CSDN通过智能技术生成

这个练习也是合并文件，但是只是合并两个文件

这是样本信息：

MegaID  Chr Star    Ref Alt Depth
170602C304  1   861302  G   A   635
170602C304  1   906303  G   T   290
170602C304  1   985841  C   T   56
170602C304  1   985866  A   G   376
170602C304  1   1120344 C   T   439
170602C304  2   21255445    G   A   1009

以下是带注释的信息，要求把注释补充到样本信息文件的后面

Chr Star    End Ref Alt Cosmic
1   69345   69345   C   A   ID=COSM911918;OCCURENCE=1(endometrium)
最低0.47元/天 解锁文章
Hao222222222222
关注
专栏目录
生信脚本练习（13）文件合并（初识numpy）
49
08-20 638
这本质上还是个文件合并的问题这是sample文件 这是database文件 要求是把database文件的第四行追加到sample文件后面。限制条件： 区间要有重合。也就是说start和end这两列的数字表示每个样本的区间 变异类型要一致，也就是type这一行要一样。 最后做出来要像这样: 最后追加一行 有符合限制条件的就追加上去，逗号隔开，没有就空着 # 循环遍历的判定方法 import
生信脚本练习（6） 求read每个位点的cg分布
49
08-07 712
这是另外一题，也是cg含量， 但是是求一个fastq文件中每个read位点的cg分布情况 这个解法是逐行读取文件，比较优雅number = {} buffer = 200 for i in range(buffer): kkk = i number[kkk] = 0 with open("Test1.fastq","r") as f: li
生信分析脚本
04-16
RNA数据格式转换脚本 在linux系统上使用python进行RNA数据格式转换
生信常用脚本
weixin_30674525的博客
09-20 844
1.用自己的gtf文件构建txdb数据结构，用于下游分析(The TxDb class is a container for storing transcript annotations) library(GenomicFeatures) #gtf、gff to txdb makeTxDbFromGFF() #other way makeTxDbFromGRanges makeT...
生信自动化流程搭建 02 | 脚本
Baimoc
08-10 1924
脚本 脚本是字符串声明，它定义了由过程执行到执行任务的命令。 一个进程仅包含一个脚本块，并且当该进程包含输入和输出声明时，它必须是最后一个语句。 输入的字符串在主机系统中作为Bash脚本执行。它可以是通常在终端 shell 程序或通用Bash脚本中使用的任何命令，脚本或它们的组合。 可以在脚本语句中使用的命令的唯一限制是目标执行系统中这些程序的可用性。 脚本可以是简单字符串或多行字符串，例如： process doMoreThings { """ blastp -db $db -query que
生信脚本练习（12）求fasta文件各序列长度并统计作图
49
08-20 6378
题目要求是要从一个fasta文件中统计出每条序列的长度分布，并作图。代码如下：import os import getpass import matplotlib.pyplot as plt usr = getpass.getuser() os.chdir('c:/Users/' + usr + '/Desktop') seq_len = {} # 把fasta文件全部读取做成字典，键是带‘>’的那
生信分析服务器配置脚本
08-13
生信分析linux环境配置，安装分析软件环境，需要root权限
生信脚本练习（7）求fastq文件质量值分布
49
08-07 926
这里是求fastq文件的质量类型和比例 注意数几种类型用了set() 的方法with open("Test1.fastq","r") as f: lines = f.readlines() qual = [] for line in lines: if lines.index(line)% 4 == 3: line= line.strip('\n')
生信脚本练习（5）求fastq文件的cg含量
49
08-07 1397
这道题要求求出一个fastq文件中序列的cg含量with open("Test1.fastq","r") as f: lines = f.readlines() seq = [] index = [] cg_n = [] cg_percentage = [] head = ['Read', 'CG'] cg = open("cg.txt","w") #这里精髓，取出每条序列的标题 for
Perl-for-Bioinformatics：尝试帮助有兴趣将Perl用于生物信息学的任何人
02-05
用Perl编写的生物信息学工具 这里的大多数脚本是在我从事不同项目时编写的，我认为这对其他人将很有用，并且可以根据需要进行扩展/修改。 IO ::常规 脚本使用自定义 Perl模块。 请通过浏览目录查看安装说明。 如果要安装lncRNApipe Pipeline，则会自动安装IO::Routine模块。 需要Bio::SeqIO模块已安装且可用。 nc lncRNA管道 从头开始提取推定的新型lncRNA的管道，其中提供了从深度测序数据（例如：RNA-Seq）和注释数据组装而成的GTF格式的转录本列表。 转到目录以获取脚本列表。 安装lncRNApipe及其所有依赖项（Mac和
Missing semicolon设置方法（去波浪线！）
08-05
Missing semicolon设置方法（去波浪线！）
R操作合并数据等
weixin_42480153的博客
06-09 224
library(data.table) setwd("/mnt/10t/mzy/48samples/16s/salmon/") fread("DNAabd138.csv")->DNA as.data.frame(DNA)->DNA row.names(DNA)<-DNA[,1] unique(DNA[,1])->namelist x<-c("") for (i in namelist){ t(data.frame(colSums(DNA[which(DNA$Name==i),3
生信脚本练习（8）合并文件 ①
49
08-07 678
这道题有难度，要把这三个文件合并到一起。文件一 1 161514631 T C|0.132632|(272,140,47,16) 2 222301193 A C|0.078624|(260,115,23,9) 3 89259567 A C|0.043716|(200,150,8,8) 4 55979552 G A|0.211921|(236,12
Annovar软件注释肿瘤基因突变：COSMIC数据库最新版下载与使用
最新发布
公众号/简说基因，知乎/简宝玉
03-10 482
COSMIC，即：Catalogue Of Somatic Mutations In Cancer。官网：https://cancer.sanger.ac.uk/cosmic注意事项：Cosmic 数据库从 v71 版本开始，学术研究免费（需要学校或研究机构的邮箱注册），而商业使用需要授权目前数据库有基于 CRCh37 和 CRCh38 两种版本，根据需要选择合适的，步骤：进入主页 > 找到...
BGI-College生信入门——6、Perl小骆驼学习笔记（一）
bgicollege的博客
07-19 369
Perl 入门
missing semicolon报错问题解决
热门推荐
Taurus_666
09-02 2万+
missing semicolon报错主要是因为eslint检测出现错误， 一般是看下提示行是否写分号，添加上就没问题了
生信技能树——GEO芯片数据的合并
weixin_48275332的博客
04-13 1万+
GSE83521和GSE89143数据合并 1.下载数据 rm(list = ls()) library(GEOquery) library(stringr) gse = "GSE83521" eSet1 <- getGEO("GSE83521", destdir = '.', getGPL = F) eSet2 <- getGEO("GSE89143", destdir = '.',
生信分析常用脚本（一）
weixin_30840573的博客
11-02 1821
1.计算reads PE长度脚本 01.Cal_PE_Depth.sh Read_count=`gzip -dc Reads1.fq.gz |wc -l |awk '{print $1/4}'` echo "read pair count : $Read_count" average_depth=`expr $Read_count \* 200 / 6000000 ` echo av...
生信自学网geo合并代码
09-26
GEO（Gene Expression ...通过生信自学网提供的GEO合并代码示例，我们可以方便地将来自GEO的数据集进行合并和分析，从而深入研究基因表达的模式和生物学机制。这对于疾病的研究以及药物研发等方面具有重要的意义。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值