文档章节

回贴基因组

xiechaoyong
 xiechaoyong
发布于 2014/02/28 18:40
字数 291
阅读 245
收藏 1
for org in `head org.v4.1 | cut -f 2`; do 
cd $org;
bash /leofs/noncode/NONCODEv4/cmds/blat.sh  $org.2bit v3.fa
cd ..;
echo $org; done

回贴基因组流程:

代码如下:

blat.sh:

twoBitDb=$1
queryFa=$2

blat $twoBitDb $queryFa $queryFa.psl
$CMD/pslReps -singleHit $queryFa.psl $queryFa.psl.1 $queryFa.psr
$CMD/pslToBed  $queryFa.psl.1 $queryFa.bed
perl /leofs/noncode/NONCODEv4/cmds/dup.pl $queryFa.bed 3 > $queryFa.bed.dup
bedToGtf.sh $queryFa.bed.dup > $queryFa.gtf
cuffcompare -C -o $queryFa.c $queryFa.gtf

cut -f 4 $queryFa.bed > ncid.Get.Blat1
perl $CMD/rmFa.pl $queryFa ncid.Get.Blat1 > $queryFa.2

blat -minScore=0 -minIdentity=90 $twoBitDb $queryFa.2 $queryFa.psl.2
$CMD/pslReps -singleHit $queryFa.psl.2 $queryFa.psl.1.2 $queryFa.psr.2
$CMD/pslToBed  $queryFa.psl.1.2 $queryFa.bed.2
perl /leofs/noncode/NONCODEv4/cmds/dup.pl $queryFa.bed.2 3 > $queryFa.bed.dup.2
bedToGtf.sh $queryFa.bed.dup.2 > $queryFa.gtf.2
cuffcompare -C -o $queryFa.c.2 $queryFa.gtf.2

gtf2Bed.pl $queryFa.c.combined.gtf > $queryFa.c.bed
gtf2Bed.pl $queryFa.c.2.combined.gtf > $queryFa.c.2.bed
perl -ne '
@t = split(/\t/,$_);
$tid = $t[0]; 
$t[4] =~ /.*?\|(.*?)\|/;
$ncid = $1;
print $tid . "\t" . $ncid . "\n"; 
' $queryFa.c.tracking > idMapping.1

perl -ne ' 
@t = split(/\t/,$_);
$tid = $t[0];
$t[4] =~ /.*?\|(.*?)\|/;
$ncid = $1;
print $tid . "\t" . $ncid . "\n";
' $queryFa.c.2.tracking > idMapping.2




####!得到NONCODE的Gtf, 然后与其它物种做cuffcompare, 将c的合并掉。
#cuffcompare -r $queryFa.gtf -o v4 $queryFa.gtf 1.gtf 2.gtf

perl $CMD/leftJoin.pl $queryFa.c.bed 4 idMapping.1 1 | awk '{print $1"\t"$2"\t"$3"\t"$14"\t"$5"\t"$6"\t"$7"\t"$8"\t"$9"\t"$10"\t"$11"\t"$12;}' > $queryFa.blat.bed
perl $CMD/leftJoin.pl $queryFa.c.2.bed 4 idMapping.2 1 | awk '{print $1"\t"$2"\t"$3"\t"$14"\t"$5"\t"$6"\t"$7"\t"$8"\t"$9"\t"$10"\t"$11"\t"$12;}' >> $queryFa.blat.bed


© 著作权归作者所有

下一篇: 生成v3.fa
xiechaoyong
粉丝 1
博文 27
码字总数 3531
作品 0
海淀
私信 提问
别再叫它“垃圾”了——这个“跳跃”基因可能就是胚胎发育的关键

  在你生命中的某个时刻,你还是个双细胞胚胎,正翻着跟头进行子宫“软着陆”。一项新的研究表明,能完成这个基础过程还要多亏了一种 DNA,而该 DNA 曾一度被诸多科学家当作“垃圾”。   ...

DeepTech深科技
2018/06/25
0
0
炸锅了!美国正式宣布人工智能+基因科学可以让人类永生!

  这是改变世界的一次握手!   世界上两大最尖端的科学:人工智能与基因科学,终于胜利会师!   刚刚,谷歌向世界宣布:谷歌拿出了一款取名为DeepVariant的人工智能,它能运用最先进的...

人工智能机器人联盟
2017/12/14
0
0
谷歌开源DeepVariant,之前的经典检测基因变异法将被颠覆

  【新智元导读】今天,谷歌开源了DeepVariant深度学习模型,这是将基因组中的识别变异(variant calling)任务转换为图像分类问题的基因组学工具,在重构基因组序列方面比以前的经典方法具...

深度学习
2017/12/06
0
0
迄今最连续人类基因组组装完成

英国《自然·生物技术》杂志29日在线发表了一项基因组学研究成果:英国科学家团队利用一个口袋大小的纳米孔测序装置,从头组装出人类基因组,这是迄今为止最连续的人类基因组组装,而且只用了...

软件达人
2018/01/30
0
0
“再生能力”背后的遗传学基础揭示

英国《自然》杂志近日发表了两篇基因学论文,欧洲两组团队分别报告了美西螈和真涡虫的基因组,揭示了神秘“再生能力”背后的遗传学基础。其中美西螈的320亿个碱基对,是目前组装出的最大基因...

软件达人
2018/01/29
0
0

没有更多内容

加载失败,请刷新页面

加载更多

nproc systemd on CentOS 7

Increasing nproc for processes launched by systemd on CentOS 7 Ask Question I have successfully increased the nofile and nproc value for the local users, but I couldn't find a p......

MtrS
今天
3
0
了解微信小程序下拉刷新功能

小程序提供了这个事件。 onPullDownRefresh() 监听用户下拉刷新事件。 如果要开启下拉刷新功能,要先到json配置: "enablePullDownRefresh":true 配置后下拉有反应了但是没有加载效果,在onP...

oixxan__
今天
2
0
springmvc java对象转json,上传下载(未完)拦截器Interceptor以及源码解析(未完待续)

package com.atguigu.my.controller;import java.util.Collection;import org.springframework.beans.factory.annotation.Autowired;import org.springframework.stereotype.Contr......

architect刘源源
今天
29
0
[日更-2019.5.24、25、26] Android系统中的Binder通信机制分析(一)--servicemanager

声明 其实对于Android系统Binder通信的机制早就有分析的想法,记得去年6、7月份Mr.Deng离职期间约定一起对其进行研究的,但因为我个人问题没能实施这个计划,留下些许遗憾... 最近,刚好在做...

Captain_小馬佩德罗
昨天
24
0
聊聊dubbo的DataStore

序 本文主要研究一下dubbo的DataStore DataStore dubbo-2.7.2/dubbo-common/src/main/java/org/apache/dubbo/common/store/DataStore.java @SPI("simple")public interface DataStore { ......

go4it
昨天
3
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部