Skip to content

Commit

Permalink
added new param
Browse files Browse the repository at this point in the history
  • Loading branch information
mahulchak committed Jul 6, 2016
1 parent ce29e56 commit c927a50
Showing 1 changed file with 3 additions and 1 deletion.
4 changes: 3 additions & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -65,7 +65,7 @@ The package contains all necessary components to run quickmerge. We also provide
```
Finally, use 'quickmerge' to merge the two assemblies (note: the order of the self and hybrid assembly is important:
```
quickmerge -d out.rq.delta -q hybrid_assembly.fasta -r self_assembly.fasta -hco 5.0 -c 1.5 -l n
quickmerge -d out.rq.delta -q hybrid_assembly.fasta -r self_assembly.fasta -hco 5.0 -c 1.5 -l n -ml m
```
Description of the parameters:

Expand All @@ -76,6 +76,8 @@ The package contains all necessary components to run quickmerge. We also provide
For both "hco" and "c", bigger the number, more stringent is the criteria for contig selection (which will lead to fewer contigs being merged). If they are too small (<1), chances of spurious merging will increase.

-l: controls the length cutoff for anchor contigs. A good rule of thumb is to start with the N50 of the 'self_assembly.fasta'. E.g. if the N50 of your self_assembly.fasta is 2Mb. Then use 2000000 as your cutoff. Lowering this value will lead to more merging but may increase the probability of mis-joins.

-ml: controls the minimum length of the alignment to be considered for merging. This is especially helpful for repeat-rich genomes. Default is 0 but 5000 is a good starting point.

4. SOME HELPFUL TIPS:

Expand Down

0 comments on commit c927a50

Please sign in to comment.