User: Michael Kosicki

Reputation:
80
Status:
New User
Location:
Wellcome Trust Sanger Institute
Last seen:
5 months ago
Joined:
6 months, 3 weeks ago
Email:
m***@sanger.ac.uk

Posts by Michael Kosicki

<prev • 8 results • page 1 of 1 • next >
3
votes
1
answer
251
views
1
answers
Answer: A: vcf files - extract SNPs
... If the first column is the "genes", then: bedtools intersect -wa -wb -a 1.vcf -b 2.vcf > shared.vcf grep -v '^#' shared.vcf | cut -f1 -d' ' shared.vcf | sort | uniq # list of unique genes Then, assuming it's the same reference (ie ref is the same on the same position) and that both vcf ...
written 6 months ago by Michael Kosicki80
0
votes
1
answer
428
views
1
answers
Comment: C: Changing the reference haplotype in VCF
... Sure, go ahead, thanks. ...
written 6 months ago by Michael Kosicki80
1
vote
1
answer
428
views
1
answers
Answer: C: Changing the reference haplotype in VCF
... Ok, I think I've found the solution myself: ## get mouse chr19 curl ftp://ftp.ensembl.org/pub/release-91/fasta/mus_musculus/dna/Mus_musculus.GRCm38.dna.chromosome.19.fa.gz > chr19.fa.gz zcat chr19.fa.gz > chr19.fa ## get BALBc SNPs curl ftp://ftp-mouse.sanger.ac.uk/cu ...
written 6 months ago by Michael Kosicki80
2
votes
1
answer
347
views
1
answers
Comment: C: Confidence Interval in R
... Look here: http://seqanswers.com/forums/showthread.php?p=161824 In short, every test functions under certain assumptions. Breaking those assumptions breaks the test. RNA-seq data breaks the assumption of the t-test that the data is drawn from a normal distribution. In practice, you'd fail to detect ...
written 6 months ago by Michael Kosicki80
2
votes
1
answer
347
views
1
answers
Comment: C: Confidence Interval in R
... Look up "edgeR" or "DESeq2" for differential gene expression testing on RNA-seq data. T-test is not appropriate in this situation, due to the way data is distributed (also, you probably need to normalize for sequencing depth between samples). P-value of a test is not a function of the confidence in ...
written 6 months ago by Michael Kosicki80
0
votes
1
answer
428
views
1
answers
Comment: C: Changing the reference haplotype in VCF
... Could you elaborate? This tool just "merges" reference fasta and a vcf file, and does not produce a new vcf file, as far as I can tell. ...
written 6 months ago by Michael Kosicki80
1
vote
2
answers
289
views
2
answers
Answer: A: Sanger Sequencing showing too much background noise?
... If I understood correctly, they have performed a PCR on the material and then Sanger sequenced the product. 1) The sequencing service provider usually specifies, whether the product needs to be purified (e.g. using beads or columns) or whether they can do it for you (usually for a fee) and how much ...
written 6 months ago by Michael Kosicki80
4
votes
1
answer
428
views
1
answer
Changing the reference haplotype in VCF
... A VCF file describes how to translate the reference genome into (one or more) alternative haplotypes. I need a VCF file that translates one of these haplotypes back into the reference. In other words, I need to "invert" the VCF file I have. If only SNPs were present, one would simply swap reference ...
vcf indel alignment snp written 6 months ago by Michael Kosicki80 • updated 6 months ago by Kevin Blighe28k

Latest awards to Michael Kosicki

Scholar 5 months ago, created an answer that has been accepted. For C: Changing the reference haplotype in VCF
Teacher 5 months ago, created an answer with at least 3 up-votes. For A: vcf files - extract SNPs
Scholar 6 months ago, created an answer that has been accepted. For C: Changing the reference haplotype in VCF

Help
Access

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