User: Paul

gravatar for Paul
Paul1.3k
Reputation:
1,280
Status:
Trusted
Location:
European Union
Last seen:
1 day, 14 hours ago
Joined:
5 years, 1 month ago
Email:
b********@gmail.com

about me

Posts by Paul

<prev • 159 results • page 1 of 16 • next >
1
vote
2
answers
159
views
2
answers
Answer: A: split vcf by individual
... Hi, another way is to use `bcftools` query and `GATK SelectVariants`. Code is here: cat split_vcf.sh Script looks: #!/bin/bash genome=~/path/to/reference/genome.fa for sample in `bcftools query -l INPUT.vcf`; do for i in *.vcf do gatk SelectVar ...
written 11 weeks ago by Paul1.3k
0
votes
5
answers
7.1k
views
5
answers
Answer: A: What Is The Fastest Method To Determine The Number Of Positions In A Bam File Wi
... For very large BAM files (like WGS) I was tested Mosdepth [Link][1]. It is very fast and there are lot of option. Mosdepth can report the mean depth in 500-base windows genome-wide info under 9 minutes of user time with 3 threads. Example of syntax: mosdepth -n --fast-mode --by 500 sample.wgs ...
written 3 months ago by Paul1.3k
0
votes
2
answers
276
views
2
answers
Comment: C: Parser for Illumina Novaseq metadata
... Glad to help you (I was working on the same task last week :-))!!! ...
written 3 months ago by Paul1.3k
1
vote
2
answers
276
views
2
answers
Comment: C: Parser for Illumina Novaseq metadata
... Just fixed my typo in example - thanks. Yes those are shell command as alternative to parsing SAV data. ...
written 3 months ago by Paul1.3k
3
votes
2
answers
276
views
2
answers
Answer: A: Parser for Illumina Novaseq metadata
... Hi, look at the Illumina InterOP [Here][1]. You can use anaconda to install [here][2]. Syntax for data extraction is very easy and then you can plot data by gnuplot or R. example to get table (csv): interop_summary runfolder > /path/to/my/output example fot plotting: interop_plot_qs ...
written 3 months ago by Paul1.3k
1
vote
1
answer
433
views
1
answers
Answer: A: Bash scripting FastQC for multiple fastq files in multiple directories
... This could work: `find -maxdepth 5 -name "*fastq.gz" | parallel fastqc {} -o qc/` , then create easy script to move output to corresponding directory. ...
written 3 months ago by Paul1.3k
2
votes
0
answers
202
views
0
answers
Comment: C: speed up bcftools annotate command
... Hi, did you try GNU parallel? [Here][1] are some examples. I have good experiences to speed up all my processing tasks with parallel. [1]: https://gist.github.com/Brainiarc7/7af2ab5e88ef238da2d9f36b4be203c0 ...
written 3 months ago by Paul1.3k
2
votes
1
answer
229
views
1
answers
Answer: A: how to set metrics file when using MarkDuplicates in Picard
... Hi, my syntax looks like: picard -Xmx"$memory"g -Djava.io.tmpdir=$PWD/tmp -XX:ParallelGCThreads=$cpu MarkDuplicates I= ${p%.bam}.bam O= ${p%.bam}_rmdp.bam METRICS_FILE= ${p%.bam}_rmdp.txt CREATE_INDEX=true USE_JDK_DEFLATER=true USE_JDK_INFLATER=true Metrics should give you an statistics about ...
written 4 months ago by Paul1.3k
2
votes
3
answers
639
views
3
answers
Answer: A: How to convert NGS .tabular file to excel file?
... Hi, here is an python script which takes as first argument TSV and second argument creates xlsx file. import csv import sys reload(sys) sys.setdefaultencoding('utf8') from xlsxwriter.workbook import Workbook # Add some command-line logic to read the file names. tsv ...
written 5 months ago by Paul1.3k
0
votes
0
answers
457
views
0
answers
Comment: C: Parsing Blast output on linux/unix
... Did you try install [bioperl][1] to your system?? Bio::SearchIO is module of Bioperl. [1]: https://bioperl.org/INSTALL.html ...
written 8 months ago by Paul1.3k

Latest awards to Paul

Good Answer 18 days ago, created an answer that was upvoted at least 5 times. For C: Is bcf just a bgzipped vcf?
Epic Question 7 weeks ago, created a question with more than 10,000 views. For bash loop for alignment RNA-seq data
Popular Question 11 weeks ago, created a question with more than 1,000 views. For How to show whole blast contig
Appreciated 12 weeks ago, created a post with more than 5 votes. For A: How To Filter Mapped Reads With Samtools
Great Question 3 months ago, created a question with more than 5,000 views. For Extract specific reads from FASTQ files based on subsequence
Teacher 3 months ago, created an answer with at least 3 up-votes. For A: Fastq Convert To Fasta
Scholar 3 months ago, created an answer that has been accepted. For A: Help With Genotyping In My Vcf File
Scholar 4 months ago, created an answer that has been accepted. For A: Help With Genotyping In My Vcf File
Popular Question 4 months ago, created a question with more than 1,000 views. For How to show whole blast contig
Student 5 months ago, asked a question with at least 3 up-votes. For How to compare assembly from Trinity and Velvet/Oases
Popular Question 8 months ago, created a question with more than 1,000 views. For Filter out specific reads from FASTQ files
Teacher 8 months ago, created an answer with at least 3 up-votes. For A: Fastq Convert To Fasta
Popular Question 8 months ago, created a question with more than 1,000 views. For Compare my FASTA with Taxonomics group
Scholar 11 months ago, created an answer that has been accepted. For A: Help With Genotyping In My Vcf File
Student 12 months ago, asked a question with at least 3 up-votes. For How to compare assembly from Trinity and Velvet/Oases
Scholar 12 months ago, created an answer that has been accepted. For A: Help With Genotyping In My Vcf File
Popular Question 12 months ago, created a question with more than 1,000 views. For Filter out specific reads from FASTQ files
Commentator 15 months ago, created a comment with at least 3 up-votes. For C: For loop script
Student 15 months ago, asked a question with at least 3 up-votes. For How to compare assembly from Trinity and Velvet/Oases
Guru 16 months ago, received more than 100 upvotes.
Appreciated 16 months ago, created a post with more than 5 votes. For A: How To Filter Mapped Reads With Samtools
Student 16 months ago, asked a question with at least 3 up-votes. For How to compare assembly from Trinity and Velvet/Oases
Appreciated 16 months ago, created a post with more than 5 votes. For A: How To Filter Mapped Reads With Samtools
Good Answer 16 months ago, created an answer that was upvoted at least 5 times. For A: Compare two matrix
Good Answer 16 months ago, created an answer that was upvoted at least 5 times. For C: Is bcf just a bgzipped vcf?

Help
Access

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