worklog——pindel炬转,GATK,delly使用

pindel:

需要文件:sort過的bam文件算灸;.fa的參考文件扼劈;config.txt文件(主要是用于標(biāo)準(zhǔn)化bam文件);bam.bai文件菲驴;.fa.fai文件將這些文件均放在一個(gè)目錄A下

步驟:

①進(jìn)入文件A荐吵,在目錄A中輸入命令:

./pindel -f xxx.fa -i xxxx_config.txt ?-c [chromosome] ?-o [outputfilename]

②則在A目錄下會(huì)出現(xiàn)輸出文件

GATK:

工作原理:

(一)GATK三大步驟:原始數(shù)據(jù)的處理——變異檢測(cè)——初步分析

(二)前期準(zhǔn)備:①軟件:須安裝picard,bwa赊瞬,htslib

②需要文件:fastq文件

(三)原始數(shù)據(jù)處理:

①用bwa mem對(duì)fastq文件進(jìn)行過濾比對(duì)——得到sam文件

java -jar /home/lsw/xijiaoda/GATK/GenomeAnalysisTK.jar -T CountReads -R exampleFASTA.fasta -I exampleBAM.bam

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末先煎,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子巧涧,更是在濱河造成了極大的恐慌薯蝎,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,997評(píng)論 6 502
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件谤绳,死亡現(xiàn)場(chǎng)離奇詭異占锯,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)缩筛,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,603評(píng)論 3 392
  • 文/潘曉璐 我一進(jìn)店門消略,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人瞎抛,你說我怎么就攤上這事艺演。” “怎么了桐臊?”我有些...
    開封第一講書人閱讀 163,359評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵胎撤,是天一觀的道長(zhǎng)。 經(jīng)常有香客問我豪硅,道長(zhǎng)哩照,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,309評(píng)論 1 292
  • 正文 為了忘掉前任懒浮,我火速辦了婚禮飘弧,結(jié)果婚禮上识藤,老公的妹妹穿的比我還像新娘。我一直安慰自己次伶,他們只是感情好痴昧,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,346評(píng)論 6 390
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著冠王,像睡著了一般赶撰。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上柱彻,一...
    開封第一講書人閱讀 51,258評(píng)論 1 300
  • 那天豪娜,我揣著相機(jī)與錄音,去河邊找鬼哟楷。 笑死瘤载,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的卖擅。 我是一名探鬼主播鸣奔,決...
    沈念sama閱讀 40,122評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼惩阶!你這毒婦竟也來了挎狸?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 38,970評(píng)論 0 275
  • 序言:老撾萬榮一對(duì)情侶失蹤断楷,失蹤者是張志新(化名)和其女友劉穎锨匆,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體脐嫂,經(jīng)...
    沈念sama閱讀 45,403評(píng)論 1 313
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡统刮,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,596評(píng)論 3 334
  • 正文 我和宋清朗相戀三年紊遵,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了账千。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 39,769評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡暗膜,死狀恐怖匀奏,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情学搜,我是刑警寧澤娃善,帶...
    沈念sama閱讀 35,464評(píng)論 5 344
  • 正文 年R本政府宣布,位于F島的核電站瑞佩,受9級(jí)特大地震影響聚磺,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜炬丸,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,075評(píng)論 3 327
  • 文/蒙蒙 一瘫寝、第九天 我趴在偏房一處隱蔽的房頂上張望蜒蕾。 院中可真熱鬧,春花似錦焕阿、人聲如沸咪啡。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,705評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)撤摸。三九已至,卻和暖如春褒纲,著一層夾襖步出監(jiān)牢的瞬間准夷,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,848評(píng)論 1 269
  • 我被黑心中介騙來泰國(guó)打工莺掠, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留冕象,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,831評(píng)論 2 370
  • 正文 我出身青樓汁蝶,卻偏偏與公主長(zhǎng)得像渐扮,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子掖棉,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,678評(píng)論 2 354

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

  • Spring Cloud為開發(fā)人員提供了快速構(gòu)建分布式系統(tǒng)中一些常見模式的工具(例如配置管理墓律,服務(wù)發(fā)現(xiàn),斷路器幔亥,智...
    卡卡羅2017閱讀 134,654評(píng)論 18 139
  • 《ijs》速成開發(fā)手冊(cè)3.0 官方用戶交流:iApp開發(fā)交流(1) 239547050iApp開發(fā)交流(2) 10...
    葉染柒丶閱讀 5,127評(píng)論 0 7
  • 第十章 使用序列數(shù)據(jù) 生物信息學(xué)的核心問題之一是處理大量的(通常定義糟糕或模糊)文件格式耻讽。久而久之,一些特定的簡(jiǎn)單...
    yangliunk1987閱讀 5,025評(píng)論 3 53
  • linux資料總章2.1 1.0寫的不好抱歉 但是2.0已經(jīng)改了很多 但是錯(cuò)誤還是無法避免 以后資料會(huì)慢慢更新 大...
    數(shù)據(jù)革命閱讀 12,160評(píng)論 2 33
  • 報(bào)怨形體不美帕棉;從零開始针肥;沒有把握住自己;圍城香伴,生機(jī)勃勃的學(xué)生過幾年變成平庸的人慰枕,不斷的設(shè)計(jì)起點(diǎn),不斷人突破圍城即纲,...
    D051飛鷹閱讀 176評(píng)論 0 0