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
                        Choosing Between NGS and qPCR
                        by seqadmin



                        Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
                        10-18-2024, 07:11 AM
                      • seqadmin
                        Non-Coding RNA Research and Technologies
                        by seqadmin




                        Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

                        Nobel Prize for MicroRNA Discovery
                        This week,...
                        10-07-2024, 08:07 AM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, 11-01-2024, 06:09 AM
                      0 responses
                      18 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 10-30-2024, 05:31 AM
                      0 responses
                      18 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 10-24-2024, 06:58 AM
                      0 responses
                      24 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 10-23-2024, 08:43 AM
                      0 responses
                      53 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X