bowtie2使用記錄

wget http://hgdownload.cse.ucsc.edu/goldenPath/mm10/bigZips/chromFa.tar.gz    下載參考基因組
tar -zxvf chromFa.tar.gz   解壓縮
cat *.fa > mm10.fa             合并文件

解壓縮后文件情況：

-rw-rw-r--. 1 yangjy yangjy 133308900 Feb 10  2012 chr10.fa
-rw-rw-r--. 1 yangjy yangjy 124524201 Feb 10  2012 chr11.fa
-rw-rw-r--. 1 yangjy yangjy 122531610 Feb 10  2012 chr12.fa
-rw-rw-r--. 1 yangjy yangjy 122830079 Feb 10  2012 chr13.fa
-rw-rw-r--. 1 yangjy yangjy 127400296 Feb 10  2012 chr14.fa
-rw-rw-r--. 1 yangjy yangjy 106124566 Feb 10  2012 chr15.fa
-rw-rw-r--. 1 yangjy yangjy 100171931 Feb 10  2012 chr16.fa
-rw-rw-r--. 1 yangjy yangjy  96887024 Feb 10  2012 chr17.fa
-rw-rw-r--. 1 yangjy yangjy  92516699 Feb 10  2012 chr18.fa
-rw-rw-r--. 1 yangjy yangjy  62660205 Feb 10  2012 chr19.fa
-rw-rw-r--. 1 yangjy yangjy 199381417 Feb 10  2012 chr1.fa
-rw-rw-r--. 1 yangjy yangjy    173142 Feb 10  2012 chr1_GL456210_random.fa
-rw-rw-r--. 1 yangjy yangjy    246592 Feb 10  2012 chr1_GL456211_random.fa
-rw-rw-r--. 1 yangjy yangjy    156713 Feb 10  2012 chr1_GL456212_random.fa
-rw-rw-r--. 1 yangjy yangjy     40149 Feb 10  2012 chr1_GL456213_random.fa
-rw-rw-r--. 1 yangjy yangjy    211123 Feb 10  2012 chr1_GL456221_random.fa
-rw-rw-r--. 1 yangjy yangjy 185755495 Feb 10  2012 chr2.fa
-rw-rw-r--. 1 yangjy yangjy 163240480 Feb 10  2012 chr3.fa
-rw-rw-r--. 1 yangjy yangjy 159638285 Feb 10  2012 chr4.fa
-rw-rw-r--. 1 yangjy yangjy     68029 Feb 10  2012 chr4_GL456216_random.fa
-rw-rw-r--. 1 yangjy yangjy    232548 Feb 10  2012 chr4_GL456350_random.fa
-rw-rw-r--. 1 yangjy yangjy     15266 Feb 10  2012 chr4_JH584292_random.fa
-rw-rw-r--. 1 yangjy yangjy    212150 Feb 10  2012 chr4_JH584293_random.fa
-rw-rw-r--. 1 yangjy yangjy    195766 Feb 10  2012 chr4_JH584294_random.fa
-rw-rw-r--. 1 yangjy yangjy      2038 Feb 10  2012 chr4_JH584295_random.fa
-rw-rw-r--. 1 yangjy yangjy 154871384 Feb 10  2012 chr5.fa
-rw-rw-r--. 1 yangjy yangjy    199935 Feb 10  2012 chr5_GL456354_random.fa
-rw-rw-r--. 1 yangjy yangjy    203378 Feb 10  2012 chr5_JH584296_random.fa
-rw-rw-r--. 1 yangjy yangjy    209914 Feb 10  2012 chr5_JH584297_random.fa
-rw-rw-r--. 1 yangjy yangjy    187895 Feb 10  2012 chr5_JH584298_random.fa
-rw-rw-r--. 1 yangjy yangjy    972095 Feb 10  2012 chr5_JH584299_random.fa
-rw-rw-r--. 1 yangjy yangjy 152731283 Feb 10  2012 chr6.fa
-rw-rw-r--. 1 yangjy yangjy 148350295 Feb 10  2012 chr7.fa
-rw-rw-r--. 1 yangjy yangjy    179510 Feb 10  2012 chr7_GL456219_random.fa
-rw-rw-r--. 1 yangjy yangjy 131989244 Feb 10  2012 chr8.fa
-rw-rw-r--. 1 yangjy yangjy 127087019 Feb 10  2012 chr9.fa
-rw-rw-r--. 1 yangjy yangjy     16631 Feb 10  2012 chrM.fa
-rw-rw-r--. 1 yangjy yangjy     40874 Feb 10  2012 chrUn_GL456239.fa
-rw-rw-r--. 1 yangjy yangjy     23450 Feb 10  2012 chrUn_GL456359.fa
-rw-rw-r--. 1 yangjy yangjy     32355 Feb 10  2012 chrUn_GL456360.fa
-rw-rw-r--. 1 yangjy yangjy     48031 Feb 10  2012 chrUn_GL456366.fa
-rw-rw-r--. 1 yangjy yangjy     42915 Feb 10  2012 chrUn_GL456367.fa
-rw-rw-r--. 1 yangjy yangjy     20629 Feb 10  2012 chrUn_GL456368.fa
-rw-rw-r--. 1 yangjy yangjy     27316 Feb 10  2012 chrUn_GL456370.fa
-rw-rw-r--. 1 yangjy yangjy     29254 Feb 10  2012 chrUn_GL456372.fa
-rw-rw-r--. 1 yangjy yangjy     32251 Feb 10  2012 chrUn_GL456378.fa
-rw-rw-r--. 1 yangjy yangjy     73849 Feb 10  2012 chrUn_GL456379.fa
-rw-rw-r--. 1 yangjy yangjy     26405 Feb 10  2012 chrUn_GL456381.fa
-rw-rw-r--. 1 yangjy yangjy     23638 Feb 10  2012 chrUn_GL456382.fa
-rw-rw-r--. 1 yangjy yangjy     39449 Feb 10  2012 chrUn_GL456383.fa
-rw-rw-r--. 1 yangjy yangjy     35961 Feb 10  2012 chrUn_GL456385.fa
-rw-rw-r--. 1 yangjy yangjy     25195 Feb 10  2012 chrUn_GL456387.fa
-rw-rw-r--. 1 yangjy yangjy     29364 Feb 10  2012 chrUn_GL456389.fa
-rw-rw-r--. 1 yangjy yangjy     25178 Feb 10  2012 chrUn_GL456390.fa
-rw-rw-r--. 1 yangjy yangjy     24118 Feb 10  2012 chrUn_GL456392.fa
-rw-rw-r--. 1 yangjy yangjy     56842 Feb 10  2012 chrUn_GL456393.fa
-rw-rw-r--. 1 yangjy yangjy     24826 Feb 10  2012 chrUn_GL456394.fa
-rw-rw-r--. 1 yangjy yangjy     21681 Feb 10  2012 chrUn_GL456396.fa
-rw-rw-r--. 1 yangjy yangjy    116758 Feb 10  2012 chrUn_JH584304.fa
-rw-rw-r--. 1 yangjy yangjy 174451931 Feb 10  2012 chrX.fa
-rw-rw-r--. 1 yangjy yangjy    343694 Feb 10  2012 chrX_GL456233_random.fa
-rw-rw-r--. 1 yangjy yangjy  93579598 Feb 10  2012 chrY.fa
-rw-rw-r--. 1 yangjy yangjy    186016 Feb 10  2012 chrY_JH584300_random.fa
-rw-rw-r--. 1 yangjy yangjy    265095 Feb 10  2012 chrY_JH584301_random.fa
-rw-rw-r--. 1 yangjy yangjy    158977 Feb 10  2012 chrY_JH584302_random.fa
-rw-rw-r--. 1 yangjy yangjy    161283 Feb 10  2012 chrY_JH584303_random.fa

