In case of using [Mauve], we need to prepare a comparison file, or a crunch file. To generate the file, one can upload the genome sequences in genbank/embl/fasta format to WebACT, or upload sequences in fasta format to Double ACT (It has been reported that genbank format didn’t work well with DoubleACT1), then wait for their notification when finished.
Alternately, one can generate a crunch file locally using blastall like this:
#!/bin/bash if [ $# -ne 4 ]; then echo -e "\nUsage: $0 [blastn/tblastx] [input_first_genome.fna] [input_second_genome.fna] [output_filename]\n" exit 1 fi # run formatdb formatdb -p F -i $2 # do blastall with 30 cores blastall -p $1 -m 8 -e 1e-5 -d $2 -i $3 -o $4 -a 30
Besides, ACT works just like Artemis, except gbk file, a gff file can be loaded as well for annotations.
To sum up, here is the steps need to be done:
- get genome sequences in fasta format
- do pairwise blastn search to generate the comparison crunch files.
- load the genome fasta and crunch files
- load the annotation files in gff format
- https://email@example.com/msg00372.html ↩