![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to use Picard's MarkDuplicates | cliff | Bioinformatics | 12 | 01-26-2015 11:56 PM |
About Picard | Jin Pan | Bioinformatics | 8 | 07-15-2014 05:55 AM |
picard error | seq_GA | Bioinformatics | 2 | 07-15-2014 05:50 AM |
Picard and GATK | poisson200 | Bioinformatics | 6 | 08-15-2012 08:04 AM |
Picard: HsMetrics | Seq84 | Bioinformatics | 0 | 04-18-2011 03:24 AM |
![]() |
|
Thread Tools |
![]() |
#1 |
Member
Location: London Join Date: Jan 2010
Posts: 65
|
![]()
Hello all,
I merged two ELAND export files by using samtools, but samtools flagstat shows that 191134982 in total 0 QC failure 0 duplicates 139825529 mapped (73.16%) 191134982 paired in sequencing 95567491 read1 95567491 read2 0 properly paired (0.00%) 115989136 with itself and mate mapped 23836393 singletons (12.47%) 0 with mate mapped to a different chr 0 with mate mapped to a different chr (mapQ>=5) Why it shows 0 properly paired? Couldn't run picard MarkDuplicates, thread exception: Exception in thread "main" java.lang.NullPointerException at net.sf.picard.sam.DiskReadEndsMap.remove(DiskReadEndsMap.java:87) at net.sf.picard.sam.MarkDuplicates.buildSortedReadEndLists(MarkDuplicates.java:289) at net.sf.picard.sam.MarkDuplicates.doWork(MarkDuplicates.java:126) at net.sf.picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:143) at net.sf.picard.sam.MarkDuplicates.main(MarkDuplicates.java:105) Any help? Many thanks. |
![]() |
![]() |
![]() |
#2 | ||
Nils Homer
Location: Boston, MA, USA Join Date: Nov 2008
Posts: 1,285
|
![]() Quote:
Quote:
Send both to samtools-help@lists.sourceforge.net |
||
![]() |
![]() |
![]() |
#3 |
Member
Location: GZ,China Join Date: Jun 2014
Posts: 42
|
![]()
I got the same error. Have you solved your problem, bair or nilshomer?
|
![]() |
![]() |
![]() |
Thread Tools | |
|
|