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 ♦ 9.2k
Reputation:
9,210
Status:
Trusted
Location:
Cambridge, MA
Website:
http://mattshirley.com/
Twitter:
mdshw5
Scholar ID:
Google Scholar Page
Last seen:
2 days, 15 hours ago
Joined:
8 years, 7 months ago
Email:
m*****@gmail.com

Posts by Matt Shirley

<prev • 708 results • page 1 of 71 • next >
1
vote
1
answer
16k
views
1
answers
Comment: C: Ibd (Identity By Descent) And The Chosen Value Of Pi_Hat
... The two measures are quantifying the same thing: the percentage of a genome that is shared by decent with an ancestor. However PIHAT and Kcoeff can be calculated using different window sizes for consecutive informative SNPS, or can be calculated with or without CNV information, which can help discri ...
written 5 months ago by Matt Shirley ♦ 9.2k
1
vote
7
answers
91k
views
7
answers
Comment: C: Correct Way To Parse A Fasta File In Python
... Thanks. I constantly forget that I'm a moderator. ...
written 8 months ago by Matt Shirley ♦ 9.2k
2
votes
7
answers
91k
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 8 months ago by Matt Shirley ♦ 9.2k
1
vote
0
answers
570
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 12 months ago by Matt Shirley ♦ 9.2k
2
votes
1
answer
540
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 12 months ago by Matt Shirley ♦ 9.2k
0
votes
4
answers
14k
views
4
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 14 months ago by Matt Shirley ♦ 9.2k
0
votes
2
answers
654
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 15 months ago by Matt Shirley ♦ 9.2k
6
votes
2
answers
654
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 15 months ago by Matt Shirley ♦ 9.2k
0
votes
12
answers
56k
views
12
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 16 months ago by Matt Shirley ♦ 9.2k
2
votes
5
answers
19k
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 16 months ago by Matt Shirley ♦ 9.2k

Latest awards to Matt Shirley

Good Answer 15 days ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Appreciated 16 days ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Scholar 8 weeks ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Teacher 8 weeks ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Appreciated 11 weeks ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Teacher 12 weeks ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Scholar 5 months ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Teacher 9 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Scholar 9 months ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Appreciated 10 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Good Answer 12 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Good Answer 12 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Scholar 14 months ago, created an answer that has been accepted. For A: How to use pygr? worldbase doesn't return anything
Good Answer 15 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Appreciated 15 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Teacher 15 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Good Answer 15 months ago, created an answer that was upvoted at least 5 times. For A: How Can I Do Principal Components Analysis ?
Teacher 15 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Great Question 17 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 18 months ago, created an answer with at least 3 up-votes. For A: What Does 2X250Bp Buy Us?
Appreciated 18 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Appreciated 19 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Appreciated 20 months ago, created a post with more than 5 votes. For A: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
Popular Question 20 months ago, created a question with more than 1,000 views. For Masking low-complexity regions using pyfaidx MutableFastaRecord
Popular Question 22 months ago, created a question with more than 1,000 views. For Masking low-complexity regions using pyfaidx MutableFastaRecord
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2346 users visited in the last hour