Closed duhuipeng closed 3 years ago
Hello @duhuipeng,
Please use the inherited hapmers instead of the full maternal / paternal dbs. See the section in here: https://github.com/marbl/merqury/wiki/1.-Prepare-meryl-dbs#3-build-hap-mer-dbs-for-trios
I'll close the issue opened in meryl repo.
Thanks, Arang
I don't understand. Do you mean my input file becomes as follows when I meryl union the patrilineal,Should I input it like this ? /home/duhuipeng/DHP/merqury-1.1/merqury.sh /home/duhuipeng/DHP/DRC/DRC.SON.k21.meryl/ union.k21.meryl ../father_purged.fasta ../purged_mo.fasta test-2 Do work like this ?It's feel strange
No, I meant to prepare the hapmer dbs first.
Make a directory and symlink your maternal and paternal meryl dbs, and run $MERQURY/trio/hapmers.sh
.
ln -s /path/to/your/DRC.mother.k21.meryl Mother.meryl
ln -s /path/to/your/DRC.father.k21.meryl Father.meryl
ln -s /path/to/your/DRC.SON.k21.meryl Son.meryl
sh $MERQURY/trio/hapmers.sh Mother.meryl Father.meryl Son.meryl
And use the inherited, filtered hapmers (they will be symlinked as Mother.hapmers.meryl and Father.hapmers.meryl).
Given the meryl db size, I assume the maternal.meryl
and paternal.meryl
were directly collected using meryl from the full sequencing data set.
The hapmers.sh
produces the inherited kmer spectrum, as described in the wiki page and in the Merqury paper.
The hapmers are haplotype specific kmers, which you'd want to evaluate the assembly against. So, after producing the hapmer.meryl dbs, run
./merqury.sh Son.meryl Mother.hapmers.meryl Father.hapmers.meryl father_purged.fasta purged_mo.fasta test
I got it ,Thanks
Great, let me know when you encounter other issues.
Dear author I guess if your place is reversed, it should correspond to my meryl file
Dear author When I run merqury, the second step: hap_blob.sh ,the problem also arises This is the code I run: /home/duhuipeng/DHP/merqury-1.1/merqury.sh Son.meryl Mother.hapmers.meryl Father.hapmers.meryl purged_mo.fasta father_purged.fassta test
This is all I have generated when running the command as above
Dear author I saw his error saying it couldn't connect Mother.hampers.meryl,saying it already exists I looked at hapmer.sh code again, it generated files should be hamper.Father.meryl files,It shouldn't be hampers.Father.meryl . I feel like my input should be ./merqury.sh Son.meryl Mother.hapmer.meryl Father.hapmer.meryl father_purged.fasta purged_mo.fasta test I see what you're giving ./merqury.sh Son.meryl Mother.hapmers.meryl Father.hapmers.meryl father_purged.fasta purged_mo.fasta test Do you feel that this input is the reason Is the error caused? Looking forward to your reply
Yes, you are correct. Use Mother.hapmer.meryl
and Father.hapmer.meryl
instead.
Apologies about the confusion.
Dear author this is my code : /home/duhuipeng/DHP/merqury-1.1/merqury.sh ../ALL.k21.meryl/ ../Mo_ALL.k21.hapmer.meryl/ ../Fa_ALL.k21.hapmer.meryl/ ../../DrcSonPhased_hap2.fasta ../../DrcSonPhased_hap1.fasta DRC But the above problem has emerged ,I guess it is insufficient memory leading to core dump. I would like to consult what parameters can be set to limit it memory size and make it work properly
Best huipengDu
When I run merqury, the second step: hap_blob.sh ,the problem arises This is the code I run: /home/duhuipeng/DHP/merqury-1.1/merqury.sh /home/duhuipeng/DHP/DRC/DRC.SON.k21.meryl/ ../../DRC.father.k21.meryl/ ../../DRC.mother.k21.meryl/ ../father_purged.fasta ../purged_mo.fasta test-2 the error as follow in hap_blob at first I thought core dump was out of memory, but I found that I couldn't run when I used the size of 3T memory. How can I solve this problem? Because I want to draw a picture of its phase separation,But this step is wrong is not drawn, I hope to get your advice.
Looking forward to your reply! Best