seqkit的安裝與使用

安裝:

conda install seqkit

ref:seqkit一個(gè)FASTA/Q序列處理神器 - 遺世獨(dú)立的愚公 - 博客園

一、序列操作:

1.取反向序列

seqkit ? seq ?test.fa ? -r ?> ?test_re.fa

2.取互補(bǔ)序列

seqkit ? seq ? test.fa ?-p ?> ?test_com.fa

3.取反向互補(bǔ)序列

seqkit ? seq ? test.fa ?-r ?-p ?> test_re_com.fa

4.DNA序列轉(zhuǎn)換為RNA序列

seqkit ? seq ? test.fa ?--nda2rna ? > ? test_rna.fa

5.RNA序列轉(zhuǎn)換為DNA序列

seqkit ? seq ?test.fa ? rna2dna ? ? > ? ?test_dna.fa

6.將序列以小寫字母的形式輸出

seqkit ?seq ?test.fa ?-l ?> ?test_lower.fa

7.將序列以大寫字母的形式輸出

seqkit ? seq ? test.fa ?-u > ?test_upper.fa

8.指定每行序列的輸出長(zhǎng)度(為0的話瞬逊,代表為一整行显歧,默認(rèn)的輸出 長(zhǎng)度是60個(gè)堿基)

seqkit ?seq ?test.fa ?-w ?10 ?> ?test_10.fa ?(指定序列的長(zhǎng)度為10)

9.將多行序列轉(zhuǎn)換為一行序列

seqkit ? seq ?test.fa ? -w ? 0 ? > ?test_w.fa

10.只輸出序列

seqkit ? seq ?test.fa ?-s ?-w 0 > test_seq.fa

11.將只輸出的序列的,指定每行輸出的堿基數(shù)

seqkit ? seq ?test_seq.fa ?-s ?-w 40 > test_seq40.fa

###注意10,11的微妙之處

###11,12也可以一步完成:

seqkit ?seq ? test.fa ? -s ?-w ?20 ?-o ?test_20.fa

二确镊、Fasta/q之間以及與tab格式互換

10.將fataq文件轉(zhuǎn)化為fasta格式.

seqkit fq2fa ? test.fq ? -o ? test.fa

11.將fasta格式轉(zhuǎn)化為tab格式

seqkit ?fx2tab ?test.fa > ?test_tab.fa (沒有seq參數(shù))

三士骤、序列信息統(tǒng)計(jì)

1.序列堿基含量

seqkit ?fx2tab ?-l ?-g ?-n ?-i ?-H ?test.fa

(這些參數(shù)組合起來比較好看)

2.序列長(zhǎng)度的整體分布統(tǒng)計(jì)

seqkit ?stat ?test.fa

四、其他用法:功能還是挺多的:

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末蕾域,一起剝皮案震驚了整個(gè)濱河市拷肌,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌旨巷,老刑警劉巖巨缘,帶你破解...
    沈念sama閱讀 217,826評(píng)論 6 506
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異采呐,居然都是意外死亡若锁,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,968評(píng)論 3 395
  • 文/潘曉璐 我一進(jìn)店門懈万,熙熙樓的掌柜王于貴愁眉苦臉地迎上來拴清,“玉大人,你說我怎么就攤上這事会通】谟瑁” “怎么了?”我有些...
    開封第一講書人閱讀 164,234評(píng)論 0 354
  • 文/不壞的土叔 我叫張陵涕侈,是天一觀的道長(zhǎng)沪停。 經(jīng)常有香客問我,道長(zhǎng)裳涛,這世上最難降的妖魔是什么木张? 我笑而不...
    開封第一講書人閱讀 58,562評(píng)論 1 293
  • 正文 為了忘掉前任,我火速辦了婚禮端三,結(jié)果婚禮上舷礼,老公的妹妹穿的比我還像新娘。我一直安慰自己郊闯,他們只是感情好妻献,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,611評(píng)論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著团赁,像睡著了一般育拨。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上欢摄,一...
    開封第一講書人閱讀 51,482評(píng)論 1 302
  • 那天熬丧,我揣著相機(jī)與錄音,去河邊找鬼怀挠。 笑死析蝴,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的绿淋。 我是一名探鬼主播嫌变,決...
    沈念sama閱讀 40,271評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼躬它!你這毒婦竟也來了腾啥?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,166評(píng)論 0 276
  • 序言:老撾萬榮一對(duì)情侶失蹤冯吓,失蹤者是張志新(化名)和其女友劉穎倘待,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體组贺,經(jīng)...
    沈念sama閱讀 45,608評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡凸舵,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,814評(píng)論 3 336
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了失尖。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片啊奄。...
    茶點(diǎn)故事閱讀 39,926評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡渐苏,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出菇夸,到底是詐尸還是另有隱情琼富,我是刑警寧澤,帶...
    沈念sama閱讀 35,644評(píng)論 5 346
  • 正文 年R本政府宣布庄新,位于F島的核電站鞠眉,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏择诈。R本人自食惡果不足惜械蹋,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,249評(píng)論 3 329
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望羞芍。 院中可真熱鬧哗戈,春花似錦、人聲如沸荷科。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,866評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽步做。三九已至副渴,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間全度,已是汗流浹背煮剧。 一陣腳步聲響...
    開封第一講書人閱讀 32,991評(píng)論 1 269
  • 我被黑心中介騙來泰國(guó)打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留将鸵,地道東北人勉盅。 一個(gè)月前我還...
    沈念sama閱讀 48,063評(píng)論 3 370
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像顶掉,于是被迫代替她去往敵國(guó)和親草娜。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,871評(píng)論 2 354

推薦閱讀更多精彩內(nèi)容