获取基因组fasta中N的位置信息

yukaiquan
yukaiquan 00后

1 人点赞了该文章 · 383 浏览

背景

在基因组fasta文件中N的分布信息

示例

test.fasta:

>chr1
ATNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNtaaattgttt
taaattgtttctgtttgcagttgacatgatctNNNNNatagaaaacacca
ataactctgccaaaaaatttagaattcataaatgaatttagtaaagttgc
>chr2
NNNNNNNNNNNNNNNTTNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNtaaattgttt
taaattgtttctgtttgcagttgacatgatcttatatatagaaaacacca
ataactctgccaaaaaatttagaattcataaatgaatttagtaaagttgc

运行代码

chr_findn.exe/chr_findn 程序

test.fasta 输入的fasta

test.bed 输出的bed文件

win:

./chr_findn.exe test.fasta test.bed

linux:

./chr_findn test.fasta test.bed

 

test.bed

chr1	3	140
chr1	183	187
chr2	1	15
chr2	18	140

程序获取:fasta_find_n

发布于 2023-05-17 15:55

免责声明:

本文由 yukaiquan 原创发布于 生信坑 ,著作权归作者所有。

登录一下,更多精彩内容等你发现,贡献精彩回答,参与评论互动

登录! 还没有账号?去注册

暂无评论

文章目录

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024