Hi, I have two fasta files that I want to align together. One (A) contains fragments of a sequence that is created by a sliding window step size of 1. The other one (B) is multiple sequences. I want to align A with B and use the output (can I only get one output?) to label the fragment in A (to show that it was mapped to a sequence in B). I tried to use bowtie 2 but apparently, I have to use a fastq file. I also wanted to try emboss water but I am not sure if I can use it for multiple sequences for both files.
Is there anyway I can achieve this?
1 answer
If the sequences are similar (and or a similar length) you could try Jim Kent's blat. It will take multiple sequences in both files and produce results for all comparisons by a single command line. You can easily parse tabular output should you choose to do do.
You can use -f option with bowtie2 to specify that input reads are fasta. You don't need to use fastq files. You can use it since you will be creating an index with one of the files (with multiple sequences) and then aligning the sequences to this index one at a time from query file.
Log in to answer this question.
-
Comparison of bed files and output matrix
written by 1934822857 •
-
sequence similarity in R
written by kng •
-
Can you force two sequences to align with pairwise alignment?
written by msening •
-
Bowtie paired-end mapping: insert size issue
written by biostart
-
Multi-fasta sequence file alignment using sliding window
written by susheelbhanu •
-
Aligning multiple Fastq files using Bowtie2
written by deepak18 •
-
How can I map reads only of a certain length?
written by a.rex
-
Sliding Window : discriminate partial and fully contained fragments (from a bam file)
written by w.elmselmi •
-
Aligning several sequences with EMBOSS water program
written by dolevrahat •
-
Align short reads to multiple reference sequence
written by Sandeep
