Hi,
I am new to this and I hope someone can help. I have pooled sequencing data that I am trying to analyse using Galaxy. I've done quite a bit of online searching and it seems that FreeBayes should be able to do this, if I select "set population", click the "Assume that samples result from pooled sequencing" option and change the ploidy to nx2 (number of alleles, where n is the number of subjects and the organism is diploid).
However, whenever I do this I get an error: usually just "Killed"
I was originally setting my polidy rather high (190 as I have 95 subjects pooled), so I wondered if this was the problem, however, it fails if I do a ploidy of only 4 too. I've tried various things to try to see where I am going wrong:
All with the same BAM file:
Set population model options: Do not set --> works
Set population model options: set, Assume that samples result from pooled sequencing: not ticked, Default ploidy for the analysis: 2 --> works
Set population model options: set, Assume that samples result from pooled sequencing: ticked, Default ploidy for the analysis: 2 --> works
Set population model options: set, Assume that samples result from pooled sequencing: ticked, Default ploidy for the analysis: 4 --> fails (killed)
Set population model options: set, Assume that samples result from pooled sequencing: ticked, Default ploidy for the analysis: 10 --> fails (killed)
Set population model options: set, Assume that samples result from pooled sequencing: ticked, Default ploidy for the analysis: 10 --> fails (killed)
It seems that it is the ploidy part that I am doing wrong, as it works with pooled data but ploidy of 2. I'm sure I have to change the ploidy though, or else how does the program know how many subjects are in the pool? Also, everywhere that I've ready says you have to change the ploidy.
I apologise if my question is naive. As I said, I am new to Galaxy and this is the first thing I am trying to do!
Any help / suggestions would be appreciated,
Thanks,
Nic