Running alignment using Bowtie by setting specific maximum number of mismatches
1
0
Entering edit mode
3 months ago

Hi everyone, I need to run blast process on a sequence file by specifying some special parameters. here are the parameters: 1- word size of 6 2- allowed maximum 2 mismatch (on each side)

I'm using ncbi blast to do the alignment using the following command: blastn -query seq.fasta -db ref -out results.txt -num_threads 8 -outfmt '6 qseqid sseqid sseq' -word_size 6

Since I couldn't define maximum mismatch for the mentioned command, I decided to use bowtie to do the alignment process. would you please tell me the bowtie command equivalent to ncbi blast command which I mentioned with considering the maximum allowed number of mismatches?

Regards.

mismatch alignment bowtie • 477 views
ADD COMMENT
0
Entering edit mode
ADD COMMENT
0
Entering edit mode

I don't have a sam or bam file! the file format is fasta.

ADD REPLY
0
Entering edit mode

but bowtie generates a sam/bam file; use samtools to filter downstream bowtie.

ADD REPLY
0
Entering edit mode

The file which I'm running alignment on it, is fasta. does bowtie support fasta files as input?

ADD REPLY

Login before adding your answer.

Traffic: 2243 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