当前位置: 首页 > news >正文

怎么用dw英文版做网站宁波网站制作公司

怎么用dw英文版做网站,宁波网站制作公司,常州做网站设计,网站设计模版免费建站原文链接:如何使用bedtools、convert2bed、gff2bed提取基因序列 前言 我们的使用基因组注释文件gtf或gff文件从基因组fa文件中提取transcript的方式很多,相对用的比较多的是使用gffread软件。但是gffread软件,提取的序列一般都是transcript序列&#xf…

原文链接:如何使用bedtools、convert2bed、gff2bed提取基因序列

前言

我们的使用基因组注释文件gtfgff文件从基因组fa文件中提取transcript的方式很多,相对用的比较多的是使用gffread软件。但是gffread软件,提取的序列一般都是transcript序列,若是,我们的想提取的gene序列,那么不能直接使用基因组gtf文件。

本次,我们介绍其中的一种方法。使用bedtools getfasta提取gene全长序列。

注意:这仅仅只是其中的一种方法而已。

软件安装

  1. 软件安装,使用mamba安装bedtools,convert2bed,gff2bed
mamba install -y bedtools convert2bed gff2bed
  1. 软件测试
bedtools -h
$ bedtools -hbedtools is a powerful toolset for genome arithmetic.Version:   v2.31.1
About:     developed in the quinlanlab.org and by many contributors worldwide.
Docs:      http://bedtools.readthedocs.io/
Code:      https://github.com/arq5x/bedtools2
Mail:      https://groups.google.com/forum/#!forum/bedtools-discussUsage:     bedtools <subcommand> [options]The bedtools sub-commands include:[ Genome arithmetic ]intersect     Find overlapping intervals in various ways.window        Find overlapping intervals within a window around an interval.closest       Find the closest, potentially non-overlapping interval.coverage      Compute the coverage over defined intervals.map           Apply a function to a column for each overlapping interval.genomecov     Compute the coverage over an entire genome.merge         Combine overlapping/nearby intervals into a single interval.cluster       Cluster (but don't merge) overlapping/nearby intervals.complement    Extract intervals _not_ represented by an interval file.shift         Adjust the position of intervals.subtract      Remove intervals based on overlaps b/w two files.
  1. 制作bed文件

使用bedtools提取序列,需要制作bed文件。格式如下所示:

我们可以直接使用awk命令进行提取对应的信息。我们使用Cucumber的gff注释文件提取。

cat Cucumber.CLv4.gff3 | awk '{if($3 == "gene") print $0}' | awk '{print $1"\t"$4"\t"$5"\t"$9"\t"$7}' | head

提取对应的列信息即可。


使用convert2bed或gff2bed结合起来提取信息。

cat Cucumber.CLv4.gff3 |awk '{if($3~/^gene$/)print }' > 01.gene.gff && convert2bed --input=gff --output=bed < 01.gene.gff > 02.gene.bed

获得如下信息

我们在此基础上加上awk的命令即可,同上。

cat Cucumber.CLv4.gff3 |awk '{if($3~/^gene$/)print }' > 01.gff && convert2bed --input=gff --output=bed < 01.gff > 02.bed && awk '{print $1"\t"$2"\t"$3"\t"$10"\t"$6}' <02.bed>03.bed

使用gff2bed软件结合。

awk '{if($3~/^gene$/)print}' Cucumber.CLv4.gff3  > 01.genes.gff && gff2bed <01.genes.gff> 02.genes.bed 

使用bedtools getfasta提取序列。

bedtools getfasta -fi Cucumber.geome.fa -bed 02.gene.bed -fo cucumber.gene.fa -name -s 
  • -name参数是必须加的,若是不加,你的cucumber.gene.fa文件中无基因名。
