trim_galore双端数据过滤nohup trim_galore -q 20 --phred33 --stringency 3 --length 20 -e 0.1 --...
写在前面 年初在汕头婚宴上,我问一个朋友,现在在做什么工作。他回复了我: 但行好事,莫问前程。 现在看来,挺好。 构建进化树的基础是序列对齐,或者说多序列比对。做序列对齐的主...
将基因fa或蛋白fa序列文件由多行转换为一行awk 'NR==1{a[NR]=$0}NR>1{if(a[NR-1] ~ ">"){print a[NR-1];a[NR]=$0...
示例fa文件: 第一种方法: 新建 script.sh,输入以下内容: 第二种方法: 只需一行awk命令(闲来无事突然想出来的骚操作): 之后运行,第一个参数为需要转换的多行...