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

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

ADD REPLYlink written 12 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: 1072 users visited in the last hour