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

Update CUT&RUN tutorial #5104

Draft
wants to merge 2 commits into
base: main
Choose a base branch
from
Draft

Conversation

lldelisle
Copy link
Collaborator

Try to make the CUT&RUN tutorial out.
I ran the tutorial on usegalaxy.eu.

  • TODO: update image

@lldelisle lldelisle requested a review from a team as a code owner June 29, 2024 18:04
@lldelisle lldelisle marked this pull request as draft June 29, 2024 18:04
> > 1. ~55% and ~57%
> > 2. 100%
> > 1. ~55% for Read 1 and ~57% for Read 2
> > 2. The last line indicates that 3.5% of pairs have been removed.
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@heylf could you confirm it makes sense (the output of bowtie2 does not have 300k reads as input).

> This is what we call dovetailing and we want to consider this pair as a valid concordant alignment.
{: .comment}


Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I decided to remove this part because the tutorial is using TrimGalore and removes adapters even if they are a single base pair.

@@ -264,7 +245,7 @@ repetitive regions but keep reads falling into regions present in alternate loci
>
> > <solution-title></solution-title>
> >
> > 41.46+57.51=98.97%
> > 36.47+62.39=98.86%
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I don't know why these number changed (maybe the version of bowtie2)...

> > > 2. 454 peaks
> > > 3. Our precision is ~ 84%. A high precision is an indicator that we can predict true binding regions with high confidence.
> > > 1. 2,865 peaks (this is the number of lines of our last output)
> > > 2. About 4000 peaks (Rep1 and Rep2 has about 6-7 thousand peaks each)
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@heylf Here I am not sure of what you consider false positive but I guess they are peaks in a single replicate which are not overlapping the second replicate.

@@ -622,7 +603,7 @@ Let's find out the sequence motifs of the TF GATA1. Studies have revealed that G
> >
> > > <solution-title></solution-title>
> > >
> > > 1. !!!! Where is this info? !!!
> > > 1. Where is this info?
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Could not find this info

@hexylena
Copy link
Member

xref #5164

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

Successfully merging this pull request may close these issues.

2 participants