Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • GATK IndelRealigner error

    I have problems with the GATK IndelRealigner. I run GATK version 1.0.5365

    Except the standard input and output files I use:
    -compress 0
    -L chr_number
    -B:indels,VCF dbSNP_indelfile.vcf

    The pipeline I use look like this:
    1. Split fastq-files into smaller parts (paired-end)
    2. Align each parts with BWA
    3. sort each part with samtools
    4. merge with samtools
    5. run GATK -T IndelRealigner on each chromosome (after RealignerTargetCreator of course)

    All chromosomes except number 11 is going ok.

    In chr11 I get this error:
    java.lang.IllegalArgumentException: Alignments added out of order in SAMFileWriterImpl.addAlignment for /glob/jonasc/alignment_test2/pipe_line_test_sorted_realign.bwa.11.bam. Sort order is coordinate. Offending records are at [11:51571970] and [11:51571888]

    Using samtools view on the region in the input bam-file I can see that the reads are sorted.

    Any help to solve this would be useful.

  • #2
    Get some satisfaction: http://getsatisfaction.com/gsa

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin




      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
      04-22-2024, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Today, 11:49 AM
    0 responses
    12 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 08:47 AM
    0 responses
    16 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    61 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    60 views
    0 likes
    Last Post seqadmin  
    Working...
    X