User: Johan Zicola
Johan Zicola • 40
- 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
... 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 Zicola • 40
0
votes
1
answer
284
views
1
answers
... 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 Zicola • 40
0
votes
1
answer
411
views
1
answers
... 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 Zicola • 40
0
votes
1
answer
411
views
1
answers
... 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 Zicola • 40
4
votes
1
answer
411
views
5 follow
1
answer
... 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 ...
written 5 months ago by
Johan Zicola • 40
• updated
3 months ago by
sgalla32 • 30
0
votes
3
answers
6.9k
views
3
answers
... 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 Zicola • 40
0
votes
3
answers
6.9k
views
3
answers
... 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 Zicola • 40
1
vote
3
answers
1.0k
views
3
answers
... 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 Zicola • 40
3
votes
1
answer
1.5k
views
1
answers
... 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 Zicola • 40
0
votes
1
answer
4.7k
views
1
answers
... 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 Zicola • 40
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.
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