IGV and IGVtools usage

Create your own .genome file

  1. Click Genomes>Create .genome File in menu bar of IGV. IGV will popup a window for entering needed information.

  2. Enter an ID and a descriptive name for the genome.

  3. Enter the path on your file system or a web URL to the FASTA file for the genome. If the FASTA file has not already been indexed, an index will be created during the import process. This will generate a file with fai suffix in the same directory as the FASTA file; Or one can index genome FASTA file using samtools faidx genome.fasta.

  4. Optionally, specify the cytoband file and the annotation (gene) file.
  5. If the sequence (chromosome) names differ between your FASTA and annotation files, you might need to create an alias file to provide a mapping between the different names.  Certain well-known aliases are built into IGV and do not require an alias file. These include mappings that involve adding or removing the prefix chr to the name, for example  1 > chr1 and chr1 > 1.  Also, NCBI identifiers of the form  gi|125745044|ref|NC_002229.3| in a FASTA file will be mapped to names of the form NC_002229.3 in the corresponding GFF file.

  6. Click Save and select the directory in which to save the genome archive (*.genome).

Transfer BAM file to tdf or wig (Index BAM first)

igvtools count --strands first --includeDuplicates rnaseq.bam rnaseq.tdf,rnaseq.wig mm9
igvtools count --strands first --includeDuplicates --pairs -w 1 rnaseq.bam rnaseq.tdf,rnaseq.wig mm9

IGV批量截图时包含红色标记区域

  • 自动截图前,先把要截图的区域导入,然后再运行batch script
CHENTONG
版权声明:本文为博主原创文章,转载请注明出处。
alipay.png WeChatPay.png

CHENTONG

CHENTONG
积微,月不胜日,时不胜月,岁不胜时。凡人好敖慢小事,大事至,然后兴之务之。如是,则常不胜夫敦比于小事者矣!何也?小事之至也数,其悬日也博,其为积也大。大事之至也希,其悬日也浅,其为积也小。故善日者王,善时者霸,补漏者危,大荒者亡!故,王者敬日,霸者敬时,仅存之国危而后戚之。亡国至亡而后知亡,至死而后知死,亡国之祸败,不可胜悔也。霸者之善著也,可以时托也。王者之功名,不可胜日志也。财物货宝以大为重,政教功名者反是,能积微者速成。诗曰:德如毛,民鲜能克举之。此之谓也。

生信宝典文章集锦

生信的作用越来越大,想学的人越来越多,不管是为了以后发展,还是为了解决眼下的问题。但生信学习不是一朝一夕就可以完成的事情,也许你可以很短时间学会一个交互式软件的操作,却不能看完程序教学视频后就直接写程序。也许你可以跟着一个测序分析流程完成操作,但不懂得背后的原理,不知道什么...… Continue reading

狗年话狗,也话GO

Published on February 20, 2018

network3D 交互式网络生成

Published on February 18, 2018