mauranolab / mapping

2 stars 1 forks source link

Normalize reads reported in counts.txt files to reads per 10mm #219

Closed cadley-nyulangone closed 4 years ago

cadley-nyulangone commented 4 years ago
new reads reported = int(10,000,000 * (number of reads in the range) / (total reads from the sample) + 0.5)
                                       ============================
                                      This is what we currently report

total reads is computed this way:

    num_bam1_reads=$(samtools view -c -F 512 ${bam1})
cadley-nyulangone commented 4 years ago

fixes #217

cadley-nyulangone commented 4 years ago

It's set in line 65 of merge_bamintersect.sh

It's passed into counts_table.sh in line 33 there.

mattmaurano commented 4 years ago

great

you need to update help message in submit_bamintersect.sh

mattmaurano commented 4 years ago

I fixed this