hisat2比對

基因組比對軟件常用bwa,轉(zhuǎn)錄組比對軟件常用bowtie2红碑、hisat2等舞吭,其中有參考基因組的常用hisat2泡垃,沒有參考基因組的常用bowtie2。

一羡鸥、建立索引

建立基因組索引

* hisat2-build –p 4 genome.fa genome

建立基因組+轉(zhuǎn)錄組+SNP索引:
bowtie2的索引只有基因組序列信息蔑穴,tophat2比對時,轉(zhuǎn)錄組信息通過-G參數(shù)指定惧浴。HISAT2建立索引時存和,就應(yīng)該把轉(zhuǎn)錄組信息加進去。
HISAT2提供兩個Python腳本將GTF文件轉(zhuǎn)換成hisat2-build能使用的文件:

extract_exons.py Homo_sapiens.GRCh38.83.chr.gtf > genome.exon
extract_splice_sites.py Homo_sapiens.GRCh38.83.chr.gtf > genome.ss
此外衷旅,HISAT2還支持將SNP信息加入到索引中捐腿,這樣比對的時候就可以考慮SNP的情況。這仍然需要將SNP文件轉(zhuǎn)換成hisat2-build能使用的文件:

extract_snps.py snp142Common.txt > genome.snp
最后柿顶,將基因組茄袖、轉(zhuǎn)錄組、SNP建立索引:

* hisat2-build -p4 genome.fa --snp genome.snp --ss genome.ss --exon genome.exon genome_snp_tran

官網(wǎng)提供了人和小鼠的索引文件下載嘁锯,壓縮包有make_grch38_tran.sh文件绞佩,詳細記錄了創(chuàng)建索引的過程。

二猪钮、運行HISAT2

* hisat2 -p10 -x ./genome -1 Sample.R1.fastq -2 Sample.R2.fastq --rna-strandness RF --fr–S Sample.sam

-p 線程數(shù)
--rna-strandness RF 鏈特異性
-x 指定基因組索引

-1 指定第一個fastq文件
-2 指定第二個fastq文件
-S 指定輸出的SAM文件

最終品山,我們需要使用samtools軟件,對sam排序得到一個sorted.bam文件烤低,用于后面的定量肘交,AS等分析;

* samtools view -uS Sample.hisat2.sam |samtools sort - -o Sample.sorted.bam && samtools index NC5.sorted.bam

官方操作手冊簡要版

