User: Johan Zicola

gravatar for Johan Zicola
Johan Zicola40
Reputation:
40
Status:
New User
Location:
Last seen:
1 week, 1 day ago
Joined:
1 year, 6 months ago
Email:
i******@hotmail.fr

Posts by Johan Zicola

<prev • 10 results • page 1 of 1 • next >
1
vote
3
answers
3.5k
views
3
answers
Answer: A: How To Find The Locations Of A Short Specific Sequence In A Genome With 1 Or 2 M
... Using **Bowtie** (for example v1.2.2 here) to find off-targets for defined CRISPR-Cas9 target sequences: Make the Bowtie index for your genome (fasta file format) bowtie-build -f genome.fa genome_prefix Search for your target sequence by allowing 1 mismatch (for your N) with the flag `-n 1` ...
written 5 weeks ago by Johan Zicola40
0
votes
1
answer
284
views
1
answers
Answer: A: How to identify which SNPs in my dataset are genome wide significant
... It seems your QQplot and Manhattan plot were generated with the R package [qqman](https://github.com/stephenturner/qqman). You can filter out the SNPs above the threshold you define, the blue and red lines being set by default at -log10(1E-5) and -log10(5E-8) (check [documentation](https://cran.r-pr ...
written 3 months ago by Johan Zicola40
0
votes
1
answer
411
views
1
answers
Comment: C: Why main diagonal of gemma relatedness matrix is not equal to 1
... I thought the terms kinship matrix and relatedness matrix could be used interchangeably. Gemma documentation should be updated to use the accurate term. I understood that a value >0.5 (self-kinship) or >1 (self-relatedness) could be explained by consanguinity but I still don't understand why I ...
written 3 months ago by Johan Zicola40
0
votes
1
answer
411
views
1
answers
Comment: C: Why main diagonal of gemma relatedness matrix is not equal to 1
... Thanks for your comment. Looking around, I could find that usually the main diagonal is made of 1 or values above 1 if consanguinity. The figure 1 of the [Bae et al 2014][1] is illustrating well what I expect from a kinship matrix. Since individuals compared to themselves are based on the same SNP c ...
written 5 months ago by Johan Zicola40
4
votes
1
answer
411
views
5 follow
1
answer
Why main diagonal of gemma relatedness matrix is not equal to 1
... I am running GWAS analyses with gemma that seem to work but when I open the centered relatedness matrix (.cXX), I do get a symmetric matrix but the main diagonal is not made of 1s ... why? The [manual][1] p.11 part 3.3.1 shows an example of matrix with a main diagonal which is not made of 1s. I al ...
gwas relatedness matrix kinship matrix gemma written 5 months ago by Johan Zicola40 • updated 3 months ago by sgalla3230
0
votes
3
answers
6.9k
views
3
answers
Comment: C: How to extract specific chromosome from vcf file
... This won't work, separating chromosome names by commas neither (though it works for `bcftools view --regions 1,2,3`). Do rather `vcftools --gzvcf --chr 1 --chr 2 [ etc. until 23] --recode --out subset_chr1-23` ...
written 7 months ago by Johan Zicola40
0
votes
3
answers
6.9k
views
3
answers
Comment: C: How to extract specific chromosome from vcf file
... Note that this method is better than grep as it includes the VCF header. However, it won't change the header of the VCF file so the unselected chromosomes will still have their ID line, e.g *##contig=*. So don't rely on `bcftools view -h subset.vcf` to verify what chromosomes are left in your VCF fi ...
written 7 months ago by Johan Zicola40
1
vote
3
answers
1.0k
views
3
answers
Answer: A: Generating random DNA sequence and paired-end alignment
... I wrote a python script with the different functions you would need to test structural variation calling on either randomly generated fastq files or fastq files generated based on a given specified fasta file. Find the script and documentation on https://github.com/johanzi/fastq_generator ...
written 7 months ago by Johan Zicola40
3
votes
1
answer
1.5k
views
1
answers
Answer: A: Does vcf-consensus support indels?
... Not sure about vcftools-consensus but indels are taken into account with "bcftools consensus" and is probably faster than vcftools: bcftools consensus --fasta-ref --sample -o See docs: https://samtools.github.io/bcftools/bcftools.html ...
written 10 months ago by Johan Zicola40
0
votes
1
answer
4.7k
views
1
answers
Comment: C: Error in sam to bam conversion after bowtie alignment
... There are 2 possibilities in Bowtie to get an output in SAM format: - Add the tag `-S ` if one wants directly the SAM file (space consuming) - In case one wants to pipe the sam output directly in samtools to get a sorted bam file (smaller), add flag --sam such as `bowtie --sam | samtools view -b ...
written 14 months ago by Johan Zicola40

Latest awards to Johan Zicola

Teacher 3 months ago, created an answer with at least 3 up-votes. For A: Does vcf-consensus support indels?
Supporter 5 months ago, voted at least 25 times.

Help
Access

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