Analysis of DNA, RNA or protein (amino acid) sequences. Here sequences are compared so analyze similarities and conserved regions. Sequence analysis is usually done by specialized software.

learn more… | top users | synonyms

0
votes
1answer
18 views

Marker validation using transcriptome and genomic sequences derived from a single cell

I conducted a SSR marker analysis using bioinformatic tool on some RNA-seq data of human tissues. Now, my supervisor believed that we have to validate these SSR on a transcriptomic and genomic ...
3
votes
0answers
51 views

Efficient algorithm to calculate various population divergence statistics

Intro and description of the data I am simulating the evolution of very long DNA sequences. The model works well, is performant and will output data in the following kind of fasta format ...
3
votes
0answers
27 views

How can I create R8 (homopolymer repeat) filter without using illumina pipeline?

illumina instruments have built-in -or online- analysis software for variant analysis (CASAVA). This software can filter out the false positive variants near the homopolymer repeats (AAAAAAAA) and ...
3
votes
0answers
129 views

Good poly-A filtering rules or tools

I am aligning a large number of ESTs. It seems poly-A tails show in many different ways. In addition to occurring at the very end, they can be flanked by the cloning sequence one one end, or have ...
1
vote
0answers
136 views

Is the default kmer 25 value in Trinity still ok for assembly of 150bp reads from Illumina HiSeq 2500?

I will soon be receiving Illumina HiSeq 2500 data (150bp PE stranded reads). In the past we have been using Trinity as our assembler of choice, but it uses a default kmer value of 25. I believe this ...
0
votes
0answers
24 views

How do I select the sequence identity threshold in cd-hit to cluster proteins for core genome comparative analysis?

Can anyone suggest a method for selecting the sequence identity threshold for protein clustering in cd-hit? If I use too high of a threshold, then the number of protein clusters shared between three ...
0
votes
0answers
7 views

Why low variant frequency in exomeseq data (ANNOVAR)?

I've got an annotated ANNOVAR file with 300.000 nucleotide variants of 4 genomes. There's one thing I don't understand when looking at the data: Why are the variant frequencies so low in comparison to ...
0
votes
0answers
43 views

Demultiplexing symmetric PacBio reads

I have PacBio CCS.h5 and the corresponding fasta and fastq files and I would like to demultiplex them. Does anyone know of how this can be done in the absence of bas.h5 files. Thanks for your help! ...