用法:
hisat2 [options]* -x <hisat2-idx> {-1 <m1> -2 <m2> | -U <r> |[-S <hit>]

主要參數(shù):
-x <hisat2-idx>
參考基因組索引文件的前綴扑馁。
-1 <m1>
雙端測序結(jié)果的第一個文件涯呻。若有多組數(shù)據(jù),使用逗號將文件分隔腻要。Reads的長度可以不一致复罐。
-2 <m2>
雙端測序結(jié)果的第二個文件。若有多組數(shù)據(jù)雄家,使用逗號將文件分隔效诅,并且文件順序要和-1參數(shù)對應(yīng)。Reads的長度可以不一致趟济。
-U <r>
單端數(shù)據(jù)文件乱投。若有多組數(shù)據(jù),使用逗號將文件分隔顷编∑蒽牛可以和-1、-2參數(shù)同時使用媳纬。Reads的長度可以不一致双肤。
–sra-acc <SRA accession number>
輸入SRA登錄號施掏,比如SRR353653,SRR353654茅糜。多組數(shù)據(jù)之間使用逗號分隔七芭。HISAT將自動下載并識別數(shù)據(jù)類型,進行比對限匣。
-S <hit>
指定輸出的SAM文件抖苦。

輸入選項:
-q:輸入文件為FASTQ格式毁菱。FASTQ格式為默認參數(shù)米死。
-qseq :輸入文件為QSEQ格式。
-f:輸入文件為FASTA格式贮庞。
-r:輸入文件中峦筒,每一行代表一條序列,沒有序列名和測序質(zhì)量等窗慎。選擇此項時物喷,–ignore-quals參數(shù)也會被選擇。
-c:此參數(shù)后是直接比對的序列遮斥,而不是包含序列的文件名峦失。序列間用逗號隔開。選擇此項時术吗,–ignore-quals參數(shù)也會被選擇尉辑。
-s/–skip <int>:跳過輸入文件中前條序列進行比對。
-u/–qupto <int>:只使用輸入文件中前條序列進行比對较屿,默認是沒有限制隧魄。
-5/–trim5 <int>:比對前去除每條序列5’端個堿基
-3/–trim3 <int>:比對前去除每條序列3’端個堿基
–phred33:輸入的FASTQ文件堿基質(zhì)量值編碼標準為phred33,phred33為默認參數(shù)隘蝎。
–phred64:輸入的FASTQ文件堿基質(zhì)量值編碼標準為phred64购啄。
–solexa-quals:將Solexa的堿基質(zhì)量值編碼標準轉(zhuǎn)換為phred。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末嘱么,一起剝皮案震驚了整個濱河市狮含,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌曼振,老刑警劉巖辉川,帶你破解...
    沈念sama閱讀 222,627評論 6 517
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件赋兵,死亡現(xiàn)場離奇詭異侣姆,居然都是意外死亡慌烧,警方通過查閱死者的電腦和手機令境,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 95,180評論 3 399
  • 文/潘曉璐 我一進店門钉赁,熙熙樓的掌柜王于貴愁眉苦臉地迎上來局嘁,“玉大人川无,你說我怎么就攤上這事喇喉。” “怎么了妆距?”我有些...
    開封第一講書人閱讀 169,346評論 0 362
  • 文/不壞的土叔 我叫張陵穷遂,是天一觀的道長。 經(jīng)常有香客問我娱据,道長蚪黑,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 60,097評論 1 300
  • 正文 為了忘掉前任中剩,我火速辦了婚禮忌穿,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘结啼。我一直安慰自己掠剑,他們只是感情好,可當我...
    茶點故事閱讀 69,100評論 6 398
  • 文/花漫 我一把揭開白布郊愧。 她就那樣靜靜地躺著朴译,像睡著了一般。 火紅的嫁衣襯著肌膚如雪属铁。 梳的紋絲不亂的頭發(fā)上眠寿,一...
    開封第一講書人閱讀 52,696評論 1 312
  • 那天,我揣著相機與錄音焦蘑,去河邊找鬼盯拱。 笑死,一個胖子當著我的面吹牛喇肋,可吹牛的內(nèi)容都是我干的坟乾。 我是一名探鬼主播,決...
    沈念sama閱讀 41,165評論 3 422
  • 文/蒼蘭香墨 我猛地睜開眼蝶防,長吁一口氣:“原來是場噩夢啊……” “哼甚侣!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起间学,我...
    開封第一講書人閱讀 40,108評論 0 277
  • 序言:老撾萬榮一對情侶失蹤殷费,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后低葫,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體详羡,經(jīng)...
    沈念sama閱讀 46,646評論 1 319
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 38,709評論 3 342
  • 正文 我和宋清朗相戀三年嘿悬,在試婚紗的時候發(fā)現(xiàn)自己被綠了实柠。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 40,861評論 1 353
  • 序言:一個原本活蹦亂跳的男人離奇死亡善涨,死狀恐怖窒盐,靈堂內(nèi)的尸體忽然破棺而出草则,到底是詐尸還是另有隱情,我是刑警寧澤蟹漓,帶...
    沈念sama閱讀 36,527評論 5 351
  • 正文 年R本政府宣布炕横,位于F島的核電站,受9級特大地震影響葡粒,放射性物質(zhì)發(fā)生泄漏份殿。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 42,196評論 3 336
  • 文/蒙蒙 一嗽交、第九天 我趴在偏房一處隱蔽的房頂上張望卿嘲。 院中可真熱鬧,春花似錦轮纫、人聲如沸腔寡。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,698評論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至忿磅,卻和暖如春糯彬,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背葱她。 一陣腳步聲響...
    開封第一講書人閱讀 33,804評論 1 274
  • 我被黑心中介騙來泰國打工撩扒, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人吨些。 一個月前我還...
    沈念sama閱讀 49,287評論 3 379
  • 正文 我出身青樓搓谆,卻偏偏與公主長得像,于是被迫代替她去往敵國和親豪墅。 傳聞我的和親對象是個殘疾皇子泉手,可洞房花燭夜當晚...
    茶點故事閱讀 45,860評論 2 361