合并文件
cat *.fa > mm10.fa
使用bowtie2 built建立索引：
bowtie2-build mm10.fa mm10
索引文件格式如下：

-rw-r--r--. 1 yangjy yangjy  888464705 May  2  2012 mm10.1.bt2
-rw-r--r--. 1 yangjy yangjy  663195880 May  2  2012 mm10.2.bt2
-rw-r--r--. 1 yangjy yangjy       6119 May  2  2012 mm10.3.bt2
-rw-r--r--. 1 yangjy yangjy  663195875 May  2  2012 mm10.4.bt2
-rw-r--r--. 1 yangjy yangjy  888464705 May  3  2012 mm10.rev.1.bt2
-rw-r--r--. 1 yangjy yangjy  663195880 May  3  2012 mm10.rev.2.bt2

處理信息（需要一定時間）

(chipseq) [yangjy@GSCG01 align]$ bowtie2-build mm10.fa mm10
Settings:
  Output files: "mm10.*.bt2"  # 輸出的文件格式：mm10.*.bt2
  Line rate: 6 (line is 64 bytes)
  Lines per side: 1 (side is 64 bytes)
  Offset rate: 4 (one in 16)
  FTable chars: 10
  Strings: unpacked
  Max bucket size: default
  Max bucket size, sqrt multiplier: default
  Max bucket size, len divisor: 4
  Difference-cover sample period: 1024
  Endianness: little
  Actual local endianness: little
  Sanity checking: disabled
  Assertions: disabled
  Random seed: 0
  Sizeofs: void*:8, int:4, long:8, size_t:8
