SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
quality problem when mapping with gMap lixiangru Bioinformatics 1 11-16-2011 06:40 PM
Problem with BWA mapping second read pbluescript Bioinformatics 3 10-12-2011 07:52 AM
Bowtie mapping problem harrike Bioinformatics 1 05-09-2011 06:29 AM
sinlge-end mapping problem cliff Bioinformatics 0 05-18-2010 10:02 AM
SNP mapping quality problem ptong7 Bioinformatics 0 12-23-2009 12:22 PM

Reply
 
Thread Tools
Old 10-05-2011, 11:33 PM   #1
StefanP
Member
 
Location: New Zealand

Join Date: Oct 2011
Posts: 10
Default deletion mapping problem

Hello,

I'm mapping 454 reads against a reference genome using bwa and then create sorted bam files using Samtools.

However, most my samples do have only one copy of a 9bp motif, whereas the reference has two. In some reads of my sequence data the motif will be followed by gaps and in some the gaps will be followed by the motif and in the end the consensus will call it as a tandem repeat.

I tried to use GATK for realignment, but it didn't work. Has anyone used GATK for this purpose successfully? Or does anyone have other ideas how to solve the problem? others than de-novo assembly?

I'm building up a fully automized pipeline so I use command line based software and scripts.

Help would be highly appreciated!

Thanks,
Stefan
StefanP is offline   Reply With Quote
Old 10-06-2011, 12:37 AM   #2
Bukowski
Senior Member
 
Location: UK

Join Date: Jan 2010
Posts: 390
Default

From : http://www.broadinstitute.org/gsa/ga...etCreator.html

"Another important note: because reads produced from the 454 technology inherently contain false indels, the realigner will not currently work with them (or with reads from similar technologies)."
Bukowski is offline   Reply With Quote
Old 10-06-2011, 12:42 AM   #3
StefanP
Member
 
Location: New Zealand

Join Date: Oct 2011
Posts: 10
Default

Thanks. But do you have another idea?

Cheers,
Stefan
StefanP is offline   Reply With Quote
Old 10-06-2011, 03:20 AM   #4
BaCh
Member
 
Location: Germany

Join Date: May 2008
Posts: 79
Default

Quote:
Originally Posted by StefanP View Post
Hello,
Help would be highly appreciated!
If the genome you are mapping to is a bacterium or lower eukaryote (30 to 100 MB), give MIRA a try. I wrote the mapping routines to behave consistently for exactly these cases.

B.
BaCh 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 03:00 PM.


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