User: Paul

gravatar for Paul
Paul1.2k
Reputation:
1,250
Status:
Trusted
Location:
European Union
Last seen:
20 hours ago
Joined:
4 years, 11 months ago
Email:
b********@gmail.com

about me

Posts by Paul

<prev • 159 results • page 1 of 16 • next >
1
vote
2
answers
86
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 18 days ago by Paul1.2k
0
votes
5
answers
7.0k
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 6 weeks ago by Paul1.2k
0
votes
2
answers
198
views
2
answers
Comment: C: Parser for Illumina Novaseq metadata
... Glad to help you (I was working on the same task last week :-))!!! ...
written 6 weeks ago by Paul1.2k
1
vote
2
answers
198
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 6 weeks ago by Paul1.2k
3
votes
2
answers
198
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 6 weeks ago by Paul1.2k
1
vote
1
answer
213
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 7 weeks ago by Paul1.2k
2
votes
0
answers
122
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 8 weeks ago by Paul1.2k
2
votes
1
answer
133
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 8 weeks ago by Paul1.2k
2
votes
3
answers
464
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 12 weeks ago by Paul1.2k
0
votes
0
answers
395
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 6 months ago by Paul1.2k

Latest awards to Paul

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

Help
Access

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