計算線粒體基因組每個堿基位點的Depth

# 建立線粒體基因組序列索引

bwa index mitochondria.fasta

#mitochondria.fasta替換成已經(jīng)組裝好的線粒體基因組序列文件

# 比對

bwa mem mitochondria.fasta read1.fastq read2.fastq > alignment.sam

#sam轉(zhuǎn)bam格式

samtools view -Sb alignment.sam > alignment.bam

#bam文件排序

samtools sort alignment.bam -o sorted_alignment.bam

#samtools計算depth

samtools depth sorted_alignment.bam > depth.txt

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末跃洛,一起剝皮案震驚了整個濱河市惜颇,隨后出現(xiàn)的幾起案子径荔,更是在濱河造成了極大的恐慌,老刑警劉巖讲婚,帶你破解...
    沈念sama閱讀 211,639評論 6 492
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異衍腥,居然都是意外死亡磺樱,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,277評論 3 385
  • 文/潘曉璐 我一進店門婆咸,熙熙樓的掌柜王于貴愁眉苦臉地迎上來竹捉,“玉大人,你說我怎么就攤上這事尚骄】椴睿” “怎么了?”我有些...
    開封第一講書人閱讀 157,221評論 0 348
  • 文/不壞的土叔 我叫張陵倔丈,是天一觀的道長憨闰。 經(jīng)常有香客問我,道長需五,這世上最難降的妖魔是什么鹉动? 我笑而不...
    開封第一講書人閱讀 56,474評論 1 283
  • 正文 為了忘掉前任,我火速辦了婚禮宏邮,結(jié)果婚禮上泽示,老公的妹妹穿的比我還像新娘缸血。我一直安慰自己,他們只是感情好械筛,可當我...
    茶點故事閱讀 65,570評論 6 386
  • 文/花漫 我一把揭開白布捎泻。 她就那樣靜靜地躺著,像睡著了一般埋哟。 火紅的嫁衣襯著肌膚如雪笆豁。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,816評論 1 290
  • 那天赤赊,我揣著相機與錄音闯狱,去河邊找鬼。 笑死抛计,一個胖子當著我的面吹牛扩氢,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播爷辱,決...
    沈念sama閱讀 38,957評論 3 408
  • 文/蒼蘭香墨 我猛地睜開眼录豺,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了饭弓?” 一聲冷哼從身側(cè)響起双饥,我...
    開封第一講書人閱讀 37,718評論 0 266
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎弟断,沒想到半個月后咏花,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 44,176評論 1 303
  • 正文 獨居荒郊野嶺守林人離奇死亡阀趴,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,511評論 2 327
  • 正文 我和宋清朗相戀三年昏翰,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片刘急。...
    茶點故事閱讀 38,646評論 1 340
  • 序言:一個原本活蹦亂跳的男人離奇死亡棚菊,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出叔汁,到底是詐尸還是另有隱情统求,我是刑警寧澤,帶...
    沈念sama閱讀 34,322評論 4 330
  • 正文 年R本政府宣布据块,位于F島的核電站码邻,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏另假。R本人自食惡果不足惜像屋,卻給世界環(huán)境...
    茶點故事閱讀 39,934評論 3 313
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望边篮。 院中可真熱鬧己莺,春花似錦因苹、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,755評論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽凶杖。三九已至胁艰,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間智蝠,已是汗流浹背腾么。 一陣腳步聲響...
    開封第一講書人閱讀 31,987評論 1 266
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留杈湾,地道東北人解虱。 一個月前我還...
    沈念sama閱讀 46,358評論 2 360
  • 正文 我出身青樓,卻偏偏與公主長得像漆撞,于是被迫代替她去往敵國和親殴泰。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 43,514評論 2 348

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