Summary the "total read count" and "read length" of each samples within
a BarcodeQcSet
object, and output a data.frame
with sample by
row and different metrics by column.
bc_summary_seqQc(x)
# S4 method for class 'BarcodeQcSet'
bc_summary_seqQc(x)
A data.frame
with 5 columns: sample_name
,
total_read
, median_read_length
, p5_read_length
and
p95_read_length.
fq_file <- dir(
system.file("extdata", "mef_test_data", package = "CellBarcode"),
full=TRUE)
bc_summary_seqQc(bc_seq_qc(fq_file))
#> sample_name total_read
#> <char> <num>
#> 1: 5290_10_BCM_195_mef_mixb_GTCATTG_S11_R1_001.fastq.gz 2000
#> 2: 5290_1_BCM_50000_mef_mixa_GTTCTCC_S2_R1_001.fastq.gz 2000
#> 3: 5290_2_BCM_50000_mef_mixb_GATGTGT_S5_R1_001.fastq.gz 2000
#> 4: 5290_3_BCM_12500_mef_mixa_TGCCTTG_S4_R1_001.fastq.gz 2000
#> 5: 5290_4_BCM_12500_mef_mixb_TAACTGC_S8_R1_001.fastq.gz 2000
#> 6: 5290_5_BCM_3125_mef_mixa_GCTTCCA_S9_R1_001.fastq.gz 2000
#> 7: 5290_6_BCM_3125_mef_mixb_TGTGAGT_S7_R1_001.fastq.gz 2000
#> 8: 5290_7_BCM_781_mef_mixa_CCTTACC_S12_R1_001.fastq.gz 2000
#> 9: 5290_8_BCM_781_mef_mixb_CGTATCC_S13_R1_001.fastq.gz 2000
#> 10: 5290_9_BCM_195_mef_mixa_GTACTGT_S14_R1_001.fastq.gz 2000
#> p5_read_length median_read_length p95_read_length
#> <num> <num> <num>
#> 1: 100 100 100
#> 2: 100 100 100
#> 3: 100 100 100
#> 4: 100 100 100
#> 5: 100 100 100
#> 6: 100 100 100
#> 7: 100 100 100
#> 8: 100 100 100
#> 9: 100 100 100
#> 10: 100 100 100
###