Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • ReorderSam.jar problem

    Hi

    I'm trying to use GATK RealignerTargetCreator but got an error message saying "Lexicographically sorted human genome sequence detected in reads."

    $ java -jar /local/software/rh53/GATK/Sting/dist/GenomeAnalysisTK.jar -T RealignerTargetCreator -R /temp/wjt/EXOME/OXC/Z1905/GATK/hg18.fa -I /temp/wjt/EXOME/OXC/Z1905/BAM/Z1905_novoalign.bam -o /temp/wjt/EXOME/OXC/Z1905/BAM/Z1905.intervals

    I've got a ref dictionary file (hg18.dict) and ran ReorderSam.jar which seemed to work ok (output below) but when rerunning GATK RealignerTargetCreator with the resorted bam I get the same error message.

    $ java -jar /local/software/rh53/GATK/Sting/dist/GenomeAnalysisTK.jar -T RealignerTargetCreator -R /temp/wjt/EXOME/OXC/Z1905/GATK/hg18.fa -I /temp/wjt/EXOME/OXC/Z1905/BAM/Z1905_novoalignRS.bam -o /temp/wjt/EXOME/OXC/Z1905/BAM/Z1905.intervals

    "Lexicographically sorted human genome sequence detected in reads."

    Thanks for any help regarding this query.

    Will

    ReorderSam.jar output:
    $ java -jar ReorderSam.jar INPUT=/temp/wjt/EXOME/OXC/Z1905/BAM/Z1905_novoalign.bam OUTPUT=/temp/wjt/EXOME/OXC/Z1905/BAM/Z1905_novoalignRS.bam REFERENCE=/temp/wjt/EXOME/OXC/Z1905/GATK/hg18.fa

    Executing as wjt@blue32 on Linux 2.6.18-164.6.1.el5 amd64; Java HotSpot(TM) 64-Bit Server VM 1.6.0_16-b01; Picard version: 1.64(1146)
    INFO 2012-03-13 15:14:41 ReorderSam SAM/BAM file
    INFO 2012-03-13 15:14:41 ReorderSam SN=%s LN=%d%nchr10135374737
    INFO 2012-03-13 15:14:41 ReorderSam SN=%s LN=%d%nchr10_random113275
    .
    .
    INFO 2012-03-13 15:14:41 ReorderSam Reordering SAM/BAM file:
    INFO 2012-03-13 15:14:41 ReorderSam Reordering read contig chr10 [index=0] to => ref contig chr10 [index=0]
    INFO 2012-03-13 15:14:41 ReorderSam Reordering read contig chr10_random [index=1] to => ref contig chr10_random [index=1]
    .
    .
    INFO 2012-03-13 15:14:41 ReorderSam Reordering read contig chrX_random [index=47] to => ref contig chrX_random [index=47]
    INFO 2012-03-13 15:14:41 ReorderSam Reordering read contig chrY [index=48] to => ref contig chrY [index=48]
    INFO 2012-03-13 15:14:41 ReorderSam Writing reads...
    INFO 2012-03-13 15:14:41 ReorderSam Processing chr10
    INFO 2012-03-13 15:15:33 ReorderSam Wrote 2877802 reads
    INFO 2012-03-13 15:15:33 ReorderSam Processing chr10_random
    .
    .
    INFO 2012-03-13 15:36:09 ReorderSam Wrote 6465436 reads
    [Tue Mar 13 15:36:09 GMT 2012] net.sf.picard.sam.ReorderSam done. Elapsed time: 21.48 minutes.
    Runtime.totalMemory()=1055588352

  • #2
    Hi

    I've got the same problem here

    When trying to use RealignerTargetCreator, GATK tells me that the sequence is in Lexicographically, so i tried to solve with ReorderSam which it works fine.
    But when I run RealignerTargetCreator on the sorted bam, i find the same error.

    Hope someone can help,
    Thanks!!

    Marco

    Comment


    • #3
      i also confront this problem, have you solve it?

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Strategies for Sequencing Challenging Samples
        by seqadmin


        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
        03-22-2024, 06:39 AM
      • seqadmin
        Techniques and Challenges in Conservation Genomics
        by seqadmin



        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

        Avian Conservation
        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
        03-08-2024, 10:41 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 06:37 PM
      0 responses
      8 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      8 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      49 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      66 views
      0 likes
      Last Post seqadmin  
      Working...
      X