Question: How to export the SNPs between whole genome alignments using the command line?
0
gravatar for Paul
7 months ago by
Paul70
India
Paul70 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 7 months ago by h.mon24k • written 7 months ago by Paul70
3
gravatar for h.mon
7 months ago by
h.mon24k
Brazil
h.mon24k 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 7 months ago by h.mon24k

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

ADD REPLYlink written 7 months ago by Paul70
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: 1020 users visited in the last hour