Moderator: Matt Shirley

gravatar for Matt Shirley
Matt Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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 Shirley8.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?
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?

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