怎么根据class_code筛选转录本

这篇文章主要介绍“怎么根据class_code筛选转录本”，在日常操作中，相信很多人在怎么根据class_code筛选转录本问题上存在疑惑，小编查阅了各式资料，整理出简单好用的操作方法，希望对大家解答”怎么根据class_code筛选转录本”的疑惑有所帮助！接下来，请跟着小编一起来学习吧！

10年积累的网站建设、网站设计经验，可以快速应对客户对网站的新想法和需求。提供各种问题对应的解决方案。让选择我们的客户得到更好、更有力的网络服务。我虽然不认识你，你也不认识我。但先网站设计后付款的网站建设流程，更有镇雄免费网站建设让你可以放心的选择与我们合作。

问题描述

链特异文库鉴定长链非编码RNA(lncRNA)的基本步骤是

hisat2将原始测序数据比对到参考基因组
samtools获得排序的bam文件
stringtie每个样本分别组装得到转录本，获得是一个gtf文件

stringtie -p 12 --rf -o transcripts.gtf sorted.bam

gffcompare合并所有样本的gtf文件

ls *.gtf > list.merged.txt
~/Biotools/gffcompare/gffcompare -r reference.gtf -i list.merged.txt -o merged

得到一个 merged.combined.gtf这个文件里给每一个转录本分配了一个class_code用来表示转录本相对于参考基因组的位置怎么根据class_code筛选转录本

以上图片来源于论文 GFF Utilities: GffRead and GffCompare

长链非编码RNA通常是选择class_code为u/x/i，比如论文 Global identification of Arabidopsis lncRNAs revealsthe regulation of MAF4 by a natural antisense RNA中就提到 only transcripts with TAIR10 annotation [Cufflinks class codes ‘u’ (intergenic transcripts),’x’ (Exonic overlap with reference on the opposite strand),’i’ (transcripts entirely within intron) were retained.

那么问题就来了，如何利用 merged.combined.gtf这个文件获得 class_code 为 u、x和i的转录本的gtf文件呢

找到了一个办法，python中有一个模块 pyGTF，github链接是https://github.com/chengcz/pyGTF

直接使用pip安装

pip install pyGTF

可以解析gft格式的注释文件

利用这个模块来写一个简单的脚本

import sys
from pyGTF import Transcript
from pyGTF import GTFReader

in_gft = sys.argv[1]
class_code = sys.argv[2]
out_gtf = sys.argv[3]

fw = open(out_gtf,'w')

with GTFReader(in_gtf,flag_stream=True) as fi:
 for i in fi:
  if i._attri['class_code'] == class_code:
   i.to_gtf(fw)
   
fw.close()

使用方法是

python 01.py in.gtf i out.gtf

####今天学到的另外一个知识点： samtools统计fasta文件序列长度，根据序列名提取序列

参考

https://www.cnblogs.com/xudongliang/p/5200655.html

使用命令

samtools faidx input.fasta

会生成一个input.fasta.fai的文件，文件的内容总共有5列第一列是序列名，第二列是序列长度，第四列是每行多少个碱基

根据序列名提取序列这里好像只能提取单条序列

samtools faidx input.fasta TCONS_00000018 > TCONS_00000018.fa

还可以加上指定的位置

samtools faidx input.fasta TCONS_00000018:1-10
>TCONS_00000018:1-10
TGGGCGAACG

到此，关于“怎么根据class_code筛选转录本”的学习就结束了，希望能够解决大家的疑惑。理论与实践的搭配能更好的帮助大家学习，快去试试吧！若想继续学习更多相关知识，请继续关注创新互联网站，小编会继续努力为大家带来更多实用的文章！

新闻名称：怎么根据class_code筛选转录本
浏览地址：http://dzwzjz.com/article/pcghsd.html

大橙子网站建设，新征程启航

怎么根据class_code筛选转录本

其他资讯