weka2019

摘要： 刚刚从发ncbi的后台下载了一套hg19的注释文件gff格式的，和ucsc上下载的refseqgene注释进行了对比，发现一个奇怪的问题: 同一个转录本的起始坐标相差10000. 花了一上午的时间对比了下才发现，原来ucsc上的基因组序列都在序列上前后各添加了10000个N. 然后又检查了genec 阅读全文

摘要： fastaread函数是matlab生物信息学工具箱内置的一个函数，给我们的使用上带来了巨大的方便。对于基因DNA序列，转录RNA序列和表达蛋白序列的读取非常方便。 使用语法为： 阅读全文

摘要： 进行没有终点，在基因位上放上尽可能多的功能，进化就是尽可能的扩展这种基因位。 matlab的很多命令其实 unix。例如mkdir, 这个函数用来创建文件夹，简直和unix终端命令一样。 类似的命令还有 cd 命令 ls 命令 dir 命令 。。。。 终于明白为什么很多老师打开matlab只保留命令 阅读全文

摘要： % 每次下载的基因组序列都很心累, 搞个自动化的clear all;clc;close all; url_base = 'http://hgdownload.soe.ucsc.edu/goldenPath/mm9/chromosomes/'; chr_list = {};for i = 1:19 c 阅读全文