Dear all,
would you please advise: beside SAMTOOLS, PICARD or SAMBAMBA, is there any tool that you would recommend in order to sort a BAM based on read name ? It is a BAM file from EGA that contains cancer sequencing data, and when I do run PICARD SortSam, I am getting the following error (below), and the file does not get sorted. Thanks !
Exception in thread "main" htsjdk.samtools.SAMFormatException: SAM validation error: ERROR: Record 957453876, Read name HWI-ST7001002R:223:C14GPACXX:3:1305:7471:56486, MAPQ should be 0 for unmapped read
•
written
by
Bogdan
1 answer
That seems to be a very strange read, how can it have a MAPQ > 0?
Anyway, you can either try samtools sort -n to sort by name, or try setting PICARD's VALIDATION_STRINGENCY to 'lenient', so for example
java -jar picard.jar SortSam I=your_input O=your_sorted_output VALIDATION_STRINGENCY=LENIENT
By default VALIDATION_STRINGENCY is set to STRICT, which is probably the reason why it stops.
Log in to answer this question.
-
SortSam Error after BWA alignment
written by loy_loy •
-
Picard CollectSequencingArtifactMetrics error
written by jrleary
-
Picard PCR duplicate gives an error
written by Learner
-
No real operator (M|I|D|N) in CIGAR
written by boymin2020 •
-
SAM validation error
written by fire_water
-
Mark duplicates fails using samtools view to remove unmapped reads
written by jsgounot
-
Errors in sam file after bwa SAMPE alignment.
written by kirannbishwa01
-
FixMateInformation issue in the process of running picard
written by zengtony743
-
Problem With Picard Dedup In Illumina Alignment Pipeline
written by tommivat
-
Gatk Realignertargetcreator Sam/Bam Header Error
written by Houkto •
