Pacbio三代全长转录组(Iso-Seq ) 进行转录本与基因组的比对

Pacbio三代全长转录组(Iso-Seq ) 进行转录本与基因组的比对

Pacbio三代全长转录组进行数据分析,得到的转录本,可以同基因组进行比对(如果有基因组的话),比对的软件很多,推荐采用GMAP。

其使用方法如下:

gmap -D /share/nas1/zhangqx/testing/pacbio/ref/gmap -d ipoBat4  -f samse -n 0 -t 16  --cross-species --max-intronlength-ends 200000 -z sense_force   hq.fastq > hq_isoforms.fasta.sam  2> hq_isoforms.fasta.sam.log

常用参数说明如下:

-D : 指定参考基因组的gmap 建库目录

-d : 指定gmap 建库中的名称

-f : 指定输出文件的格式

-n:  输出结果中允许path 最多的数目

-t : 线程数量

--cross-species:以更加灵敏的算法进行可剪切比对,特别适合在物种间比对 

--max-intronlength-ends: 第一个和最后一个内含子的最大间距

-z: 指定转录本的方向



  • 发表于 2018-08-10 15:40
  • 阅读 ( 6523 )
  • 分类:三代测序

0 条评论

请先 登录 后评论
microRNA
microRNA

115 篇文章

作家榜 »

  1. omicsgene 700 文章
  2. 安生水 348 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. 红橙子 78 文章
  7. rzx 75 文章
  8. CORNERSTONE 72 文章