Moderator: Ian

gravatar for Ian
Ian4.9k
Reputation:
4,860
Status:
Trusted
Location:
University of Manchester, UK
Website:
http://personalpages.m...
Twitter:
IanDcalling
Scholar ID:
Google Scholar Page
Last seen:
5 days, 12 hours ago
Joined:
6 years, 12 months ago
Email:
d***********@gmail.com

Member of the Faculty of Life Science's bioinformatics core facility.

.

http://www.youtube.com/user/ManchesterBCF/videos?view=1&flow=grid for video tutorials.

.

I specialise in ChIP-seq analysis, but also run RNA-seq and other NGS applications. Data comes, almost exclusively from illumina sequencers, but I used to mostly analyse SOLiD sequence.

.

I am interested in the analysis, organisation and presentation of genomics data. I have a special place in my heart for the hunting of transcription factor binding sites.

Posts by Ian

<prev • 412 results • page 1 of 42 • next >
1
vote
1
answer
85
views
1
answers
Comment: C: How to choose Trimmomatic's parameter 'MINLEN '?
... I would not go down to 10bp. Remember the shorter the read length the greater the chance of a false positive match. You could go to 25bp if you are despirate. However, looking at the numbers of reads per length I think the majority of reads are >35bp anyway. Below 35bp numbers are only in the ...
written 12 days ago by Ian4.9k
0
votes
1
answer
95
views
1
answers
Answer: A: Differential chip-seq peak calling using macs2
... If you have paired reads then you need to use -f BAMPE, otherwise only the 5' read of each pair will be used and cross-correlation will be used to estimate fragment length, as opposed to using the actual fragment lengths between pairs. Have you compared the 'd' value for each analysis? I haven't ...
written 12 days ago by Ian4.9k
2
votes
1
answer
85
views
1
answers
Answer: A: How to choose Trimmomatic's parameter 'MINLEN '?
... The lowest I select is MINLEN:35 as this was the read length Illumina sequence for a long time. The important thing is ensure the base quality is good, especially at the 3' end. I usually use SLIDINGWINDOW:4:20. I personally do not use LEADING, as it has never been a problem. Also make sure you ...
written 12 days ago by Ian4.9k
2
votes
1
answer
153
views
1
answers
Answer: A: ChIP-seq Input replicate correlation
... Well, they are by definition your replicates, even if they don't look particularly similar. The difference between the two inputs might be down to coverage, unless they are already normalised. You could run S1_input (as a pseudo ChIP sample) against the S2_input (as an input), and visa versa, in ...
written 5 weeks ago by Ian4.9k
1
vote
3
answers
259
views
3
answers
Comment: C: Peak calling for broad histone-modification regions
... Looks interesting; a way of complementing traditional peak calling methods - thanks. ...
written 5 weeks ago by Ian4.9k
2
votes
1
answer
194
views
1
answers
Answer: A: Filter MACS peaks by enrichment/FDR?
... Have you considered that there are many binding regions for the protein you are investigating, some factors are ubiquitous. The FDR cut off sounds reasonable, although MACS2 suggests a q-value of 0.01 initially. I usually use FE >=10 as a cut off, but it is arbitrary. The best way is to observ ...
written 6 weeks ago by Ian4.9k
0
votes
1
answer
225
views
1
answers
Comment: C: prepare pair-ended reads for peak calling
... What 'd' values do you get? The mean fragment length for single end reads is determined by cross-correlation, and for paired-end the mean of the known fragment lengths is used. You could use the same value of 'd' by using --nomodel --extsize 200. ...
written 9 weeks ago by Ian4.9k
0
votes
1
answer
225
views
1
answers
Comment: C: prepare pair-ended reads for peak calling
... A BEDPE file created from the bowtie2 paired BAM file is a necessary intermediate from which you can create your BED file, as suggested above. ...
written 12 weeks ago by Ian4.9k
0
votes
6
answers
13k
views
6
answers
Comment: C: Extract Sequence From Fasta File Using Ids From A Separate File
... http://hgdownload.cse.ucsc.edu/admin/exe/ ...
written 3 months ago by Ian4.9k
0
votes
1
answer
173
views
1
answers
Comment: C: How to extend variable length intervals to the same final length?
... Thanks for the addition. After discussing this with a colleague this morning it was pointed out that finding the mid-point of each region and then extending out works equally well. I knew I had missed something! ...
written 3 months ago by Ian4.9k

Latest awards to Ian

Teacher 4 weeks ago, created an answer with at least 3 up-votes. For A: Download Genomic Ranges From A List Of Coordinates Via Web
Good Answer 7 weeks ago, created an answer that was upvoted at least 5 times. For A: Resource For Gff-Files?
Great Question 9 weeks ago, created a question with more than 5,000 views. For How To Plot Read Coverage Over Many Different Genomic Regions?
Appreciated 3 months ago, created a post with more than 5 votes. For A: Resource For Gff-Files?
Appreciated 4 months ago, created a post with more than 5 votes. For A: Resource For Gff-Files?
Popular Question 4 months ago, created a question with more than 1,000 views. For Masc - A Tool For Calculating Mean Fragment Size For Chip-Seq Peak Calling Analysis
Epic Question 4 months ago, created a question with more than 10,000 views. For What Is The Difference Between Samtools Mpileup And Pileup
Appreciated 5 months ago, created a post with more than 5 votes. For A: Resource For Gff-Files?
Great Question 5 months ago, created a question with more than 5,000 views. For Extracting Paired-End Reads Sitting In Different Chromosomes
Appreciated 6 months ago, created a post with more than 5 votes. For A: Resource For Gff-Files?
Popular Question 6 months ago, created a question with more than 1,000 views. For Mitochondrial Genome In Chip-Seq?
Popular Question 6 months ago, created a question with more than 1,000 views. For Mitochondrial Genome In Chip-Seq?
Popular Question 7 months ago, created a question with more than 1,000 views. For Mitochondrial Genome In Chip-Seq?
Popular Question 9 months ago, created a question with more than 1,000 views. For Mitochondrial Genome In Chip-Seq?
Teacher 9 months ago, created an answer with at least 3 up-votes. For A: Download Genomic Ranges From A List Of Coordinates Via Web
Appreciated 11 months ago, created a post with more than 5 votes. For A: Resource For Gff-Files?
Great Question 11 months ago, created a question with more than 5,000 views. For Extracting Paired-End Reads Sitting In Different Chromosomes
Student 13 months ago, asked a question with at least 3 up-votes. For Illumina'S New Reduced Resolution Q (Phred) Scores - Good Or Bad?
Great Question 13 months ago, created a question with more than 5,000 views. For Extracting Paired-End Reads Sitting In Different Chromosomes
Teacher 13 months ago, created an answer with at least 3 up-votes. For A: Download Genomic Ranges From A List Of Coordinates Via Web
Popular Question 13 months ago, created a question with more than 1,000 views. For Masc - A Tool For Calculating Mean Fragment Size For Chip-Seq Peak Calling Analysis
Great Question 14 months ago, created a question with more than 5,000 views. For Extracting Paired-End Reads Sitting In Different Chromosomes
Popular Question 14 months ago, created a question with more than 1,000 views. For Regular expression matching with Python and biopython SeqIO
Popular Question 14 months ago, created a question with more than 1,000 views. For Masc - A Tool For Calculating Mean Fragment Size For Chip-Seq Peak Calling Analysis

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 603 users visited in the last hour