SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BBDUK and BBMERGE, which goes first chloe1005 Bioinformatics 1 04-12-2018 03:43 AM
NexteraXT Bionalayzer - undertagmentation? Meyana Illumina/Solexa 15 12-14-2017 07:23 PM
bbduk.sh barcode filter lamon Bioinformatics 8 11-18-2017 03:51 AM
Removing contamination with BBDUK glfrey Bioinformatics 2 10-02-2017 04:08 AM
Trimmomatic vs bbduk.sh mslider Bioinformatics 1 04-18-2017 10:10 AM

Reply
 
Thread Tools
Old 05-07-2018, 10:07 AM   #1
arivers
Junior Member
 
Location: Gainesville, FL

Join Date: May 2018
Posts: 2
Default BBduk for NexteraXT

I'm using BBduk 38.00 for adaptor trimming on NexteraXT libraries multiplexed with i5 and i7 indices and sequenced on an NextSeq500 2x150:

Code:
bbduk.sh in=sample.fq.gz ref=adapters.fa bhist=bhist.txt k=23 hdist=1 ktrim='r' overwrite=t mink=8 tbe tbo ftm=5
When I plot the bhist output from BBduk (A, C, G, T, N % at each position on read one then read2), I get strong GC and AT precent divergence over the first 16-18 bp of read 1 and read 2.

It looks like I'm missing some forward primer or transposase sequence and, in fact, the most common bases are similar to the read 2 Nextera transposase sequence.

Here is an example plot:
https://photos.app.goo.gl/hH2UjJ9iErPbC5nP9
arivers is offline   Reply With Quote
Old 05-07-2018, 10:22 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,731
Default

Have you come across this blog post?

What you are observing is the "bias" observed similar to random primed libraries that is present in nextera libraries.
GenoMax is offline   Reply With Quote
Old 05-07-2018, 10:43 AM   #3
arivers
Junior Member
 
Location: Gainesville, FL

Join Date: May 2018
Posts: 2
Default

I knew that was an issue for random hexamer primed RNAseq libraries and I've heard that tagmentation sites were non-random. I usually work with Truseq/kapa data so this is my first time actually seeing how non-random things are.

Using BBduk to left trim with a hamming distance of 2 or 3 reduces, or at least masks the positional sequence bias.

Glad to know the consensus is that it's not a problem.
arivers is offline   Reply With Quote
Old 05-07-2018, 10:52 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,731
Default

You should align the data without doing any additional manipulations beyond removal of adapters. Data should align normally.
GenoMax is offline   Reply With Quote
Reply

Tags
bbduk, gc divergence, nexteraxt

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 06:59 PM.


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