人类、小鼠、新冠病毒(SARS-CoV-2)参考序列下载

为了学习jimmy大神的视频,一直在努力重复他视频里的每一个步骤,结果linux都快学完了,才发现😅

image.png

哎,反正我已经折腾的差不多了,就这样吧,jimmy大神的东西真的很全,也很用心,但也正因为太全了,太用心了,每到一个点就恨不得把相关的都告诉你,再加上我是按照他的他视频的顺序(第一个为‘生信人应该这样学linux’,弹幕也有吐槽逻辑有时候有些混乱,后来发现还要学点他的别的视频)学习,并且这是他很早以前录制的,可能那时候还没有完全弄好学习的顺序,尤其是对于小白,我就一步一步为了实现他视频里所有的东西,用各种办法,加上自己的想法来处理,吃了很多苦。但,我觉得这些苦没有白吃,也收获了很多,等我他的第一部分彻底弄完(已经进行到P11),我会做个jimmy哥的整理版给大家,是可以真正重复他的视频的,希望对大家有帮助。
因为我觉得git是可以觉得,但是纯正的Linux不论是未来你做生信,还是想自己做点什么都是离不开的好工具(这是我吃了很多苦,同时也感觉到很多快乐),今天先给大家提供几个数据的下载,大家最好养成建立目录的习惯,不要乱存数据,否则后面很难找,也很难重复别人的视频,比如参考序列尽量放在一个reference的目录下,然后再逐级建立目录。

下载gh38

cd ~/reference
mkdir -p genome/hg38  && cd genome/hg38
nohup wget http://hgdownload.cse.ucsc.edu/goldenPath/hg38/bigZips/hg38.fa.gz  & 
gunzip hg38.fa.gz

命令说明

cd ~/reference #先进入reference目录,当然你没有的话cd reference建立一个就可以了

mkdir -p genome/hg38 && cd genome/hg38 #mkdir 是建立目录的意思,-P 是逐级建立,也就是再genome下再建立一个hg38;&& 是运行完前面的命令再运行后面的

nohup wget http://hgdownload.cse.ucsc.edu/goldenPath/hg38/bigZips/hg38.fa.gz & # nohup是你关闭这个输入终端,仍然运行的命令;wget是ubuntu的一个下载命令;网址是下载地址;&这个是不显示前端

gunzip hg38.fa.gz #这是解压缩.gz的命令

下载 hg19

cd ~/reference
mkdir -p genome/hg19  && cd genome/hg19
nohup wget ttp://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/chromFa.tar.gz &
tar zvfx chromFa.tar.gz
cat *.fa > hg19.fa
rm chr*.fa

命令说明

这个和上面的区别是压缩文件不一样,这个是tar.gz,所以需要解压gz和同时处理tar
命令就是tar zvfx;解压缩完是好多个chr*.fa;cat有很多用处,在这里是讲多个文件合并为一个;rm是删除之前解压出来的fa

下载小鼠

cd ~/reference

mkdir -p  genome/mm10  && cd genome/mm10

nohup wget http://hgdownload.cse.ucsc.edu/goldenPath/mm10/bigZips/chromFa.tar.gz  &

tar zvfx chromFa.tar.gz

cat *.fa > mm10.fa

rm chr*.fa

命令说明

没啥说明的,和hg19一个意思

下载SARS-CoV-2基因组

多说一句,SARS-CoV-2这个是病毒名(这个是国际病毒分类委员会干的)COVID-19这个是疾病名(这个是WHO干的)不要搞混。当然命名是否合理是另一回事,确实容易有奇异,但也是有一定道理的,至少没像MERS(ME代表Middle East),来个WHRS这样就已经很庆幸了,这种利用最早发病地区命名的方式,是应该舍弃了,这种影响是无法估量的,但在数据库还是叫做wuhCor1,这个是2月7日上传的,而WHO命名是2月11日命名的,不能怪他。所以先按照他的方式下载下来,我们再自己改名字。

image.png

image.png

image.png

这个是我自己根据上面的方式(jimmy哥的帖子),自己尝试的也欢迎大家指正,反正是下来了,也做了bwa的index,后面的再说。

cd ~/reference
mkdir SARS-CoV-2 && cd SARS-CoV-2
nohup wget https://hgdownload.soe.ucsc.edu/goldenPath/wuhCor1/bigZips/chromFa.tar.gz &
tar zvfx chromFa.tar.gz

今天先到这里

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 219,539评论 6 508
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,594评论 3 396
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 165,871评论 0 356
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,963评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,984评论 6 393
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,763评论 1 307
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,468评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,357评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,850评论 1 317
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,002评论 3 338
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,144评论 1 351
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,823评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,483评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,026评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,150评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,415评论 3 373
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,092评论 2 355

推荐阅读更多精彩内容