diffrence between number of reads in fastqc reports and collectalignmentsummarymetrics of Picard
hi everyone
i have 4 fastq files (two lane correspond to one organism), after filtration of reads, FastQC tell me that total number of reads is 628748302, i mapped them to reference genome and merged resulted two bam file and then sorted it. Now when i run collectalignmentsummarymetrics of Picard to get total and mapped read number, Picard tell me that total number of reads is 628601638. also when i used samtools flagstat, number of total reads counted 629770210, which number is correct?
thanks
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Samtools flagstat result
written by kiran •
-
Is there any difference in the output files when you make BAM files of a paired-end file separately…
written by soheil •
-
Understanding MarkDuplicates and featureCounts results
written by xosugarxspicexo •
-
samtools flagstat calculation
written by Nitha •
-
samtools stats result different from samtools flagstat
written by snishtala03 •
-
Percent reads mapped in Picard
written by phylofun
-
Counting paired-end sequencing mapped reads
written by heir_of_isildur88 •
-
samtools flagstat result is different from samtools view result
written by Louis Kok •
-
Number of mapped reads from BAM file
written by Prakki Rama
-
Using Ngs, What Is The Best Mirna Expression Normalization Method?
written by Doctoroots
See past threads on this question:
How come the QC-passed reads from samtools flagstat is different from FASTQ (R1+R2) read count?
Samtools flagstat number of reads do not match actual total number of reads.