SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Trimming left end (5') of reads?? blindtiger454 RNA Sequencing 23 01-11-2018 12:34 AM
Is left end (_1.fq) of normally mapped reads always positive strand? bioiion Illumina/Solexa 2 02-15-2012 07:19 AM
500 million reads needed for RNA-Seq?! epistatic RNA Sequencing 6 10-31-2011 03:53 PM
Tophat - Left Kept Reads jkozubek Bioinformatics 2 07-27-2011 05:40 AM
De novo assembler for 300 million Solexa reads odysseus Bioinformatics 10 11-05-2010 10:19 AM

Reply
 
Thread Tools
Old 08-04-2010, 04:26 AM   #1
Manu
Junior Member
 
Location: Freiburg, Germany

Join Date: May 2010
Posts: 4
Default What's left of 5 million reads

Hi,

I'm working with Illumina GAII data (5 million reads) from a 5' cap capture library. I did the mapping against a transcriptome with ELAND_standalone and now I'm left with 500000 reads that were mapped to a Gene-ID (after filtering). Most of the reads had more than one hit or had one or two mismatches. I used the default values for the filtering (not sure if I can change these with the standalone).

Now I just want to know if you get the same numbers (only 10 % mappable reads) with your data or if that number is too low.
Manu is offline   Reply With Quote
Old 08-04-2010, 06:24 AM   #2
NextGenSeq
Senior Member
 
Location: USA

Join Date: Apr 2009
Posts: 482
Default

10% is very low. BLAST or BLAT some of the unmapped reads to see if they are real or artifacts. If real obviously its a software processing error, if an artifact obviously a library problem.
NextGenSeq is offline   Reply With Quote
Old 08-05-2010, 07:48 AM   #3
mattanswers
Member
 
Location: Boston

Join Date: Oct 2009
Posts: 65
Default

Hi Manu,

I would try to download and use bowtie. It is much more flexible as to what bases you use to align as well as much faster. For me, it was a matter of hours for an alignment with ELAND as opposed to a couple of minutes using bowtie. So, with bowtie it is much easier to try different alignment parameters such as trimming 5' or 3' bases from your sequences. You can also print out a separate file of non-aligned sequences. It does not take too long to familiarize yourself with bowtie and in the long run seems well worth it.

Also, fastqc, http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/ can be used with s_*_sequence.txt files and will quickly provide you with some quality assessment. It will tell you about the base qualities, GC content, and also provide an idea if there were adapter seqs or many duplicate reads as well as other stats.
mattanswers is offline   Reply With Quote
Old 08-18-2010, 03:47 AM   #4
Manu
Junior Member
 
Location: Freiburg, Germany

Join Date: May 2010
Posts: 4
Default

I tried bowtie and bwa and both got better results. I think the ELAND filtering is very strict and I didn't find any information how to change the filtering parameter with the ELAND_standalone.
Manu is offline   Reply With Quote
Reply

Tags
eland, eland_standalone, mapping, mapping quality

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:31 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO