=======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