SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Illumina HiSeq2500 data stormin Bioinformatics 3 08-27-2014 08:00 PM
Data Processing Pipeline question mboursnell Bioinformatics 10 02-28-2013 06:03 PM
Illumina Genome in a day on HiSeq2500 low99 Illumina/Solexa 4 02-16-2012 09:10 PM
Pipeline Illumina question bioinfogirl Bioinformatics 4 11-24-2008 04:00 AM

Reply
 
Thread Tools
Old 11-24-2016, 01:45 PM   #1
hamcan
Member
 
Location: Toronto

Join Date: Nov 2016
Posts: 19
Default Illumina HiSeq2500 Pipeline Question

Hi there,
I have just received my Illumina HiSeq data back. I have 2x250bp runs (so forward and reverse reads).

I ran a FASTQC and found that I have an adaptor sequence present in my data. Should I first remove the adaptors then merge my forward and reverse reads using FLASh? Or should I merge first, then remove the adaptors and trim at the same time (using trimmomatic).
Thank you!
hamcan is offline   Reply With Quote
Old 11-24-2016, 02:09 PM   #2
mastal
Senior Member
 
Location: uk

Join Date: Mar 2009
Posts: 667
Default

Remove the adapters first.

Note that for the reads that have adapters, the fragment length is less than 250 bp, and your forward and reverse reads, after adapter removal, should be reverse complements of each other.
mastal is offline   Reply With Quote
Old 11-24-2016, 02:12 PM   #3
hamcan
Member
 
Location: Toronto

Join Date: Nov 2016
Posts: 19
Default

Quote:
Originally Posted by mastal View Post
Remove the adapters first.

Note that for the reads that have adapters, the fragment length is less than 250 bp, and your forward and reverse reads, after adapter removal, should be reverse complements of each other.

So should I first use trimmomatic then to remove the adaptors and to overall trim the sequences at once? Then merge?

Thank you!
hamcan is offline   Reply With Quote
Old 11-24-2016, 02:48 PM   #4
Zapages
Member
 
Location: NJ

Join Date: Oct 2012
Posts: 94
Default

Merging the reads into one file depends on the task that you are working on.

and also do some quality control on the sequences.

Please FastQC to check the quality of after each trimming and/or quality control trimming.
Zapages is offline   Reply With Quote
Old 11-24-2016, 07:23 PM   #5
hamcan
Member
 
Location: Toronto

Join Date: Nov 2016
Posts: 19
Default

Quote:
Originally Posted by Zapages View Post
Merging the reads into one file depends on the task that you are working on.

and also do some quality control on the sequences.

Please FastQC to check the quality of after each trimming and/or quality control trimming.
Okay thanks!!
I am planning on doing cutadapt to remove the adaptor sequences. Then I will merge my sequences and then Trim them using trimmomatic so that I can blast them against the NCBI-nr database. This is an environmental metagenomics sample and I'm looking to see what bacterial species are present in certain environments.
hamcan is offline   Reply With Quote
Old 03-10-2017, 06:56 AM   #6
lwebs
Junior Member
 
Location: MA

Join Date: Mar 2017
Posts: 7
Default

I used FLASH v1.2.11 to merge forward and reverse reads from a HiSeq4000 run on environmental DNA (estimated frag length 280-700 bp, read length 150 nt) and am obtaining really low read overlap ~18.40% (counting both innie and outie pairs). Is this to be expected? Any help/insight would be appreciated!
lwebs is offline   Reply With Quote
Old 03-10-2017, 07:00 AM   #7
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,951
Default

Quote:
Originally Posted by lwebs View Post
I used FLASH v1.2.11 to merge forward and reverse reads from a HiSeq4000 run on environmental DNA (estimated frag length 280-700 bp, read length 150 nt) and am obtaining really low read overlap ~18.40% (counting both innie and outie pairs). Is this to be expected? Any help/insight would be appreciated!
Looking at your fragment lengths it sounds reasonable. Once your fragment size goes over ~ 300 bp R1/R2 will no longer overlap at 150 bp length.

You could try BBMerge from bbmap suite to see if that does any better in terms of getting better merging.
GenoMax is offline   Reply With Quote
Reply

Tags
flash, merge, pair end sequencing, trim adapter, trimmomatic

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 03:40 AM.


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