統計tophat map上的read數量

2021-09-07 02:32:43 字數 1007 閱讀 8972

samtools flagstat/sra111111/srr111222/accepted_hits.bam

78406056 + 0 paired in sequencing (3)

39915264 + 0 read1 (4)

38490792 + 0 read2 (5)

4805744 + 0 singletons (6.13%:-nan%) (8)

(2)=(7)+(8)

(3)=(4)+(5)

usage: samtools flagstat 

$ samtools flagstat example.bam

11945742 + 0

in total (qc-passed reads + qc-failed reads) #總共的reads數

0 + 0

duplicates

11945742 + 0 paired in

sequencing #有多少reads是屬於paired reads

5972871 + 0

read1 #reads1中的reads數

5972871 + 0

read2 #reads2中的reads數

6412042 + 0 properly paired (53.68%:-nan%) #完美匹配的reads數:比對到同一條參考序列,並且兩條reads之間的距離符合設定的閾值

6899708 + 0

636656 + 0 singletons (5.33%:-nan%) #單獨一條匹配到參考序列上的reads數,和上乙個相加,則是總的匹配上的reads數。

469868 + 0

) #同上乙個,只是其中比對質量》=5的reads的數量

samtools view  ./accepted_hits.bam  | cut -f1 | sort | uniq | wc -l

ref:

規則和統計(上)

科學家們原本以為隨著nl語法概括得越來越全面,同時計算機計算能力的提高,這種方法可以逐步解決nl理解的問題。但是很快遇到了麻煩,要想通過文法規則覆蓋哪怕20 的真實語句,文法規則的數量也至少是幾萬條。語言學家幾乎已經來不及寫了,而且這些規則到後來甚至會出現矛盾,為了解決矛盾,還要說明各個規則的使用環...

github上的提交數統計

之前在github上使用一直沒有注意乙個問題,就是我最近頻繁的提交,但是github上沒有顯示我提交頻繁程度。我明明每天都有提交,這個是什麼原因?於是我查閱了一下,github上對於貢獻統計的要求,一般來說都有以下幾點了 1.fork是不能算的,你必須是你提交的這個專案中的成員 2.你提交的pull...

關於Jedis無法連線上Linux上Redis問題

環境 centos7 redis 主要解決jedis客戶端無法連線linux上redis服務問題 1 修改redis目錄下的redis.conf配置檔案 如果未設定redis認證密碼,則需要設定保護模式為no protected mode no 2 停止centos7防火牆或設定規則 這裡使用停止方...