BOL

MIRA follows for the ref based assembly (like mapping the reads on the reference with xxx, etc.) and if its possible to have a figure with the workflow.

Hi BOL community,

I, Geetanjali, recently completed Msc in bioinformatics. I would like to know more about job offers in Bioinformatics.

Please help.

I am using this perl code right nowadays.

#!/usr/bin/perl
use warnings;
use strict;
open(FASTA, "<", seq.fa) or die("Can't open\n");
my %singleLineSequences;
my $sequence_id;
while(){
    my $line = $_; chomp($line);
    if ($line =~ m/^>(\S+)/){
        $sequence_id = $1;  
        $singleLineSequences{$sequence_id} = ""; 
        }
    else {
        $singleLineSequences{$sequence_id} = $singleLineSequences{$sequence_id} . $line;
        }
   }


foreach my $sequence_entry (keys %singleLineSequences){
    my $currentSequence = $singleLineSequences{$sequence_entry};
    my $lengthSequence = length($currentSequence);
    print $sequence_entry . "," . $lengthSequence . "\n";
}

I am doing megablast of a contig against NR database, but unable to understand this

-outfmt '6 qseqid staxids qstart qend sstart send qseq sseq evalue length'

Can you please help me to undestand all these.

Until there is a radical change in the way that academic credit is given, the principal record of scientific research is still the peer-reviewed publication. Given that software is a fundamental part of doing science in the digital age, the question we are often asked is: where can I publish papers which are primarily focused on my scientific software?

Is there any database exists which provides RNAseq data for different cancers stages?

I run SatsumaSynteny on my server using following command:

[jit@hm satsuma-code-0]$ ./SatsumaSynteny -q Genome/renamedG.fa -t Genome/genome_v4.fasta -o Genome/OutFile -m 128 -ni 10 -n 256 -chain_only

But it kill the program with followng error: 

./SatsumaSynteny: /lib64/libc.so.6: version `GLIBC_2.14' not found (required by ./SatsumaSynteny)

Can anyone please help to resolve this.

My server detail:

[jit@hm satsuma-code-0]$ lscpu
Architecture: x86_64
CPU op-mode(s): 32-bit, 64-bit
Byte Order: Little Endian
CPU(s): 12
On-line CPU(s) list: 0-11
Thread(s) per core: 1
Core(s) per socket: 6
Socket(s): 2
NUMA node(s): 2
Vendor ID: AuthenticAMD
CPU family: 16
Model: 8
Stepping: 1
CPU MHz: 2600.186
BogoMIPS: 5200.03
Virtualization: AMD-V
L1d cache: 64K
L1i cache: 64K
L2 cache: 512K
L3 cache: 5118K
NUMA node0 CPU(s): 0,2-6
NUMA node1 CPU(s): 1,7-11

Is there any easiest known way to extract organism's chromosome length?

I received a new genome to annotated. I need your help to quickly count the occurance of all the letters in fasta file.