SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Quantitative Linkage Analysis 454 Amplicon Sequencing sdl11 Bioinformatics 0 01-26-2012 01:55 PM
small rna sequencing sample quality control entomology RNA Sequencing 2 01-16-2012 03:59 PM
RNA-Seq: htSeqTools: High-Throughput Sequencing Quality Control, Processing and Visua Newsbot! Literature Watch 0 12-27-2011 02:00 AM
The quality control index or parameter of the whole genome de novo sequencing wave001 Bioinformatics 0 06-25-2010 08:28 AM
Volunteers wanted! Sequencing Quality Control Project (SEQC) Joann Events / Conferences 2 10-09-2009 03:24 AM

Reply
 
Thread Tools
Old 02-21-2012, 06:32 AM   #1
Mucki0815
Junior Member
 
Location: SouthernGermany

Join Date: Sep 2011
Posts: 2
Question Amplicon Sequencing Quality control and downstream analysis

Hello everyone,

I'm posting this thread because I'm having a hard time trying to do some quality control on our sequencing data.

We are doing amplicon sequencing on an Ion PGM using a multiplex PCR for amplicon generation. I tried using FastQC for a fast and easy assassment of run quality but it doesn't seem to be suitable for this kind of application.

The program tells me that more than 80% of my reads are duplicates which in my opinion is an artifact created by the fact that only the ends of the reads are compared to each other. These ends are characterised by the PCR primer regions used in the multiplex PCR and are therefore (largely) identical. Are there any tools for duplicate detection that ignore the outermost base pairs?

As our reads have a length between 30 and 200 bp (mean 72) I am currently using bwasw for alignment. Base quality values are pretty good up to bp 100 but than they start to slope down rapidly. As I wasn't able to find this somewhere else I would like to ask here if anybody knows what happens to these low quality bases at the end of the reads in the bwasw algorythm? In the standard bwa algorythm a -q trim command can be used but there is no such thing for bwasw, does anybody know why?`

Thanks beforehand for any help provided.

Greatings,
Mucki
Mucki0815 is offline   Reply With Quote
Old 10-03-2012, 11:05 PM   #2
jjmmii
Junior Member
 
Location: Hong Kong

Join Date: Oct 2012
Posts: 4
Default

I would also like to know why as well. Any help would be appreciated.

Jamie =]
jjmmii 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 10:42 PM.


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