售后

编辑本页后请点击“保存”。请参阅 syntax 了解维基语法。只有在您能改进该页面的前提下才编辑它。如果您想尝试一些东西，请先到 playground 热身。

=======hisat2 mapping参数=======

转录有参、医口流程hisat2 mapping命令如下：

hisat2 -x <hisat2 index> -p 4 --dta -t --phred33 -1 sample_1.clean.fq.gz -2 sample_2.clean.fq.gz --un-conc-gz sample.unmap.fq.gz 2> sample_align.log | samtools sort -O BAM --threads 4 -o sample.bam -
  
  均为默认参数，hisat2版本version 2.0.5

=======hisat2 mapping log=======
  Time loading forward index: 00:00:31
  Time loading reference: 00:00:06
  Multiseed full-index search: 00:13:05
  21501684 reads; of these:
    21501684 (100.00%) were paired; of these:
      2675430 (12.44%) aligned concordantly 0 times
      15080209 (70.14%) aligned concordantly exactly 1 time
      3746045 (17.42%) aligned concordantly >1 times
      ----
      2675430 pairs aligned concordantly 0 times; of these:
        234828 (8.78%) aligned discordantly 1 time
      ----
      2440602 pairs aligned 0 times concordantly or discordantly; of these:
        4881204 mates make up the pairs; of these:
          3031909 (62.11%) aligned 0 times
          1359660 (27.86%) aligned exactly 1 time
          489635 (10.03%) aligned >1 times
  92.95% overall alignment rate
  Time searching: 00:13:12
  Overall time: 00:13:43

总的mapping计算如下
    92.95% =((15080209+3746045)*2 + 234828*2 + 1359660 + 489635)/(21501684*2)
  软件输出的unmap reads为 2675430 (12.44%) aligned concordantly 0 times

编辑摘要

当您选择开始编辑本页，即寓示你同意将你贡献的内容按下列许可协议发布： CC Attribution-Share Alike 4.0 International

售后

用户工具

站点工具

页面工具