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
                        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
                      67 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X