Hi Gary,

narrowPeak files can be visualized directly on the UCSC Genome Browser as a custom track, no need to convert to bigBed first. However if you wanted to, example 3 on this page: http://genome.ucsc.edu/goldenPath/help/bigBed.html is a good example of how to convert a non-standard bed file to a bigBed file, in that you need to supply bedToBigBed with an autoSql file that describes your data. For more information see this question from our mailing list archive: https://groups.google.com/a/soe.ucsc.edu/d/msg/genome/9PXjH2mlqrE/MrBs3pZ9WLEJ

That being said, I'm slightly confused by the bedToBigBed command you're running, what is FemaleOnly_H3K4me3.bed? Did that come from macs2? According to the macs2 documentation here: https://github.com/taoliu/MACS it appears that this file can be loaded directly onto the Genome Browser. Have you tried that already? I also found this issue: https://github.com/taoliu/MACS/issues/123, which seems to be unresolved, perhaps you could ask there for more information? Lastly, since the number in the score field only determines the darkness of the peak in the genome browser, perhaps you could write a little awk script to check if the values are over 1000 and if so, convert them to 1000.

It would be nice if you could send follow up questions to our mailing list at: https://groups.google.com/a/soe.ucsc.edu/forum/#!forum/genome, this way our entire team will see your question and can help you out, in addition to all the people that follow that list and may not see the question here.

Thanks,

ChrisL from UCSC