minor modification in mapping_qc

wbaopaul · Mar 4, 2022 · b6b5566 · b6b5566
1 parent 2545e8e
commit b6b5566
Show file tree

Hide file tree

Showing 4 changed files with 14 additions and 8 deletions.
diff --git a/README.md b/README.md
@@ -155,7 +155,7 @@ Step by step guide to running scATAC-pro
 
 -   **IMPORTANT**: you can run scATAC-pro sequentially. The input of a later analysis module is the output of the previous analysis modules. The following tutorial uses fastq files downloaded from [PBMC10k 10X Genomics](https://support.10xgenomics.com/single-cell-atac/datasets/1.1.0/atac_v1_pbmc_10k?) 
 
--   <u>Run scATAC-pro sequentially (specifyi PEAK_CALLER = MACS2 and CELL_CALLER = FILTER or other values in the configure_user.txt file) </u>
+-   <u>Run scATAC-pro sequentially (specify PEAK_CALLER = MACS2 and CELL_CALLER = FILTER or other values in the configure_user.txt file) </u>
 
 ```
     $ scATAC-pro -s demplx_fastq 

diff --git a/scATAC-pro.wiki b/scATAC-pro.wiki
diff --git a/scripts/cell_mapping_qc.sh b/scripts/cell_mapping_qc.sh
@@ -46,9 +46,12 @@ tmp_bam_file=${output_dir}/tmp.bam
 
 
 if [[ $MAPPING_METHOD == bwa ]]; then
-    ${SAMTOOLS_PATH}/samtools view -@ $ncore -q 5 -f $flag0 -b ${input_pre}.bam > $tmp_bam_file
-   total_uniq_mapped=$( ${SAMTOOLS_PATH}/samtools view -c $tmp_bam_file )  ## number of unique mapped reads
-    rm $tmp_bam_file
+    #${SAMTOOLS_PATH}/samtools view -@ $ncore -q 5 -f $flag0 -b ${input_pre}.bam > $tmp_bam_file
+    #total_uniq_mapped=$( ${SAMTOOLS_PATH}/samtools view -c $tmp_bam_file )  ## number of unique mapped reads
+    #rm $tmp_bam_file
+
+    ## alternatively
+    total_uniq_mapped=$( ${SAMTOOLS_PATH}/samtools view -q 1 -@ $ncore -f $flag0 ${input_pre}.bam | grep -v XA: | wc -l )
 else
     ${SAMTOOLS_PATH}/samtools view -@ $ncore -q 5 -f $flag0 ${input_pre}.bam > $tmp_sam_file
    total_uniq_mapped=$( grep -E "@|NM:" $tmp_sam_file | grep -v "XS:" | wc -l )

diff --git a/scripts/mapping_qc.sh b/scripts/mapping_qc.sh
@@ -41,9 +41,12 @@ tmp_sam_file=${output_dir}/tmp.sam
 tmp_bam_file=${output_dir}/tmp.bam
 
 if [ $MAPPING_METHOD == 'bwa' ]; then
-    ${SAMTOOLS_PATH}/samtools view -@ $ncore -q 5 -f $flag0 -b ${input_pre}.positionsort.bam > $tmp_bam_file
-   total_uniq_mapped=$( ${SAMTOOLS_PATH}/samtools view -c $tmp_bam_file )  ## number of unique mapped reads
-    rm $tmp_bam_file
+    #${SAMTOOLS_PATH}/samtools view -@ $ncore -q 5 -f $flag0 -b ${input_pre}.positionsort.bam > $tmp_bam_file
+    #total_uniq_mapped=$( ${SAMTOOLS_PATH}/samtools view -c $tmp_bam_file )  ## number of unique mapped reads
+    #rm $tmp_bam_file
+
+    ## alternatively
+    total_uniq_mapped=$( ${SAMTOOLS_PATH}/samtools view -q 1 -@ $ncore -f $flag0 ${input_pre}.positionsort.bam | grep -v XA: | wc -l )
 else
     ${SAMTOOLS_PATH}/samtools view -@ $ncore -q 5 -f $flag0 ${input_pre}.positionsort.bam > $tmp_sam_file
    total_uniq_mapped=$( grep -E "@|NM:" $tmp_sam_file | grep -v "XS:" | wc -l )