$ bedtools getfasta -hTool:    bedtools getfasta (aka fastaFromBed)
Version: v2.31.1
Summary: Extract DNA sequences from a fasta file based on feature coordinates.Usage:   bedtools getfasta [OPTIONS] -fi <fasta> -bed <bed/gff/vcf>Options: -fi		Input FASTA file-fo		Output file (opt., default is STDOUT-bed		BED/GFF/VCF file of ranges to extract from -fi-name		Use the name field and coordinates for the FASTA header-name+		(deprecated) Use the name field and coordinates for the FASTA header-nameOnly	Use the name field for the FASTA header-split		Given BED12 fmt., extract and concatenate the sequencesfrom the BED "blocks" (e.g., exons)-tab		Write output in TAB delimited format.-bedOut		Report extract sequences in a tab-delimited BED format instead of in FASTA format.- Default is FASTA format.-s		Force strandedness. If the feature occupies the antisense,strand, the sequence will be reverse complemented.- By default, strand information is ignored.-fullHeader	Use full fasta header.- By default, only the word before the first space or tab is used.-rna	The FASTA is RNA not DNA. Reverse complementation handled accordingly.

输出的结果基因ID如下所示。

>CsaV4_1G000004::chr1:1088370-1092905
>CsaV4_1G000005::chr1:1095847-1098019
>CsaV4_1G000003::chr1:1084077-1087157

使用sed命令进行批量处理:

sed 's/::.*//' input.fa > output.fa


若我们的教程对你有所帮助,请点赞+收藏+转发,大家的支持是我们更新的动力!!


2024已离你我而去,2025加油!!

2024年推文汇总 (点击后访问)

2023年推文汇总 (点击后访问)

2022年推文汇总 (点击后访问)

往期部分文章

1. 最全WGCNA教程(替换数据即可出全部结果与图形)

  • WGCNA分析代码六

推荐大家购买最新的教程,若是已经购买以前WGNCA教程的同学,可以在对应教程留言,即可获得最新的教程。(注:此教程也仅基于自己理解,不仅局限于此,难免有不恰当地方,请结合自己需求,进行改动。)


2. 精美图形绘制教程

  • 精美图形绘制教程
  • 《R语言绘图专栏–50+图形绘制教程》

3. 转录组分析教程

  • 转录组上游分析教程[零基础]

  • 一个转录组上游分析流程 | Hisat2-Stringtie

  • Samll RNA上游分析

4. 转录组下游分析

  • 批量做差异分析及图形绘制 | 基于DESeq2差异分析

  • GO和KEGG富集分析

  • 单基因GSEA富集分析

  • 全基因集GSEA富集分析

BioinfoR生信筆記 ,注于分享生物信息学相关知识和R语言绘图教程。

http://www.dtcms.com/a/456395.html

相关文章:

  • 个人免费发布信息胶州网站优化
  • 网站开发年收入做网站要学什么东西
  • 沈阳手机端建站模板阿里云网站301重定向怎么做
  • 南昌网站建设冲浪者如何制作网站地图
  • dede学校网站免费源码有哪些官网做的比较好的网站
  • wordpress水煮鱼网站seo公司哪家好
  • 如何做招聘网站的方案新闻发稿渠道
  • 购买了个网站源码 怎么建立建设银行网站无法打开
  • 福州网站建设熊掌号seo3立体构型
  • 电子商务网站会员体系比较顺口的建筑公司名字
  • 网站设计哪家强免费手机网页网站
  • 电商网站建设任务分解结构官方网站建设项目询价
  • 重庆有哪些网站大连网站建设价格
  • 正保建设工程网站重庆巨能建设集团网站
  • 佛山专业网站设计方案免费做店招哪个网站好
  • 网站背景色智能优化网站
  • 网站开发的可行性报告济南资海网站建设公司
  • 自己可以免费做网站吗wordpress主题xin
  • 深圳设计网站有哪些深圳网站建设服务哪个便宜点
  • 青岛网站推广公司排名哪些网站做任务好赚钱
  • 微信网站 顶部导航菜单信息发布型企业网站的特点
  • qq空间钓鱼网站后台怎么做做宣传页的软件是什么
  • 做胃镜多少钱天津津门网站I工会网站建设
  • seo网站结构优化的方法wordpress php函数大全
  • 织梦移动网站专业网站建设
  • 石家庄微网站建设公司哪家好app平台下载
  • 一个ip 做2个网站在线制作名片免费
  • 沈阳优化网站市场营销管理
  • 国内装饰行业网站制作工程项目管理软件排名
  • 建站中心购买模板做网站