Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ROSE usage for mm39 genome #2

Open
Brahmandam-Gayatri opened this issue Nov 12, 2024 · 0 comments
Open

ROSE usage for mm39 genome #2

Brahmandam-Gayatri opened this issue Nov 12, 2024 · 0 comments

Comments

@Brahmandam-Gayatri
Copy link

Brahmandam-Gayatri commented Nov 12, 2024

Hello,

Thank you for ROSE for calling super-enhancers. I have a query: How can I make the refseq table for mm39 and use it for stitching the super-enhancers of my ChIP-seq data?

In addition, I am getting the following error:

MAKING TSS COLLECTION
12900
9000
13000
13100
13200
13400
9100
13300
9200
9000
13500
9300
9400
9100
13600
13700
9200
13800
9300
13900
9400
14000
13400
14100
13500
13600
13700
13800
13900
14000
14100
9500
9600
9700
9800
14200
14300
14400
14500
14600
14200
14300
14400
14500
14600
9500
9600
9700
9800
Traceback (most recent call last):
  File "/path/to/file/ROSE/bin/ROSE_geneMapper.py", line 299, in <module>
    main()
  File "path/to/file/ROSE/bin/ROSE_geneMapper.py", line 282, in main
    enhancerToGeneTable,geneToEnhancerTable,withGenesTable = mapEnhancerToGene(annotFile,enhancerFile, uniqueGenes=True, byRefseq=options.refseq, subtractInput=options.control, transcribedFile=transcribedFile)
                                                             ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/path/to/file/ROSE/bin/ROSE_geneMapper.py", line 140, in mapEnhancerToGene
    closestGene = allEnhancerGenes[distList.index(min(distList))]
                  ~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
IndexError: list index out of range

Command used:

ROSE_main.py -g MM39 -i /path/to/file/treated1.bed -r /path/to/file/treated1.sorted.bam -c /path/to/file/treated1_input.sorted.bam -o . -s 12500 -t 2500

I have tried to make a refseq file for mm39 genome as specified in the instructions. Hence any help in this regard is highly appreciated.

Thanks & regards,
Gayatri.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant