Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Use VALIDATION_STRINGENCY=LENIANT

    That way, it will report all the errors it sees to STDIN, but it will finish the job anyway.

    Comment


    • #17
      Originally posted by swbarnes2 View Post
      Use VALIDATION_STRINGENCY=LENIANT

      That way, it will report all the errors it sees to STDIN, but it will finish the job anyway.
      That's true, it doesn't complain, but the size of output file is 0 !
      Is it something, the sequencing vendor should fix at their end ?

      Comment


      • #18
        any hint ? I have asked this question in picard, but havent heard from them either;
        I appreciate any hint,

        Comment


        • #19
          does spelling LENIENT correctly help?

          Comment


          • #20
            Originally posted by gringer View Post
            does spelling LENIENT correctly help?
            My spelling in the command line is already correct,
            Here, I just Quoted what "swbarnes2" has suggested ...

            Comment


            • #21
              Originally posted by Bruins View Post
              Well... the problem is that after reading this thread and some googling, I get the feeling that SILENT (or IGNORE=) is the only solution, other than running FixMates.

              That takes me back to my question one: what would be the consequence of using SILENT? Drio would you like to comment on that some more?

              I'm planning to run some tests, I'll report back later.

              Is anyone still interested in this question? I recently noticed the consequence of using the lenient stringency in Picard. When AddOrReplace read groups, the strict mode will generate some column like "RG:Z:***" in the reads, while the lenient mode will not, although both modes will generate bam files with @RG in the head. This will cause a further problem when calling SNPs in GATK, which just assume that you did not add read groups, while reading groups is required by GATK.

              I spotted the problem, but havenot figured out how to solve it. Could anyone help me here?

              Comment


              • #22
                I am having this problem now with sam created by bwa-0.6.2

                java -Xmx5g -jar ../picard/picard-tools-1.39/SortSam.jar I=P6T.sam SO=coordinate O=P6Tf.bam VALIDATION_STRINGENCY=LENIENT CREATE_INDEX=true

                Ignoring SAM validation error due to lenient parsing:
                Error parsing text SAM file. MAPQ should be 0 for unmapped read.; Line 10820584
                Ignoring SAM validation error due to lenient parsing:
                Error parsing text SAM file. CIGAR should have zero elements for unmapped read.;

                I am getting about 30 of these errors among my 39mil reads. Presumably because of the LENIENT flag, the program finishes with a decent sized bam. Is it ok?

                Is this fixed in the bwa-0.7.5a and picard-1.94 combo??

                Comment


                • #23
                  Interesting. With bwa-0.7.5a and picard-1.94 combo, CIGAR errors are gone but MAPQ errors remain.

                  Comment


                  • #24
                    Hey guys,

                    Truly use VALIDATION_STRINGENCY=LENIENT OR SILENT. when u run most of the picard commands.

                    That will definitely help. I am experiencing it personally.

                    Good luck ahead to all NGS data analysing friends.

                    Regards,
                    Vishnu.

                    Comment


                    • #25
                      Yes, it complains less ! but because we shut it down ; It does not necessarily solve the issue. So I would like to add, watch out the output when you are walking with VALIDATION_STRINGENCY=LENIENT OR SILENT

                      Comment

                      Latest Articles

                      Collapse

                      • seqadmin
                        Recent Innovations in Spatial Biology
                        by seqadmin


                        Spatial biology is an exciting field that encompasses a wide range of techniques and technologies aimed at mapping the organization and interactions of various biomolecules in their native environments. As this area of research progresses, new tools and methodologies are being introduced, accompanied by efforts to establish benchmarking standards and drive technological innovation.

                        3D Genomics
                        While spatial biology often involves studying proteins and RNAs in their...
                        Yesterday, 07:30 PM
                      • seqadmin
                        Advancing Precision Medicine for Rare Diseases in Children
                        by seqadmin




                        Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
                        12-16-2024, 07:57 AM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, 12-30-2024, 01:35 PM
                      0 responses
                      21 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 12-17-2024, 10:28 AM
                      0 responses
                      41 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 12-13-2024, 08:24 AM
                      0 responses
                      55 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 12-12-2024, 07:41 AM
                      0 responses
                      40 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X