SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Detect large deletions srr8623 Bioinformatics 2 11-28-2012 02:18 PM
gatk countcovariate ERROR MESSAGE Read starting with deletion. Cigar: 1D22M1I71M8S aituka Bioinformatics 2 09-13-2012 05:07 AM
Illumina insertions deletions andreitudor Bioinformatics 2 02-14-2012 06:20 AM
Coverage and deletions Hena Bioinformatics 0 01-05-2011 10:05 AM
SOAP2 and deletions/insertions joa_ds Bioinformatics 7 09-21-2010 03:59 AM

Reply
 
Thread Tools
Old 04-07-2013, 07:05 PM   #1
JezSupreme
Junior Member
 
Location: Perth, Western Australia

Join Date: Mar 2013
Posts: 6
Default GATK -rf BadCigar (Reads starting with deletions)

Please excuse any naive or ignorant questions here, but I'm still very new to bioinformatics. I'm trying to run some illumina sequence data through the GATK pipeline and have been running into problems with reads beginning with deletions. I've been running some initial commands on my .bam files to make sure GATK is OK with the files and while I can run CountReads successfully, I get a BadCigar error return when I run the CountLoci command. Of the 3 data files I'm looking at, all three return BadCigar errors, and all 3 errors are related to reads starting with deletions.

While I can run these files successfully by using the '-rf BadCigar' flag, this use raises a couple of questions:

1.) Why does GATK consider these reads to be in error when we have yet to map the InDels and conduct a realignment?

2.) The process of using the BadCigar flag seems to ignore the entire read, which is not ideal for my purposes. Is it possible to force the analysis to run and simply permit the bad cigars rather than ignore the entire read?

3.)The GATK doesn't seem to report which reads have been ignored (beyond the initial run of the file not using the flag, in which it error's out at the first bad cigar). Is it possible to have the GATK report each bad cigar as it is ignored in the analysis?

Thanks to anyone who can shed some light on these questions.
JezSupreme is offline   Reply With Quote
Reply

Tags
badcigar, bioinformatics, gatk

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:18 PM.


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