使用i**檢視bam檔案的時候，發現i**有時候並不能顯示所有的reads（因為chr13:32907421，顯示有ins，但是沒有找到這個ins，懷疑i**沒有顯示所有reads）
所以想要檢視比對這個位點的所有reads，然後再看這個位點ins資訊。
方法一：
samtools view -h samplename_t.realn.bam "chr13:32907420-32907421" > chr13_32907421.sam方法二 :
samtools tview samplename_t.realn.bam -p "chr13:32907420-32907421" -d t > chr13_32907421.text將找到ins 的這些 reads name 去i**裡檢視，確實沒有相關reads，也是很奇怪了，以後找個時間研究一下 i**的設定再看看。
處理bam檔案提取資訊
一般來說,乙個bam檔案通常只包含乙個樣本的資訊,最多需要進行染色體位置的處理,samtools也提供了簡單的處理方式,比如要提取 chr1的reads,只需要 samtools view input.bam ch1 這幾天遇到了10x genomics的bam結果,發現單細胞的reads全包含在乙...
從檔案中提取數字
程式的功能是從指定的檔案中提取所有的整數，並把它們依次存入到乙個字串中 void iostest char a 50 char b 50 istrstream sin a ostrstream sout b,sizeof b ifstream file1 w1.dat ios in ios nocr...
MATLAB之從fig檔案中提取資料
open figname.fig 開啟檔案 obj get gca,children 資料儲存在obj裡面 data get obj,cdata open figname.fig h line get gca,children get linehandles xdata get h line,xda...