Does Alignment Require All The Reads With The Same Length?
3
0
Entering edit mode
11.3 years ago
Tonyzeng ▴ 310

HI, My question is as my title, also is that any tools to help me get the same length of reads?

alignment • 2.6k views
ADD COMMENT
2
Entering edit mode
11.3 years ago
Ido Tamir 5.2k

It depends on the aligner, but the popular aligners I know (bowtie, bwa) don't require same length of reads (an early version of tophat required this). Maybe eland (CASAVA) with its hashing of reads requires the same length. For processing of reads there are many tools available e.g. http://hannonlab.cshl.edu/fastx_toolkit/commandline.html#fastx_trimmer_usage fastx trimmer

ADD COMMENT
0
Entering edit mode

Also, if you feel uncomfortable using command line tool, you can trimm reads using online tools like galaxy http://galaxyproject.org/

ADD REPLY
0
Entering edit mode
11.3 years ago
Tonyzeng ▴ 310

Actually, I tried Galaxy to cut my overrepresent sequence, but my sequence is on 5' of read and reported as Illumina Single End PCR Primer 1. So I used Galaxy fastx-clipped but it can only clip adapt or primer sequence in 3'. So I changed to use cutadapt tool and cut 5' of read, it works well. However, I got a wide variety of length of reads after that from 0 bp to 88 bp

ADD COMMENT
1
Entering edit mode

Galaxy can indeed do trimming from the 5'end. Just search more for the right tool.

ADD REPLY
0
Entering edit mode

This is not an answer. Please add this clarification to your original question and delete this "answer". It will help us all keep track of how your question is answered. Thanks.

ADD REPLY
0
Entering edit mode
11.3 years ago
Josh Herr 5.8k

If you're aligning your reads for phylogenetic analysis you'll have to have a matrix that is all the same length, so in that case, you'll have to trim. As Ido Tamir recommended, the FASTX TRIMMER from the Hannon Lab is very good.

ADD COMMENT

Login before adding your answer.

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