Hi everyone, I need to align scaffolds, CDS and transcripts (local alignments) and see the aligned/nonaligned and the coding and noncoding regions. I have a fasta file with one transcript, CDS and multiple scaffolds and I need to allign all of those. Could you please suggest me any program that I can use? I heard bowtie can be used in this case?? Thanks a bunch!
•
written
by
MAPK
2 answers
Hello,
I suppose you have your reasons for aligning CDS features to scaffolds.
I would suggest bwa. I am currently working with it, and find the results acceptable.
There is a short intro by a Broad member arguing BWA being the balance between performance and accuracy. https://www.broadinstitute.org/files/shared/mpg/nextgen2010/nextgen_li.pdf
I suggest indexing your scaffolds, and aligning CDS features to scaffolds as opposed to the other way around.
Adrian
For multisequence alignment, just have a look at this page for a number of options.
I don't think I've ever seen multiple sequence alignment programs use global alignment, it wouldn't make much sense. I'm assuming that you know what multiple sequence alignment actually is and how it differs from what bowtie/bwa/etc. do, though perhaps that's an incorrect assumption.
Log in to answer this question.
-
Can you suggest tools or script to get canonical sequences from a cds fasta file?
written by Yuto •
-
Program to align scalffold to scalffold?
written by carlosgonzalezcruz327 •
-
can lnc-RNA target coding RNA?
written by ta_awwad
-
How do I extract aligned reads from BAM file?
written by MAPK
-
Aligning reads to genome/transcriptome
written by EVR
-
Gene prediction tool
written by MAPK
-
How Do I Align More Than Five Long Dna Scaffolds (Contigs) To The Reference Gene?
written by biostar •
-
Bowtie Assembly Parameters
written by biolab
-
How One Can Separate Noncoding From Coding Ucsc And Ensembl Transcripts
written by biorepine
-
Blast
written by pjyoti29 •
Are you trying to get a multisequence alignment? Is there a reason you mentioned bowtie? It's not normally used to align contigs (that's not really its forte).
@dpryan79, Thanks for replying to my thread. Yes I am trying to do multisequence alignment and I just mentioned Bowtie because someone told me to try Bowtie. Is there any program you know I could use in this case?
Are you trying to align short read data (RNA-Seq, etc.) to a reference sequence? What exactly are you trying to do? It sounds like you are confused about what you are doing.