SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Pacific Biosciences



Similar Threads
Thread Thread Starter Forum Replies Last Post
P6C4 has higher error rate than P5C3? ymc Pacific Biosciences 4 10-07-2015 05:30 PM
454 Roche. homopolymer error rubyb 454 Pyrosequencing 1 06-17-2015 07:54 AM
Homopolymer run errors, polyA bias sulfobus 454 Pyrosequencing 6 04-09-2010 04:15 AM
454 homopolymer error rate joa_ds 454 Pyrosequencing 10 12-18-2009 10:33 AM
Homopolymer error in 454 sequencing data Salmon 454 Pyrosequencing 3 02-06-2009 09:08 AM

Reply
 
Thread Tools
Old 04-25-2017, 05:24 AM   #1
cklopp
Member
 
Location: Toulouse France

Join Date: Sep 2009
Posts: 12
Default PacBio P6C4 homopolymer error bias

Using the alignment of the raw pacbio reads on a reference genome it is possible to count the number of insertions and deletions at each position (from the mpileup file) and to sum these figures for homopolymers. It is then possible to calculate the mean error rates per variation type and homopolymer length.

The results for one of our genomes of interest is presented hereunder.
We clearly see that the insertion proportion decreases with the homopolymer length and on the contrary the deletion proportion increases.

http://genoweb.toulouse.inra.fr/~klo...ion_pacbio.png

NB1 : this genome having only few homopolymers larger than 7 base pairs the corresponding figures have to be taken with caution.

NB2 : We have the same shapes for other genomes.

PacBio data is meant to have random errors which does not seem to be true looking at these figures.

I'm interested in any finding or comment about this issue.
cklopp is offline   Reply With Quote
Old 04-26-2017, 12:21 AM   #2
dcopetti
Junior Member
 
Location: Arizona

Join Date: Nov 2011
Posts: 2
Default

interesting.
Does the bias disappear after Pilon polishing with Illumina reads, for example?
It would be also nice to see a similar analysis about Nanopore data.
good point, we need to be aware of that
dcopetti is offline   Reply With Quote
Old 04-26-2017, 12:40 AM   #3
cklopp
Member
 
Location: Toulouse France

Join Date: Sep 2009
Posts: 12
Default

Unfortunately I do not have Illumina reads for this genome, I will therefore not be able to answer your pilon related request.
For Oxford Nanopore we are aware of the evident homopolymer bias.
cklopp 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 11:20 PM.


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