밤파일 만들때 레퍼런스랑 지금 돌리는 레퍼런스랑 다를 경우 발생되는 에러이다.
<에러메세지>
org.broadinstitute.sting.utils.exceptions.ReviewedStingException: Invalid sequence number 24
at org.broadinstitute.sting.gatk.datasources.reads.GATKBAMIndex.<init>(GATKBAMIndex.java:100)
at org.broadinstitute.sting.gatk.datasources.reads.BAMSchedule.<init>(BAMSchedule.java:105)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.getNextOverlappingBAMScheduleEntry(BAMScheduler.java:205)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.advance(BAMScheduler.java:108)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.next(BAMScheduler.java:79)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.next(BAMScheduler.java:47)
at net.sf.picard.util.PeekableIterator.advance(PeekableIterator.java:71)
at net.sf.picard.util.PeekableIterator.next(PeekableIterator.java:57)
at org.broadinstitute.sting.gatk.datasources.reads.LowMemoryIntervalSharder.next(LowMemoryIntervalSharder.java:61)
at org.broadinstitute.sting.gatk.datasources.reads.LowMemoryIntervalSharder.next(LowMemoryIntervalSharder.java:36)
at org.broadinstitute.sting.gatk.datasources.reads.LocusShardStrategy.next(LocusShardStrategy.java:142)
at org.broadinstitute.sting.gatk.datasources.reads.LocusShardStrategy.next(LocusShardStrategy.java:42)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator$1.doBody(GenomeLocProcessingTracker.java:258)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator$1.doBody(GenomeLocProcessingTracker.java:250)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$WithLock.run(GenomeLocProcessingTracker.java:423)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator.next(GenomeLocProcessingTracker.java:250)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator.next(GenomeLocProcessingTracker.java:213)
at org.broadinstitute.sting.gatk.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:53)
at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:235)
at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:117)
at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:221)
at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:87)
==>
This is a terrible error message with what I suspect is a fairly simple cause: I think you're running with a different reference (or different version of the same reference) than the reference to which your BAMs were initially aligned. Can you check to make sure that the version of the reference you're using matches the BAMs exactly?
(http://getsatisfaction.com/gsa/topics/unifiedgenotyper_problem_with_nt_option)
<에러메세지>
org.broadinstitute.sting.utils.exceptions.ReviewedStingException: Invalid sequence number 24
at org.broadinstitute.sting.gatk.datasources.reads.GATKBAMIndex.<init>(GATKBAMIndex.java:100)
at org.broadinstitute.sting.gatk.datasources.reads.BAMSchedule.<init>(BAMSchedule.java:105)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.getNextOverlappingBAMScheduleEntry(BAMScheduler.java:205)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.advance(BAMScheduler.java:108)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.next(BAMScheduler.java:79)
at org.broadinstitute.sting.gatk.datasources.reads.BAMScheduler.next(BAMScheduler.java:47)
at net.sf.picard.util.PeekableIterator.advance(PeekableIterator.java:71)
at net.sf.picard.util.PeekableIterator.next(PeekableIterator.java:57)
at org.broadinstitute.sting.gatk.datasources.reads.LowMemoryIntervalSharder.next(LowMemoryIntervalSharder.java:61)
at org.broadinstitute.sting.gatk.datasources.reads.LowMemoryIntervalSharder.next(LowMemoryIntervalSharder.java:36)
at org.broadinstitute.sting.gatk.datasources.reads.LocusShardStrategy.next(LocusShardStrategy.java:142)
at org.broadinstitute.sting.gatk.datasources.reads.LocusShardStrategy.next(LocusShardStrategy.java:42)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator$1.doBody(GenomeLocProcessingTracker.java:258)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator$1.doBody(GenomeLocProcessingTracker.java:250)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$WithLock.run(GenomeLocProcessingTracker.java:423)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator.next(GenomeLocProcessingTracker.java:250)
at org.broadinstitute.sting.utils.threading.GenomeLocProcessingTracker$OwnershipIterator.next(GenomeLocProcessingTracker.java:213)
at org.broadinstitute.sting.gatk.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:53)
at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:235)
at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:117)
at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:221)
at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:87)
==>
This is a terrible error message with what I suspect is a fairly simple cause: I think you're running with a different reference (or different version of the same reference) than the reference to which your BAMs were initially aligned. Can you check to make sure that the version of the reference you're using matches the BAMs exactly?
(http://getsatisfaction.com/gsa/topics/unifiedgenotyper_problem_with_nt_option)
'Bioinformatics > Biological data analysis' 카테고리의 다른 글
[bwa] BWA 페어링에서 속도가 현저히 늦을 경우 (0) | 2012.02.20 |
---|---|
[hotspotter] hotspotter-1.2.1 설치 (0) | 2011.12.13 |
[base quality] Base Quality in BAM (0) | 2011.08.23 |
[java] 자바에서 시스템상의 명령어를 직접 호출할 경우 (0) | 2011.08.11 |
[picard] CollectGcBiasMetrixcs.jar (0) | 2011.08.10 |