SAM/BAM文件處理
當測序得到的fastq文件map到基因組之后，我們通常會得到一個sam或者bam為擴展名的文件。SAM的全稱是sequence alignment/map format。而BAM就是SAM的二進制文件(B取自binary)。 那么SAM文件的格式是什么樣子的呢？如果你想真實地了解SAM文件 ...
原文：Pysam 處理bam文件
Pysam可用來處理bam文件 安裝： 用 pip 或者 conda即可 使用： Pysam的函數有很多，主要的讀取函數有： AlignmentFile：讀取BAM CRAM SAM文件 VariantFile：讀取變異數據 VCF或者BCF TabixFile：讀取由tabix索引的文件 FastaFile：讀取fasta序列文件 FastqFile：讀取fastq測序序列文件 一般常用的是第一 ...
2019-12-05 17:03 0 469 推薦指數：
相關推薦
pysam - 多種格式基因組數據（sam/bam/vcf/bcf/cram/…）讀寫與處理模塊（python）
在開發基因組相關流程或工具時，經常需要讀取、處理和創建bam、vcf、bcf文件。目前已經有一些主流的處理此類格式文件的工具，如samtools、picard、vcftools、bcftools，但此類工具集成的大多是標准功能，在編程時如果直接調用的話往往顯得不夠靈活。 本文介紹的是一個處理 ...
SAMTOOLS使用 SAM BAM文件處理
【怪毛匠子 整理】 samtools學習及使用范例,以及官方文檔詳解 #第一步：把sam文件轉換成bam文件,我們得到map.bam文件 system"samtools view -bS map.sam > map.bam"; #第二步：sort 一下 BAM ...
處理bam文件提取信息
一般來說,一個bam文件通常只包含一個樣本的信息,最多需要進行染色體位置的處理, samtools也提供了簡單的處理方式,比如要提取 chr1的reads, 只需要: samtools view input.bam ch1 這幾天遇到了10x genomics的bam結果,發現 ...
pysam操作sam文件
pysam模塊 因為要分析sam文件中序列的情況，因此要對reads進行細分，所以之前想用數據庫將sam文件信息存儲，然后用sql語句進行分類。后來發現很麻煩，pysam就是一個高效讀取存儲在SAM / BAM / CRAM格式文件中的映射短讀序列數據信息的python模塊，可以輕松 ...
bam文件格式說明
bam文件說明 bam文件和sam文件內容其實是一樣的，只是bam是二進制的壓縮文件，需要通過特定的軟件來進行查看，bam文件通常可以理解為12個字段組成 BAM格式分為header section（頭部分，注釋信息，以@開頭，可有可無）和alignment section（比對 ...
使用Python處理BAM的方法
在上一篇的文章里我詳細介紹了BAM（SAM/CRAM）的格式和一些需要注意的細節，還說了該如何使用samtools在命令行中對其進行操作。但是很多時候這些操作是不能滿足我們的實際需要的，比如統計比對率、計算在某個比對質量值之上的read有多少，或者計算PE比對的插入片段長度分布，甚至需要 ...
bam/sam 文件格式詳解
sam/bam 是一種序列比對格式標准，由sanger制定，是以TAB為分割符的文本格式。主要應用於測序序列mapping到基因組上的結果表示，當然也可以表示任意的多重比對結果。通常是把FASTQ文件格式的測序數據比對到對應的參考基因組版本得到的。 header 部分 sam 分為兩部分，注釋 ...