Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
All the CIGAR strings of reads mapped by Bowtie are "I" 610617109 Bioinformatics 1 11-18-2015 12:42 AM
Trim off variable length 'N' strings at the end of the read illinu Bioinformatics 8 05-16-2014 07:41 AM
Accessing SAM CIGAR strings in ShortRead (Bioconductor) mixter Bioinformatics 1 11-26-2012 08:34 AM
2NA 4NA search strings HMorrison Bioinformatics 0 12-13-2011 12:41 PM
CIGAR strings and 'N' symbols pieffe Bioinformatics 4 03-11-2010 12:08 PM

Thread Tools
Old 06-23-2016, 11:31 PM   #1
Location: Germany

Join Date: Sep 2014
Posts: 12
Default consensus of CIGAR strings

Hi guys,

I am writing an tool for amplicon sequencing, which needs to condense reads mapping to the same position. At the moment, I just take the most frequent read, which is obviously far from perfect, especially in low quality or low coverage situations.
While it is easy to calculate the consensus of the base string and to approximate the quality, the CIGAR string is much harder. Does anyone has an algorithm or even source code to do this?


Flo1000 is offline   Reply With Quote
Old 06-24-2016, 09:24 AM   #2
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

One solution is to do multiple sequence alignment... I'm sure there are tools for this exact problem, but I don't know what they are, specifically. I think the Mothur package has something for that, though, specifically for creating a consensus of reads from a given amplicon.

Another approach would be to gather all the reads mapped to the same position, and assemble them. This might be computationally more efficient than multiple sequence alignment. Possibly, an OLC assembler like MIRA would be best.
Brian Bushnell 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 03:16 AM.

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