User: haibol

gravatar for haibol
haibol90
Reputation:
90
Status:
Trusted
Location:
Last seen:
1 month ago
Joined:
8 months, 3 weeks ago
Email:
h*****@iastate.edu

Posts by haibol

<prev • 9 results • page 1 of 1 • next >
2
votes
3
answers
157
views
3
answers
Answer: A: Error in running ATACseqQC
... 1. In out-pdf functions such as fragSizeDist, PTscore, only when I exit R session can I get non-empty pdf document, and if I don't exit R, the out-pdf is always empty, so that run a function, I have to go into R, and then exit to get the target file. Is this the way it was done? Or there other way ...
written 4 weeks ago by haibol90
1
vote
0
answers
114
views
0
answers
Comment: C: ATACseqQC function bamSplit returns error when input .bam is from mouse
... Hello, If you check the document for the function splitBam, you will see: splitBam(bamfile, tags, index = bamfile, outPath = NULL, txs, genome, conservation, positive = 4L, negative = 5L, breaks = c(0, 100, 180, 247, 315, 473, 558, 615, Inf), labels = c("NucleosomeFree", "inter1", "mononu ...
written 3 months ago by haibol90
0
votes
3
answers
207
views
3
answers
Comment: C: ATAC-QC shiftGAlignmentsList function
... Shuang, did you use the mm10 reference genome from the same source, such as NCBI, UCSC Genome Browser or Ensembl, to do alignment and the ATACseqQC analysis? I don't think ATACseqQC is species-aware. Haibo ...
written 3 months ago by haibol90
0
votes
3
answers
207
views
3
answers
Comment: C: ATAC-QC shiftGAlignmentsList function
... Shuang, did you use the mm10 reference genome from the same source, such as NCBI, UCSC Genome Browser or Ensembl, to do alignment and the ATACseqQC analysis? I don't think ATACseqQC is species-aware. Haibo ...
written 3 months ago by haibol90
1
vote
1
answer
1.2k
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Hello Gary, You new plots based on filtered BAM looks better: many spiky peaks due to duplicates are gone. But you still missing the 38-100bp fragments, which was permanently removed by previous wet lab steps. A good tutorial for analyzing ATAC-seq data is here: https://informatics.fas.harvar ...
written 8 months ago by haibol90
2
votes
1
answer
1.2k
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Hello Gary, You can find the scripts for filtering the unwanted reads from the supplementary file 1. Here are some script excerpts: #!/bin/bash #BSUB -n 8 # minimal numbers of processors required for a parallel job #BSUB -R rusage[mem=16000] # ask for memory 16G #BSUB - ...
written 8 months ago by haibol90
1
vote
1
answer
1.2k
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Hello Gary, Yes, I agree with you that there is no or very few mitochondrial DNA sequence in the Chinese bulbul genome assembly. Haibo ...
written 8 months ago by haibol90
1
vote
1
answer
1.2k
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Gary, Your command did not remove mitochondrial read alignment; it only output mitochondrial read alignments. One thin you can do is to do blast by using the black-capped bulbul mitochondrial genome as query to find out whether there is mitochondrial genome assembly in the draft genome assembly. ...
written 8 months ago by haibol90
1
vote
1
answer
1.2k
views
1
answers
Comment: C: How about our ATAC-Seq data quality analyzed by ATACseqQC
... Hello Gary, I am one of the co-developer of the ATACseqQC package. Your library fragment size plot is not optimal. It lacks fragments less than 100bp, which should be from OPEN chromatin regions. On the other hand your plots did not show a apparent ladder pattern of mono-,bi-, tri-,...nucleos ...
written 8 months ago by haibol90

Latest awards to haibol

No awards yet. Soon to come :-)

Help
Access

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