Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In

Moderator: Matt Shirley

gravatar for Matt Shirley
Matt Shirley ♦ 8.9k
Reputation:
8,900
Status:
Trusted
Location:
Cambridge, MA
Website:
http://mattshirley.com/
Twitter:
mdshw5
Scholar ID:
Google Scholar Page
Last seen:
16 hours ago
Joined:
7 years, 11 months ago
Email:
m*****@gmail.com

Posts by Matt Shirley

<prev • 707 results • page 1 of 71 • next >
1
vote
7
answers
79k
views
7
answers
Comment: C: Correct Way To Parse A Fasta File In Python
... Thanks. I constantly forget that I'm a moderator. ...
written 16 days ago by Matt Shirley ♦ 8.9k
1
vote
7
answers
79k
views
7
answers
Comment: C: Correct Way To Parse A Fasta File In Python
... The example is *now* wrong, but almost 9 years ago when this was written it was perfectly valid. The Bio.Seq.Seq.tostring() method was removed in the latest Biopython release, and was deprecated a bit before: https://github.com/biopython/biopython/blob/654309121f2cc0c01dfff73cd3dec3a435d76fc2/DEPREC ...
written 19 days ago by Matt Shirley ♦ 8.9k
1
vote
0
answers
366
views
0
answers
Comment: C: RNA seq - TMM values
... I think the question was why didn't you go ahead and use edgeR/limma for the statistics part? Most people will use the linear modeling in limma/DESeq2 for statistical testing. ...
written 4 months ago by Matt Shirley ♦ 8.9k
2
votes
1
answer
303
views
1
answers
Answer: A: Get sequences from a BAM file as strings spaced out as described by cigar
... You're looking for a gapped representation of the alignment. I'm not sure how to accomplish this in R, but I addressed this issue in the python simplesam package, where you can call the [`gapped` method of a `Sam` object](https://simplesam.readthedocs.io/en/latest/#simplesam.Sam.gapped) to get the s ...
written 4 months ago by Matt Shirley ♦ 8.9k
0
votes
3
answers
10k
views
3
answers
Comment: C: Get chromosome sizes from fasta file
... `faidx` was failing because the cli script expects UCSC-style regions encoded as `contig:start-end`, where `start-end` is optional. Since the HLA alt contains contain `:` characters you ran into a parsing issue. If you want to use `pyfaidx` in the same way as biopython you can do: from pyfaidx ...
written 6 months ago by Matt Shirley ♦ 8.9k
0
votes
2
answers
480
views
2
answers
Comment: C: is this a potential mosaicism?
... It's a little more interesting to see this in blood tissue, since there's definitely more mosaicism in blood. If you have access to assays that sequence cell-free DNA, it might be interesting to see if there's an enrichment of the allele in the cfDNA, since this may be shed from some other maternal ...
written 7 months ago by Matt Shirley ♦ 8.9k
6
votes
2
answers
480
views
2
answers
Answer: A: is this a potential mosaicism?
... For the mother's mosaic allele to be passed to the child you she would need to have the mosaicism present in **her** germline cells. What tissue are you sequencing? If it's a buccal sample or skin punch then it's unlikely that the same developmental progenitor cell that gave rise to the sample you s ...
written 7 months ago by Matt Shirley ♦ 8.9k
0
votes
11
answers
49k
views
11
answers
Comment: C: How To Extract A Sequence From A Big (6Gb) Multifasta File ?
... By default the bash [read](http://tldp.org/LDP/Bash-Beginners-Guide/html/sect_08_02.html) command splits input on newline characters. So in this case `p` will be each incoming line, determined by the `\n` (newline) character. You can change the internal file separator to something like spaces or wha ...
written 8 months ago by Matt Shirley ♦ 8.9k
2
votes
5
answers
17k
views
5
answers
Comment: C: How Do I Convert Fa Files To Bed Format?
... You can accomplish this using `faidx -i bed genome.fa > out.bed`. For more details you can check out the documentation: https://github.com/mdshw5/pyfaidx#cli-script-faidx ...
written 8 months ago by Matt Shirley ♦ 8.9k
0
votes
4
answers
11k
views
4
answers
Comment: C: Converting Gvcf Files Into Vcf
... I would try Chris's suggestion, which should really be the accepted answer. ...
written 11 months ago by Matt Shirley ♦ 8.9k

Latest awards to Matt Shirley

Teacher 28 days ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Scholar 6 weeks ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Appreciated 12 weeks ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Good Answer 4 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Good Answer 4 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Scholar 6 months ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Good Answer 7 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Appreciated 7 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Teacher 7 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Good Answer 7 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Teacher 7 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Great Question 9 months ago, created a question with more than 5,000 views. For Pyfaidx: Efficient, "Pythonic" Random Access To Fasta Files Using Samtools-Compatible Indexing
Teacher 10 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Appreciated 10 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Appreciated 11 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Appreciated 11 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Popular Question 12 months ago, created a question with more than 1,000 views. For Masking low-complexity regions using pyfaidx MutableFastaRecord
Popular Question 13 months ago, created a question with more than 1,000 views. For Masking low-complexity regions using pyfaidx MutableFastaRecord
Scholar 14 months ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Appreciated 14 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Appreciated 15 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Good Answer 15 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Commentator 16 months ago, created a comment with at least 3 up-votes. For C: What Does 2X250Bp Buy Us?
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1758 users visited in the last hour