• 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 

  • 相关阅读:
    在 MAC 下配置 Nginx
    Color Schema 配色随笔
    .Net与 WebAssembly 随笔
    关于Xamarin、Qml、数据绑定、MVC、MVVM 相关的散讲
    用Nuget部署程序包
    Qt3D
    Qt3D Shader
    Qt QML 2D shader
    LearnOpenGL
    Qt3D 5.9 and future
  • 原文地址:https://www.cnblogs.com/renping/p/7852048.html
Copyright © 2020-2023  润新知