Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

What is the difference between rmdup_sam and rmdup_sam2 files? #3

Open
chatterjee89 opened this issue Feb 5, 2021 · 0 comments
Open

Comments

@chatterjee89
Copy link

chatterjee89 commented Feb 5, 2021

Hi Dr. Cao,

I had some questions regarding the duplicate removal steps. With the following code, I understand that we are removing duplicates based on identical UMI, tagmentation site, reverse transcription index and ligation adaptor index...

`echo "Start removing duplicates..."
echo input folder: $filtered_sam_folder
echo output folder: $rmdup_sam_folder
mkdir -p $rmdup_sam_folder
module unload python

bash_script=$script_path/sci3_rmdup_nomismatch.sh
Rscript $R_script $bash_script $filtered_sam_folder $sample_ID $rmdup_sam_folder $core $mismatch

echo
echo "Start removing duplicates..."
echo input folder: $all_output_folder/rmdup_sam
echo output folder: $all_output_folder/rmdup_sam_2
mkdir -p $all_output_folder/rmdup_sam_2
module unload python

bash_script=$script_path/sci3_rmdup_nomismatch.sh
filtered_sam_folder=$all_output_folder/rmdup_sam
rmdup_sam_folder=$all_output_folder/rmdup_sam_2
Rscript $R_script $bash_script $filtered_sam_folder $sample_ID $rmdup_sam_folder $core $mismatch`

My questions are:
1. What exactly is the difference between the sam files in rmdup_sam and rmdup_sam_2 folders?
2. Why are two rounds of the same script being run post barcode attachment?
3. The previous step is barcode+UMI attachment but the way I read it, we are only attaching the barcode, does that sound right?
I don't understand where the UMI is being defined. Am I missing something?

I am trying to understand which step is doing exactly what within sci3_main.sh and hence the questions.
Thank you!

Regards,
Deep

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant