Question: How to export the SNPs between whole genome alignments using the command line?
0
gravatar for Paul
19 months ago by
Paul80
India
Paul80 wrote:

I am using Mauve to align two the whole genome sequences, one of the sequence is in the form of contigs. Another one is in the form of a complete genome.

I could align the two sequences using the following command

progressiveMauve ref.fa contigs.fa --output alignment.xmfa

And to extract the SNPs, I tried the solution from https://bioinformatics.stackexchange.com/questions/3658/can-i-export-the-snps-between-whole-genome-alignments-using-the-command-line/3664#3664

But while exporting the classpath path it shows an error

export CLASSPATH="$(find "$MAUVE_DIR" -name \*.jar -print0 | tr '\0' :)$CLASSPATH"
find: ‘‹/home/me/Downloads/mauve_snapshot_2015-02-13/›’: No such file or directory

I have downloaded the mauve file and extracted it in the following location

/home/me/Downloads/mauve_snapshot_2015-02-13/

I am using a linux machine. How can I solve this issue?

ADD COMMENTlink modified 19 months ago by h.mon29k • written 19 months ago by Paul80
3
gravatar for h.mon
19 months ago by
h.mon29k
Brazil
h.mon29k wrote:

How did you define $MAUVE_DIR? The error message says:

find: ‘‹/home/me/Downloads/mauve_snapshot_2015-02-13/›’: No such file or directory

Which means you are defining $MAUVE_DIR as ‹/home/me/Downloads/mauve_snapshot_2015-02-13/›, but it should be /home/me/Downloads/mauve_snapshot_2015-02-13/ instead.

ADD COMMENTlink written 19 months ago by h.mon29k

Thanks a lot @h.mon :) it worked...

ADD REPLYlink written 19 months ago by Paul80
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1646 users visited in the last hour