User: apa@stowers

gravatar for apa@stowers
apa@stowers340
Reputation:
340
Status:
Trusted
Location:
Kansas City
Last seen:
3 months, 2 weeks ago
Joined:
6 years, 2 months ago
Email:
a**@stowers.org

Posts by apa@stowers

<prev • 60 results • page 1 of 6 • next >
0
votes
3
answers
3.6k
views
3
answers
Answer: A: Checking Version of Bioconductor
... As of about R 3.2.0, the way to get the associated Bioconductor version is to run the non-exported built-in function: \> tools:::.BioC_version_associated_with_R_version() [1] '3.1' Mind the '.' in front of the 'BioC'. That said, it may be possible for some installed packages to have been sourc ...
written 3 months ago by apa@stowers340
0
votes
4
answers
3.2k
views
4
answers
Answer: A: bedtools coverage gives error received illegal bin number, but that number isn't
... This can also happen if you have double-0 intervals, AND are using bedTools > 2.24.0. Basically if one of your bed files has a line like "chrX 0 0 ...", intersectBed v2.24.0 will run it successfully, but intersectBed v2.26.0 will throw the error "Received illegal bin number 4294967295", yes w ...
written 9 months ago by apa@stowers340
0
votes
2
answers
819
views
2
answers
Comment: C: tophat aligns to wrong strand
... When switching strands in the bam file, you would need to switch all strands, globally -- otherwise you are right, mate pairs would suddenly appear as discordant. IGV, in particular will color discordant pairs differently, but that is all. Other software may or may not pay attention to the concord ...
written 16 months ago by apa@stowers340
0
votes
2
answers
819
views
2
answers
Comment: C: tophat aligns to wrong strand
... > However, it's not correct that non-Tuxedo pipelines will read > alignments backwards. What I meant was that it reads the *correct* strand of the read backwards, sorry for my semantic glitch. The true strand of the RNA fragment is being reported wrong, because the library was made from rev- ...
written 16 months ago by apa@stowers340
0
votes
1
answer
377
views
1
answers
Comment: C: What is the difference between structural RNA-derived reads and non-structural R
... If anyone still does duplex-specific nuclease library normalization, it can have an effect too. It supposedly doesn't degrade sRNA species with hairpins or other secondary structure, but when we tried it, appeared to do it anyway. ...
written 16 months ago by apa@stowers340
0
votes
3
answers
1.5k
views
3
answers
Comment: C: RPKM and FPKM directly comparable?
... You can convert them, in principle. The only difference between F and R is that F=R/N where N is the number of reads per fragment. So for single-end data, RPKM and FPKM are identical. For paired-end, FPKM = RPKM/2. In principle. If you used an FPKM estimator like Cufflinks or RSEM, then it will ...
written 16 months ago by apa@stowers340
1
vote
1
answer
871
views
1
answers
Comment: C: HT-seq Counts Very High Number of No Feature
... Ah hold on, I was only partially right. Even with "--library-type fr-firststrand" I think the strand will still be reversed in the bit flags. Tophat will just add the "XS" tag for strand -- which will contradict the bit flags -- but only Cufflinks reads this tag. Other software will still use the ...
written 16 months ago by apa@stowers340
0
votes
2
answers
819
views
2
answers
Answer: A: tophat aligns to wrong strand
... If I recall, use of "--library-type fr-firststrand" does not actually change the SAM bit flags. Instead, it adds the XS strand tag to the SAM record -- basically for Cufflinks to use. Non-Cufflinks programs will still read the alignments backwards. To really fix it, you'll need to edit the bit fl ...
written 16 months ago by apa@stowers340
0
votes
2
answers
2.1k
views
2
answers
Answer: C: Error in read.table(file = file, header = header, sep = sep, quote = quote, :
... It is probably a problem with your header row. If the first row in the file has fewer fields than the second row, read.delim() (and all other functions built on it) will automatically assign column 1 to be rownames. ...
written 16 months ago by apa@stowers340
1
vote
1
answer
1.2k
views
1
answers
Answer: A: Data resources for Chromatin accessibility (DNAse, FAIRE, ATAC - seq) and Transc
... You are basically looking for a public repository (google "sequencing public repository", for instance) Here are a few more big ones which handle next-gen sequencing data: - GEO ([https://www.ncbi.nlm.nih.gov/geo/][1]) - SRA ([https://www.ncbi.nlm.nih.gov/sra][2]) - ArrayExpress ([https://www. ...
written 16 months ago by apa@stowers340

Latest awards to apa@stowers

Scholar 16 months ago, created an answer that has been accepted. For A: lift coordinates mapping to_alt chromosomes in hg38
Supporter 16 months ago, voted at least 25 times.
Scholar 16 months ago, created an answer that has been accepted. For A: lift coordinates mapping to_alt chromosomes in hg38
Good Answer 6.2 years ago, created an answer that was upvoted at least 5 times. For A: Is There Any Reason To Do De Novo Transcript Assembly If A Reference Is Availabl
Teacher 6.2 years ago, created an answer with at least 3 up-votes. For A: Is There Any Reason To Do De Novo Transcript Assembly If A Reference Is Availabl

Help
Access

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