我們可以先用公司提供的bam檔案提取出mapq為0的reads[裡面包含大多數 multiple reads],看一下大致的分布,後面再用我自己比對得到的bam作進更加準確的分析。
mateinfo:
mate 序列所在參考序列的名稱;mate 序列在參考序列上的位置;估計出的片段的長度,當mate 序列位於本序列上游時該值為負值。
samtools index in.bam # 生成in.bam的索引檔案in.bam.bai
samtools view in.bam chr22 # 跳轉到chr22染色體
samtools view in.bam chr22:16050103 # 跳轉到chr22:16050103位置
samtools view in.bam chr22:16050103-16050103 # 只檢視該位置
samtools view -h in.bam chr22:16050103-16050203 | samtools view -sb - > small.bam
samtools tview --reference hg38.fa in.bam #該模式下,按下鍵盤『g』後,會跳出乙個goto框
生物資訊學軟體 自學生物資訊學
我是生物工程專業出身,在大三保研時選擇了生物資訊的道路,到現在為止已經在行業裡摸爬滾打了6年的時間,在這6年的學習之路上疑惑過,也迷茫過,特此把我學習的過程以及遇到的問題總結出來以讓大家避免出現同樣的問題。在我學習生物資訊過程的基礎上帶著大家順暢的走一遍。在學習生物資訊學之前,我們先來了解一下什麼是...
生物資訊學(Bioinformatics)
生物資訊學 bioinformatics 這個名詞有許多不同的定義。從字面上來看,生物資訊學是將資訊科學應用於生物學。生物資訊學廣義的概念是指應用資訊科學研究生物體系和生物過程中資訊的存貯 資訊的內涵和資訊的傳遞,研究和分析生物體細胞 組織 器官的生理 病理 藥理過程的中各種生物資訊,或者說是生命科...
生物資訊學緒論
hgp第乙個五年總結報告 生物資訊學是一門交叉學科 它包含了生物資訊的獲取 加工 儲存 分配 分析 解釋等在內的所有方面,它運用數學 電腦科學和生物學的各種工具闡明和解釋大量資料所包含的生物學意義。美國喬治亞理工大學 生物資訊學是採用數學 統計學和電腦科學,分析生物學 生物化學和生物物理學資料的一門...