SHZ66 / zinba

Automatically exported from code.google.com/p/zinba
1 stars 1 forks source link

Mappability files for other read lengths? #19

Closed GoogleCodeExporter closed 8 years ago

GoogleCodeExporter commented 8 years ago
Hi, 

this is not really an error/issues report, but more like a question. 

Seeing as how all the mappability files you've provided are for read lengths 36 
and 50, it would appear that it is for Illumina GA2 generated datasets only? 
What if I have HiSeq data (read length 100), or other types of sequencing data 
of other read lengths? 
Is there a way I can generate my own mappability files so that ZINBA can be 
used on these datasets?

Thanks

Original issue reported on code.google.com by angela...@gmail.com on 10 Apr 2012 at 7:42

GoogleCodeExporter commented 8 years ago
It is possible to create your own mappibility files but they take quite a while 
to generate (especially if your genome is large), here are some details: 
http://www.biostars.org/post/show/14092/generating-mappability-files/

Original comment by daiyi...@gmail.com on 27 Apr 2012 at 1:27

GoogleCodeExporter commented 8 years ago
Yes, it is possible to generate your own mappability maps using code from the 
peakSeq software from the gerstein lab available here: 
http://info.gersteinlab.org/PeakSeq .  Essentially the output that you need is 
the *b.out files for each chromosome placed in a directory of your choice with 
no *b.out files from other genomes, etc.

Original comment by homer...@gmail.com on 28 May 2012 at 6:26

GoogleCodeExporter commented 8 years ago

Original comment by homer...@gmail.com on 28 May 2012 at 6:27