User: vmicrobio

gravatar for vmicrobio
vmicrobio230
Reputation:
230
Status:
Trusted
Location:
France
Last seen:
7 hours ago
Joined:
1 year, 11 months ago
Email:
v**************@gmail.com

Posts by vmicrobio

<prev • 92 results • page 1 of 10 • next >
0
votes
1
answer
167
views
1
answers
Comment: C: How can I make a .bed file from my fasta?
... are your sequences located on the genbank? if so, you can retrieve genebank files and execute a script such as : #!/usr/bin/env python # encoding: utf-8 """ bed_from_genbank.py grab the gene records from a genbank file (edit for other record types). - requires ...
written 4 weeks ago by vmicrobio230
0
votes
1
answer
149
views
1
answers
Comment: C: idxstats non mapped reads
... that's pretty clear, thanks! ...
written 3 months ago by vmicrobio230
0
votes
1
answer
350
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 3 months ago by vmicrobio230
3
votes
1
answer
149
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 3 months ago by vmicrobio230 • updated 3 months ago by Devon Ryan86k
0
votes
11
answers
8.0k
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 5 months ago by vmicrobio230
0
votes
1
answer
396
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 5 months ago by vmicrobio230
0
votes
4
answers
275
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 5 months ago by vmicrobio230
0
votes
3
answers
681
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 7 months ago by vmicrobio230
0
votes
3
answers
681
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 7 months ago by vmicrobio230
1
vote
7
answers
1.4k
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 9 months ago by vmicrobio230

Latest awards to vmicrobio

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