24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1609  |  回复: 1

Rhea_nyx

新虫 (初入文坛)

[交流] 用matlab从基因组(fasta格式)中提取特定序列已有1人参与

现在有一个fasta格式的基因组序列,怎样从中提取特定位置,比如5000-5500的序列?
【想问一下代码里可以直接打开fa文件就是写成  fid=fopen('*****.fa','r')  这样子吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainbow3311

铁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
data = fastaread('sequence.fasta');
da2=data.sequence;
first_500_bases = seqdisp(da2(5000:5500));
2楼2020-11-21 17:12:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Rhea_nyx 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见