epigen-UCSD / epigen_ucsd_django

1 stars 1 forks source link

[bug] shallow seq's saturation is not correct. #503

Closed biomystery closed 3 years ago

biomystery commented 3 years ago

{'estimated_number_of_cells': 7058, 'mean_reads_per_cell': 4431, 'median_genes_per_cell': 827, 'number_of_reads': 31279107, 'valid_barcodes': 0.968, 'sequencing_saturation': 0.9, 'q30_bases_in_barcode': 0.979, 'q30_bases_in_rna_read': 0.932, 'q30_bases_in_sample_index': None, 'q30_bases_in_UMI': 0.974, 'reads_mapped_to_genome': 0.948, 'reads_mapped_confidently_to_genome': 0.932, 'reads_mapped_confidently_to_intergenic_regions': 0.63, 'reads_mapped_confidently_to_intronic_regions': 0.0, 'reads_mapped_confidently_to_exonic_regions': 0.869, 'reads_mapped_confidently_to_transcriptome': 0.538, 'reads_mapped_antisense_to_gene': 0.306, 'frac_reads_in_cells': 0.817, 'total_genes_detected': 21329, 'median_UMI_counts_per_cell': 1469}

(django) [zhc268@epigenomics epigen_ucsd_django]$ cat ~/data/outputs/scRNA/JB_262/outs/metrics_summary.csv Estimated Number of Cells,Mean Reads per Cell,Median Genes per Cell,Number of Reads,Valid Barcodes,Sequencing Saturation,Q30 Bases in Barcode,Q30 Bases in RNA Read,Q30 Bases in UMI,Reads Mapped to Genome,Reads Mapped Confidently to Genome,Reads Mapped Confidently to Intergenic Regions,Reads Mapped Confidently to Intronic Regions,Reads Mapped Confidently to Exonic Regions,Reads Mapped Confidently to Transcriptome,Reads Mapped Antisense to Gene,Fraction Reads in Cells,Total Genes Detected,Median UMI Counts per Cell "7,058","4,431",827,"31,279,107",96.8%,9.0%,97.9%,93.2%,97.4%,94.8%,93.2%,6.3%,0.0%,86.9%,53.8%,30.6%,81.7%,"21,329","1,469"

Updated

{'estimated_number_of_cells': 7058, 'mean_reads_per_cell': 4431, 'median_genes_per_cell': 827, 'number_of_reads': 31279107, 'valid_barcodes': 0.968, 'sequencing_saturation': 0.09, 'q30_bases_in_barcode': 0.9790000000000001, 'q30_bases_in_rna_read': 0.932, 'q30_bases_in_sample_index': None, 'q30_bases_in_UMI': 0.9740000000000001, 'reads_mapped_to_genome': 0.948, 'reads_mapped_confidently_to_genome': 0.932, 'reads_mapped_confidently_to_intergenic_regions': 0.063, 'reads_mapped_confidently_to_intronic_regions': 0.0, 'reads_mapped_confidently_to_exonic_regions': 0.8690000000000001, 'reads_mapped_confidently_to_transcriptome': 0.5379999999999999, 'reads_mapped_antisense_to_gene': 0.306, 'frac_reads_in_cells': 0.8170000000000001, 'total_genes_detected': 21329, 'median_UMI_counts_per_cell': 1469}

biomystery commented 3 years ago

correct exampe;

(django) [zhc268@epigenomics epigen_ucsd_django]$ cat ~/data/outputs/scRNA/AVD_580/outs/metrics_summary.csv Estimated Number of Cells,Mean Reads per Cell,Median Genes per Cell,Number of Reads,Valid Barcodes,Sequencing Saturation,Q30 Bases in Barcode,Q30 Bases in RNA Read,Q30 Bases in Sample Index,Q30 Bases in UMI,Reads Mapped to Genome,Reads Mapped Confidently to Genome,Reads Mapped Confidently to Intergenic Regions,Reads Mapped Confidently to Intronic Regions,Reads Mapped Confidently to Exonic Regions,Reads Mapped Confidently to Transcriptome,Reads Mapped Antisense to Gene,Fraction Reads in Cells,Total Genes Detected,Median UMI Counts per Cell "5,658","15,998",735,"90,519,258",88.0%,10.0%,89.8%,74.6%,87.2%,87.3%,83.0%,75.9%,7.3%,7.0%,61.6%,53.7%,3.8%,28.0%,"18,958","1,230"

biomystery commented 3 years ago

6aaa6f8