Moderator: Matt Shirley

gravatar for Matt Shirley
Matt Shirley8.9k
Reputation:
8,930
Status:
Trusted
Location:
Cambridge, MA
Website:
http://mattshirley.com/
Twitter:
mdshw5
Scholar ID:
Google Scholar Page
Last seen:
1 day, 1 hour ago
Joined:
8 years ago
Email:
m*****@gmail.com

Posts by Matt Shirley

<prev • 707 results • page 1 of 71 • next >
1
vote
7
answers
81k
views
7
answers
Comment: C: Correct Way To Parse A Fasta File In Python
... Thanks. I constantly forget that I'm a moderator. ...
written 6 weeks ago by Matt Shirley8.9k
2
votes
7
answers
81k
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 6 weeks ago by Matt Shirley8.9k
1
vote
0
answers
401
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 5 months ago by Matt Shirley8.9k
2
votes
1
answer
335
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 5 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 7 months ago by Matt Shirley8.9k
0
votes
2
answers
508
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 8 months ago by Matt Shirley8.9k
6
votes
2
answers
508
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 8 months ago by Matt Shirley8.9k
0
votes
11
answers
50k
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 9 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 9 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 12 months ago by Matt Shirley8.9k

Latest awards to Matt Shirley

Teacher 8 weeks ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Scholar 10 weeks ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Appreciated 3 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Good Answer 5 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Good Answer 5 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Scholar 7 months ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Good Answer 8 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Appreciated 8 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Teacher 8 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Good Answer 8 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Teacher 8 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Teacher 11 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
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 12 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Appreciated 12 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Popular Question 13 months ago, created a question with more than 1,000 views. For Masking low-complexity regions using pyfaidx MutableFastaRecord
Popular Question 14 months ago, created a question with more than 1,000 views. For Masking low-complexity regions using pyfaidx MutableFastaRecord
Scholar 15 months ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Teacher 15 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Teacher 15 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Appreciated 15 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Appreciated 16 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Good Answer 16 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Commentator 17 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: 1431 users visited in the last hour