Input files DNA, FASTA:
  mm10.fa
Building a SMALL index
Reading reference sizes
  Time reading reference sizes: 00:00:11
Calculating joined length
Writing header
Reserving space for joined string
Joining reference sequences
  Time to join reference sequences: 00:00:09
bmax according to bmaxDivN setting: 663195875
Using parameters --bmax 497396907 --dcv 1024
  Doing ahead-of-time memory usage test
  Passed!  Constructing with these parameters: --bmax 497396907 --dcv 1024
Constructing suffix-array element generator
Building DifferenceCoverSample
  Building sPrime
  Building sPrimeOrder
  V-Sorting samples
  V-Sorting samples time: 00:00:56
  Allocating rank array
  Ranking v-sort output
  Ranking v-sort output time: 00:00:14
  Invoking Larsson-Sadakane on ranks
  Invoking Larsson-Sadakane on ranks time: 00:00:25
  Sanity-checking and returning
Building samples
Reserving space for 12 sample suffixes
Generating random suffixes
QSorting 12 sample offsets, eliminating duplicates
QSorting sample offsets, eliminating duplicates time: 00:00:00
Multikey QSorting 12 samples
  (Using difference cover)
  Multikey QSorting samples time: 00:00:00
Calculating bucket sizes
Splitting and merging
  Splitting and merging time: 00:00:00
Avg bucket size: 2.65278e+09 (target: 497396906)
Converting suffix-array elements to index image
Allocating ftab, absorbFtab
Entering Ebwt loop
Getting block 1 of 1
  No samples; assembling all-inclusive block
  Sorting block of length 2652783500 for bucket 1
  (Using difference cover)

