Benjamin Redelings
unread,Jun 19, 2020, 1:21:49 PM6/19/20Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to Benjamin Schlossman, bali-ph...@googlegroups.com
Hi Ben,
BTW, would you mind sending questions to the google group, so that other
people can benefit from the answers?
1. I think you should probably use the alignment-max command instead of
the alignment-consensus command.
2. You can do something like this, using pipes:
cut-range --skip=1000 run-1/C1.P1.fastas run-2/C1.P1.fastas |
alignment-chop-internal -N 25 | alignment-max > p1-wpd.fasta
The cut-range command is almost like "cat", but it allows you to skip
the first (say) 1000 iterations. The alignment-chop-internal command
removes internal sequences, but you have to specify the number of leaf
sequences to keep.
3. If you just run `bp-analyze run-1 run2`, you should get a nice
report. It should include an alignment based on samples from both runs,
constructed using alignment-max. If you look in Results/command.log, it
records the command that are run. However, with pipes, you might only
see the last command in the pipeline.
-BenRI