SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 11-27-2017, 06:35 AM   #1
Razi
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, 06:41 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,578
Default

This question is cross-posted and has been answered: https://www.biostars.org/p/282328/#286154
GenoMax is offline   Reply With Quote
Old 12-01-2017, 11:52 PM   #3
[email protected]
Member
 
Location: Shenzhen, China

Join Date: Aug 2015
Posts: 15
Default

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: http://opengene.org/fastp/fastp.html
JSON report: http://opengene.org/fastp/fastp.json

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

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


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