User: haro

gravatar for haro
haro230
Reputation:
230
Status:
Trusted
Location:
France
Last seen:
4 days, 8 hours ago
Joined:
1 year, 8 months ago
Email:
v**************@gmail.com

Posts by haro

<prev • 91 results • page 1 of 10 • next >
0
votes
1
answer
68
views
1
answers
Comment: C: idxstats non mapped reads
... that's pretty clear, thanks! ...
written 7 days ago by haro230
0
votes
1
answer
264
views
1
answers
Comment: C: Create a custom exome .genome for IGV
... so in order to do high throughput there is no ways to do it using command lines? ...
written 7 days ago by haro230
3
votes
1
answer
68
views
1
answer
idxstats non mapped reads
... Hi, I'm doing alignment using bwa, trying to determine if sequences are mapped or not to a reference. I tried to compare fastq files to one reference and obtained my output in csv using samtools idxstats. I got this output: AE005174.2 5528445 2080495 6121 * 0 0 92628 and do not understand ...
idxstats alignment samtools written 7 days ago by haro230 • updated 7 days ago by Devon Ryan82k
0
votes
11
answers
7.1k
views
11
answers
Answer: A: Downsampling dataset with more than 60 million reads
... you may try something in awk : cat file.fastq | awk '{ printf("%s",$0); n++; if(n%4==0) {printf("\n");} else { printf("\t");} }' | awk -v k=8000 'BEGIN{srand(systime() + PROCINFO["pid"]);}{s=x++ "downsamp-file.fastq"}' ...
written 9 weeks ago by haro230
0
votes
1
answer
138
views
1
answers
Comment: C: conversion of GFF3 formate to BED format
... try this : [https://www.biostars.org/p/61386/][1] [1]: https://www.biostars.org/p/61386/ ...
written 9 weeks ago by haro230
0
votes
4
answers
178
views
4
answers
Answer: A: database for reference genome
... Hello, You may retrieve your reference genomes using esearch, with something like this: .../edirect/esearch -db nucleotide \ -query "name_of_your_organism[organism] AND genome[title]" \ | .../edirect/efetch -format fasta > references.fasta ...
written 9 weeks ago by haro230
0
votes
3
answers
323
views
3
answers
Comment: C: Read a fasta file | Java
... you can create a class FastaSequence containing the code above, add a 'getHeader' and 'getSequence' and then return only the sequence for your use ...
written 4 months ago by haro230
0
votes
3
answers
323
views
3
answers
Answer: A: Read a fasta file | Java
... Hi ozdavidd, you may try this : private void readFastaFile(File fastaFile) { InputStream flux; String line; try { flux = new FileInputStream(fastaFile); InputStreamReader lecture = new InputStreamReader(flux); BufferedReader buff = ne ...
written 4 months ago by haro230
1
vote
7
answers
1.1k
views
7
answers
Answer: A: shannon entropy score
... Give a try to biojava: import java.util.*; import org.biojava.bio.dist.*; import org.biojava.bio.seq.*; import org.biojava.bio.symbol.*; public class Entropy { public static void main(String[] args) { Distribution dist = null; try { //crea ...
written 6 months ago by haro230
1
vote
1
answer
184
views
1
answers
Answer: A: Mean coverage in Ion Torrent PGM
... Hi Abdul Rafay Khan, You can use samtools depth: samtools depth -b file.bed -a file.bam > coverage.csv Also, a plugin named 'coverageAnalysis' exists directly on the Proton machine ...
written 6 months ago by haro230

Latest awards to haro

Popular Question 4 weeks ago, created a question with more than 1,000 views. For 16S rRNA seq software
Scholar 10 months ago, created an answer that has been accepted. For A: What is the best way to start learning R?
Scholar 12 months ago, created an answer that has been accepted. For A: What is the best way to start learning R?
Supporter 16 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: 1475 users visited in the last hour