運行bowtie2 獲取 SAM 文件
bowtie2 -p 6 -3 5 --local -x mm10 -1 example_1.fastq -2 example_2.fastq -S example.sam

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 人面猴
    序言：七十年代末撞芍，一起剝皮案震驚了整個濱河市，隨后出現(xiàn)的幾起案子跨扮，更是在濱河造成了極大的恐慌序无，老刑警劉巖，帶你破解...
    沈念sama閱讀 206,126評論 6贊 481
  • 死咒
    序言：濱河連續(xù)發(fā)生了三起死亡事件衡创，死亡現(xiàn)場離奇詭異帝嗡，居然都是意外死亡，警方通過查閱死者的電腦和手機璃氢，發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,254評論 2贊 382
  • 救了他兩次的神仙讓他今天三更去死
    文/潘曉璐 我一進店門哟玷，熙熙樓的掌柜王于貴愁眉苦臉地迎上來，“玉大人一也，你說我怎么就攤上這事巢寡。” “怎么了椰苟？”我有些...
    開封第一講書人閱讀 152,445評論 0贊 341
  • 道士緝兇錄：失蹤的賣姜人
    文/不壞的土叔 我叫張陵抑月，是天一觀的道長。 經(jīng)常有香客問我舆蝴，道長谦絮，這世上最難降的妖魔是什么？ 我笑而不...
    開封第一講書人閱讀 55,185評論 1贊 278
  • ?港島之戀（遺憾婚禮）
    正文 為了忘掉前任洁仗，我火速辦了婚禮层皱，結(jié)果婚禮上，老公的妹妹穿的比我還像新娘赠潦。我一直安慰自己叫胖，他們只是感情好，可當(dāng)我...
    茶點故事閱讀 64,178評論 5贊 371
  • 惡毒庶女頂嫁案：這布局不是一般人想出來的
    文/花漫 我一把揭開白布祭椰。 她就那樣靜靜地躺著臭家，像睡著了一般。 火紅的嫁衣襯著肌膚如雪方淤。 梳的紋絲不亂的頭發(fā)上钉赁，一...
    開封第一講書人閱讀 48,970評論 1贊 284
  • 城市分裂傳說
    那天，我揣著相機與錄音携茂，去河邊找鬼你踩。 笑死，一個胖子當(dāng)著我的面吹牛，可吹牛的內(nèi)容都是我干的带膜。 我是一名探鬼主播吩谦，決...
    沈念sama閱讀 38,276評論 3贊 399
  • 雙鴛鴦連環(huán)套：你想象不到人心有多黑
    文/蒼蘭香墨 我猛地睜開眼，長吁一口氣：“原來是場噩夢啊……” “哼膝藕！你這毒婦竟也來了式廷？” 一聲冷哼從身側(cè)響起，我...
    開封第一講書人閱讀 36,927評論 0贊 259
  • 萬榮殺人案實錄
    序言：老撾萬榮一對情侶失蹤芭挽，失蹤者是張志新（化名）和其女友劉穎滑废，沒想到半個月后，有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體袜爪，經(jīng)...
    沈念sama閱讀 43,400評論 1贊 300
  • ?護林員之死
    正文 獨居荒郊野嶺守林人離奇死亡蠕趁，尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,883評論 2贊 323
  • ?白月光啟示錄
    正文 我和宋清朗相戀三年，在試婚紗的時候發(fā)現(xiàn)自己被綠了辛馆。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片俺陋。...
    茶點故事閱讀 37,997評論 1贊 333
  • 活死人
    序言：一個原本活蹦亂跳的男人離奇死亡，死狀恐怖昙篙，靈堂內(nèi)的尸體忽然破棺而出腊状，到底是詐尸還是另有隱情，我是刑警寧澤苔可，帶...
    沈念sama閱讀 33,646評論 4贊 322
  • ?日本核電站爆炸內(nèi)幕
    正文 年R本政府宣布寿酌，位于F島的核電站，受9級特大地震影響硕蛹，放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜硕并，卻給世界環(huán)境...
    茶點故事閱讀 39,213評論 3贊 307
  • 男人毒藥：我在死后第九天來索命
    文/蒙蒙 一法焰、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧倔毙，春花似錦埃仪、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,204評論 0贊 19
  • 一樁弒父案卵蛉，背后竟有這般陰謀
    文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至么库，卻和暖如春傻丝，著一層夾襖步出監(jiān)牢的瞬間，已是汗流浹背诉儒。 一陣腳步聲響...
    開封第一講書人閱讀 31,423評論 1贊 260
  • 情欲美人皮
    我被黑心中介騙來泰國打工葡缰， 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留，地道東北人。 一個月前我還...
    沈念sama閱讀 45,423評論 2贊 352
  • 代替公主和親
    正文 我出身青樓泛释，卻偏偏與公主長得像滤愕，于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子怜校，可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 42,722評論 2贊 345

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