标签：剪切 junction read isoform reads 可视化 exon chr19 bam

根据bulk RNA-seq或者single-cell RNA-seq来找isoform的类型。

 

难点：二代NGS的reads长度较短，可能无法成功组装出full-length的isoform。

 

我们只需要关注junction read，基数即可。

 

junction read：在成熟mRNA中，测序的reads如果同时跨越了两个exon，则为junction read。

如何提取出junction read？

 

/home/lizhixin/softwares/regtools/build/regtools junctions extract -r chr19:797075-812327 -s 0 UE.bam

　　

自己写一个程序吧，用bedtools coverage这个功能。

bedtools coverage -a regions.bed -b reads.bam

　　

spliceGraph

 

核心就是count reads，同时比对到两个点上。【这个不难，但是并不能说明问题】

samtools view -u -@ 2 UE.bam chr19:803636-803636 | less -S

　　

 必须鉴定出拼接两个exon的read，这样的才是有意义的证据，单纯的连接的数据并不实用【中间可能插入了很多个exon，连接≠拼接】。

 

cat gencode.v27.annotation.gtf | grep "\"PKM\"" | awk '{if($3=="exon")print $0}' > PKM.txt

　　

 

samtools view -b -@ 2 23c9.bam chr19:797075-812327 > PTBP1.23c9.bam

　　

 

参考：

Question: Count Junctions In A Sam/Bam File

junctions extract - RegTools

Question: How Many Reads In A Bam File Are Aligned To a Specific Region

 

标签：剪切,junction,read,isoform,reads,可视化,exon,chr19,bam
来源： https://www.cnblogs.com/leezx/p/14355867.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。