igvteam / igv-reports

Python application to generate self-contained pages embedding IGV visualizations, with no dependency on original input files.

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

igv-reports doesn't display bam but igv does

azbard opened this issue · comments

Absolute noob here - be gentle.

Is there any conceptual reason why a bam/bed combo would display fine in igv proper (java) but only the bed shows up in a report generated as follows?

create_report bedfile.bed hg19.fasta --tracks bamfile.bam

These two made with the same files...
Screenshot 2021-12-14 151244
image

[edited to show matching locations in igv and igv-reports]

In your igv-reports screenshot I see chr20 is displayed, but in the IGV desktop screenshot chr5. If you are displaying the chr5 data at the chr20 location then you won't see any alignments. I think this is the root cause, but its difficult to say how this could occur without having the files. If you're able to share the report (generated html) you can email it privately to igvteam (at) broadinstitute.org and I will have a look.

Thanks for the data. I think your bed file might be malformed, based on the table column names in your screenshot above. Could you zip and attach it (just the bedfile).

Actually the table column names is a report bug, it would still be useful to see the bed file. You can zip and attach it here, or email to the igvteam email you sent the report to.

I figured out the report not showing the tracks - my pysam library was causing the problem - it had been altered for another project. Ugh - this was all my own fault.

The bed though...
bedfile-copy.zip

Thank you for your support and prompt response to our issue. We were able to get igv-reports to function as desired. Really love this tool.
Feel free to close this issue.