This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Run The Soapaligner?

In the SOAP website (http://soap.genomics.org.cn/soapaligner.html#evalu2) there is the following path to run data:

Format reference sequence: <ExecutablePath>/2bwt-builder <FastaPath YourFasta="">

./2bwt-builder ~/human_genome.fa

Then under the directory there will be 13 index files, all their prefixes are yourfasta file name with “.index” added, e.g. humangenome.fa.index. The suffixes include *.amb, *.ann, *.bwt, *.fmv, *.hot, *.lkt, *.pac, *.rev.bwt, *.rev.fmv, *.rev.lkt, *.rev.pac, *.sa, and *.sai.

Alignment quick start

For alignment of single-end reads:

./soap –a <reads_a> -D index.files> -o output/output>

I got the 13 index files as described above (1), but I am not being able to set up the second command. Do I put all the 13 index files in <index.files> or that is actually the path for those files?, Do I randomly create my output files? I have 10 files for each sample, can I put them one after another in <reads_a> ? so that it runs all of them at once?.

Sorry if these are simple questions.

deleted-post

0 answers

No answers yet.

Log in to answer this question.