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
                        Advanced Tools Transforming the Field of Cytogenomics
                        by seqadmin


                        At the intersection of cytogenetics and genomics lies the exciting field of cytogenomics. It focuses on studying chromosomes at a molecular scale, involving techniques that analyze either the whole genome or particular DNA sequences to examine variations in structure and behavior at the chromosomal or subchromosomal level. By integrating cytogenetic techniques with genomic analysis, researchers can effectively investigate chromosomal abnormalities related to diseases, particularly...
                        09-26-2023, 06:26 AM
                      • seqadmin
                        How RNA-Seq is Transforming Cancer Studies
                        by seqadmin



                        Cancer research has been transformed through numerous molecular techniques, with RNA sequencing (RNA-seq) playing a crucial role in understanding the complexity of the disease. Maša Ivin, Ph.D., Scientific Writer at Lexogen, and Yvonne Goepel Ph.D., Product Manager at Lexogen, remarked that “The high-throughput nature of RNA-seq allows for rapid profiling and deep exploration of the transcriptome.” They emphasized its indispensable role in cancer research, aiding in biomarker...
                        09-07-2023, 11:15 PM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, Yesterday, 09:36 AM
                      0 responses
                      8 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 10-02-2023, 07:14 AM
                      0 responses
                      15 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 09-29-2023, 09:38 AM
                      0 responses
                      16 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 09-27-2023, 06:57 AM
                      0 responses
                      16 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X