trim_galore雙端數(shù)據(jù)過(guò)濾nohup trim_galore -q 20 --phred33 --stringency 3 --length 20 -e 0.1 --...
寫在前面 年初在汕頭婚宴上,我問(wèn)一個(gè)朋友,現(xiàn)在在做什么工作。他回復(fù)了我: 但行好事,莫問(wèn)前程。 現(xiàn)在看來(lái),挺好。 構(gòu)建進(jìn)化樹的基礎(chǔ)是序列對(duì)齊,或者說(shuō)多序列比對(duì)。做序列對(duì)齊的主...
將基因fa或蛋白fa序列文件由多行轉(zhuǎn)換為一行awk 'NR==1{a[NR]=$0}NR>1{if(a[NR-1] ~ ">"){print a[NR-1];a[NR]=$0...
示例fa文件: 第一種方法: 新建 script.sh,輸入以下內(nèi)容: 第二種方法: 只需一行awk命令(閑來(lái)無(wú)事突然想出來(lái)的騷操作): 之后運(yùn)行,第一個(gè)參數(shù)為需要轉(zhuǎn)換的多行...