zoukankan      html  css  js  c++  java
  • 2、Tophat align_summary.txt and samtools flagstat accepted_hits.bam disagree

    ###https://www.biostars.org/p/195758/

    Left reads:
              Input     :  49801387
               Mapped   :  46258301 (92.9% of input)
                of these:   1703360 ( 3.7%) have multiple alignments (103886 have >20)
    Right reads:
              Input     :  49801387
               Mapped   :  45542088 (91.4% of input)
                of these:   1680180 ( 3.7%) have multiple alignments (103887 have >20)
    92.2% overall read mapping rate.
    
    Aligned pairs:  44100265
         of these:   1622622 ( 3.7%) have multiple alignments
                      261400 ( 0.6%) are discordant alignments
    88.0% concordant pair alignment rate.
    1、解释计算原理:(44100265 - 261400) / 49801387 = 88%` ;
    1622622/1622622 = 3.7% ;
    2、If I needed to report a number, I'd report the 92.2% overall alignment rate. You can probably get something similar with samtools view -c -F 260 accepted_hits.bam.(得到输入left和
    right的比对数目)/
    49801387 *2
    3、samtools flagstat accepted_hits.bam 

  • 相关阅读:
    elasticsearch api
    kaili camera
    mysql create db utf8 character
    npm run-script
    d-link kvm 关闭声音
    setInterval js
    jpa datasource config
    mvn添加本地jar
    Sublime Text 2 中文包
    初遇ping++
  • 原文地址:https://www.cnblogs.com/renping/p/7852048.html
Copyright © 2011-2022 走看看