Problem in alignment with soap
0
0
Entering edit mode
6.1 years ago

Hi all

I'm trying to align probe sequence of rat230_2 micoarray to the transcriptome of rattus with utmost 6 mismatch. first I've build the index file from the transcriptome with the following command:

./Aligning/2bwt-builder ./Aligning/transcriptom.fasta

then I've tried to align the probe sequence to the index file with following command:

./Aligning/soap -a ./Aligning/probes.fasta -D ./Aligning/transcriptom.fasta.index -o output.txt -v 6 -p 4 -r 2

but i got the following output in the command line:

length y < 0, skip anyway ...

and the percent of alignment is 0.

where is the problem?

thanks,

Mansoor,

soap • 1.9k views
ADD COMMENT
0
Entering edit mode

I think you misuse "./"

You only need to use this for an execution file, e.g.:

./Aligning/2bwt-builder /Aligning/transcriptom.fasta

./Aligning/soap -a /Aligning/probes.fasta -D /Aligning/transcriptom.fasta.index -o output.txt -v 6 -p 4 -r 2
ADD REPLY
0
Entering edit mode

yes you are right but the problem is still exist.

ADD REPLY
0
Entering edit mode

Are you sure soap works with FASTA instead of FASTQ?? I am not sure, I have never tried it.

Try to see if you get unmapped reads with -u option.

ADD REPLY
0
Entering edit mode

Yes it is work. I found the problem. it doesn't work for read with length less than 27. so I used bwa aligner.

ADD REPLY

Login before adding your answer.

Traffic: 927 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6