|
From: Sébastien B. <seb...@ul...> - 2011-07-13 20:29:58
|
I think you generated too few reads. That is why it fails. 1 302 2 422 3 982 4 1550 5 1804 ********** peakCoverage 6 1596 7 1220 8 1004 9 666 10 344 11 202 12 126 13 84 14 20 You need to get the peak coverage at least above 20. How did you generate the reads ? Estimated genome length: 5161 nucleotides is correct > Okay, so I've downloaded sebhtml/ray 1.6.1-rc3, and have got as far as the > first smoothing bug I found. It doesn't seem to work properly with any Ray > that I have tested. > > Coverage: http://pastebin.com/egRf1dzj > > Note that the analysis doesn't match the coverage, and doesn't agree with > itself (min/max coverage is 1, but Percentage of vertices with coverage 1: > 2.92579 %). > > Input (exact reads generated from the phiX genome): > http://pastebin.com/yae0xNKC > > Output: http://pastebin.com/tjrj5W2q > > -- David Eccles (gringer) Sébastien |
