Go Back   SEQanswers > Bioinformatics > Bioinformatics

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

Thread Tools
Old 04-07-2013, 08:05 PM   #1
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

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 04:24 AM.

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