Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
paired-end adapter trimming vinay052003 Bioinformatics 16 05-02-2017 07:58 PM
Fixed length trimming of paired end reads for MISO sjuzhet Bioinformatics 1 03-06-2016 08:17 PM
Mapping paired-end reads, include also the SE-reads after trimming? JonB Bioinformatics 2 12-10-2015 02:03 AM
Aligning paired end read files that have different numbers of reads after trimming d1antho RNA Sequencing 2 01-21-2013 09:41 AM
paired-end reads mapped to genome.. gene with only one direction of paired-end reads? danwiththeplan Bioinformatics 2 09-22-2011 02:06 AM

Thread Tools
Old 11-27-2017, 05:35 AM   #1
Junior Member
Location: karaj

Join Date: Nov 2017
Posts: 1
Default Trimming the primers from paired end reads

Hi every one

I received some paired end reads containing specific Lexogen primers. So before assembly the primers must be removed from the reads. According to Lexogen protocol I must remove 9 nucleotides from 5' side of R1 and 6 nucleotides from 3' side of of R2. Is there any especial command in BBMap or any other software to accomplish this kind of trimming? In BBMap I just found command for removal specific number of nucleotides from tail of a read but for 5' side of the read I could not find a command as well. Would you please guide me to trim my reads?
Razi is offline   Reply With Quote
Old 11-27-2017, 05:41 AM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,655

This question is cross-posted and has been answered:
GenoMax is offline   Reply With Quote
Old 12-01-2017, 10:52 PM   #3
[email protected]
Location: Shenzhen, China

Join Date: Aug 2015
Posts: 15

I suggest fastp to do automatic adapter trimming, read filtering and quality control. fastp is ultra-fast since it's developed in C++ and with multi-threading support.

fastp has following features:
1, filter out bad reads (too low quality, too short, or too many N...)
2, cut low quality bases for per read in its 5' and 3' by evaluating the mean quality from a sliding window (like Trimmomatic but faster).
3, trim all reads in front and tail
4, cut adapters. Adapter sequences can be automatically detected,which means you don't have to input the adapter sequences to trim them.
5, correct mismatched base pairs in overlapped regions of paired end reads, if one base is with high quality while the other is with ultra low quality
6, preprocess unique molecular identifer (UMI) enabled data, shift UMI to sequence name.
7, report JSON format result for further interpreting.
8, visualize quality control and filtering results on a single HTML page (like FASTQC but faster and more informative).
9, split the output to multiple files (0001.R1.gz, 0002.R1.gz...) to support parallel processing. Two modes can be used, limiting the total split file number, or limitting the lines of each split file.
10, support long reads (data from PacBio / Nanopore devices).

fastp creates reports in both HTML and JSON format.

HTML report:
JSON report:

fastp is an open source project at github:
OpenGene(Libraries and tools for NGS data analysis),AfterQC(Fastq Filtering and QC)
FusionDirect.jl( Detect gene fusion), SeqMaker.jl(Next Generation Sequencing simulation) is offline   Reply With Quote

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:32 AM.

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