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
4 views

Properties that can be derived from amino acid sequences

What are the properties that can be derived from an amino acid sequence apart from those mentioned in the website? In the mentioned website, the properties that can be calculated from amino acid ...
3
votes
0answers
52 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
28 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
137 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
25 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! ...