Is it possible to extract Coverage track data to a table?

58 views
Skip to first unread message

jspai...@otogenetics.com

unread,
Jul 5, 2018, 1:30:10 PM7/5/18
to igv-help
I am trying to extract the Coverage track data including total count, breakdown of nucleotides, and forward/reverse read breakdown from the IGV interface.

This information appears in a pop up box when you hover over a column in the coverage track.

I am using v2.4.10 of the interface on a Windows 10 OS.  I have not used the command line version but i am not adverse to doing so if this would help.

I have tried using the count command in IGVtools with a  window size of one. This gives me a .wig file with the total count information per bp with out any of the breakdown. 

I have also tried using both samtools and bcftools with no success but any suggestions would be very helpful.

Is there another way to export this information into a table?  Is it also possible to take a larger file and export select ranges of the coverage information by referencing a .bed file?

Thank you, 
-John Spainhour

James Robinson

unread,
Jul 5, 2018, 3:17:11 PM7/5/18
to igv-help
Look at the documentation for the command line igvtools,  there are flags for the count command to create ".wig" type files with the multiple columns.

--

---
You received this message because you are subscribed to the Google Groups "igv-help" group.
To unsubscribe from this group and stop receiving emails from it, send an email to igv-help+unsubscribe@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/igv-help/b6d6f9c1-01e9-45fb-81ac-d8e45dc6f017%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply all
Reply to author
Forward
0 new messages