Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [alignment]

These questions are about sequence alignment.

2
votes
1answer
36 views

Alignment for predicting DNA hybridization?

I am currently working on a Computer Science project where we are trying to build a large set of orthogonal single-stranded DNA sequences. The goal would be to ensure that when put in solution, the ...
2
votes
2answers
133 views

Frequency of specific viral sequence in .BAM or .fastq

I was wondering if anyone had any experience 'counting' the frequency of a particular viral sequence in an individual's fastq sequence. So basically, counting the number of occurrences of a particular ...
2
votes
1answer
97 views

How do I find split reads?

How can I detect a split read in a BAM file? Is there any sign in the CIGAR string that describes split read?
3
votes
1answer
59 views

Pipeline for extracting gene from multiple genomes for use in HyPhy selection analyses?

I have been trying to obtain some preliminary data from HyPhy selection analyses to inform a larger project. I have obtained a number of assembled mammalian genomes from NCBI with the initial goal of ...
6
votes
2answers
85 views

Definition of “seed” in sequence alignment

I would like to know what is meant by "seed" for various sequence aligners. How is it important?
0
votes
0answers
32 views

What Ensembl genome version of wheat should I use for alignments? [duplicate]

There are many genome files available from Ensembl. Which one is the best to use/download for wheat RNA seq differential gene expression analysis? Also, tell me which tool or aligner or mapper to use? ...
1
vote
0answers
51 views

Aligning nucleotide sequences in APE software

I am very new to APE software. I am trying to align complementary oligo-DNA strands using APE on macOS. I have both the forward and reverse sequence of an oligo-DNA insert in two separate files. But I ...
3
votes
3answers
71 views

Alignment with arbitrary number of mismatches or gaps

I have 23bp long reads and want to find all possible alignments of them to the human genome (hg19, hg38) for an arbitrary number of mismatches (<7), possibly also small indels. I've read in ...
3
votes
1answer
27 views

Determining Read Groups

Which Read Groups are correct: ...
4
votes
1answer
216 views

How to get fasta alignment file from SAM/BAM file?

I am not talking about consensus sequence, I know how to get consensus sequence using mpileup in samtools/bcftools. As I understand , SAM/BAM files are basically sequence alignment format so it's ...
3
votes
0answers
39 views

convert supplementary reads to primary in sam or bam

I have a problem with a tool, that possibly ignores supplementary reads. I want to find out a bit how this tool works by converting all supplementary reads to primary reads, and then change the names ...
4
votes
1answer
135 views

Can blat use more than one core/CPU to speed up the alignment?

I am using BLAT to align two versions of the genome of C. elegans. I can see in the Activity Monitor of my Mac Book Pro High Sierra that blat is using 100% of a CPU....
3
votes
1answer
52 views

Multiple sequence pairwise global alignment against reference sequence

Background I have a fasta file with protein sequences. The first sequence in the fasta file is the reference sequence and every subsequent sequence is something I want to align to the reference. The ...
5
votes
5answers
160 views

Identifying Indels from Chromatograms

I have around 100 chromatograms (.ab1 files) from Sanger sequencing a genome at loci believed to have an indel. I'm new to interpreting this kind of data in ...
5
votes
1answer
81 views

Why is a PacBio read length larger than the aligned reference region?

I recently had some Iso-Seq sequencing done on my organism catfish on the new Sequel platform and got weird alignments for a size selected 4 Kilobase and up fraction after running the isoseq3 pipeline....
2
votes
2answers
158 views

How to count the number of mapped read in 100-bp window from a BAM/SAM file

Although I know how to get total number of mapped read using samtools flagstat (samtools flagstat file_sorted.bam) but I want to count total number of mapped read ...
5
votes
1answer
97 views

How to output all sequences with bwa mem, not `*`?

I've been running bwa mem -a for alignment, using the -a flag---this will output all alignments for SE or unpaired PE I've ...
1
vote
0answers
13 views

Finding the conserved region of protein in a given set of Kingdom

We have to find the conserved region of proteins in a given set of kingdom and we have to give more weightage to distantly related organism. What approach I should take to solve this problem?
5
votes
1answer
73 views

Does the “.full.aln” file produced by snippy-core contain all bases of my input sequences aligned to the reference genome?

I have a number of sequences and a reference genome. I used snippy to align each individual sequence with the reference genome. I then used ...
5
votes
3answers
233 views

Total reads aligning to each reference within a bam file

I have two PCR amplicons that have been multiplexed and sequenced using the nanopore minion. I have aligned the fastq reads using minimap2 with a reference file containing both amplicon sequences and ...
3
votes
1answer
62 views

cDNA and alignment mapping

I am confused with RNA seq alignment. My understanding is that after Mature mRNA is isolated from the cell, it is then fragmented and using reverse transcriptase enzyme a cDNA copy is created which is ...
4
votes
2answers
78 views

How to align output of grep --color=always? (To QC fasta/fastq files)

Grepping out short sequences from a fasta or fastq file is a really useful way to look at sequencing data. Using the option --color=always makes this even more ...
0
votes
1answer
66 views

How to calculate Gene Ontology terms in python

I am testing different Protein-Protein Interaction networks alignment tools. I observed that all alignment tools show some ...
3
votes
0answers
46 views

How to assign the best gap penalty and gap extension penalty using BLOSUM65

For an assignment I must do a pairwise optimal local alignment using BLOSUM65 and five protein sequences. The algorithm I want to use is the Smith-Waterman. Context protein sequencing using Blastp: ...
6
votes
2answers
235 views

Sequence alignment using Markov Model

I am learning about applying Markov model to sequence alignment. The prof says that the transition probabilities from a gap-residue alignment to a residue-gap alignment and vice versa are both 0. Is ...
4
votes
1answer
104 views

How to map short sequences to long reads, recovering all multiply-mapped high-quality matches

The dilemma: I have a problem where I have around one million short sequences (21 bp to several 100s of basepairs) for which I need to identify all occurrences of in 20-30x coverage noisy long reads (...
5
votes
2answers
90 views

Counting letters in phylip alignment columns with Biopython

I have been using python 3.6 and biopython 1.72 to work with protein data files. I am using a protein sequence file (phylip format), for example: ...
4
votes
1answer
73 views

Efficiently aligning a lot of reads on the same small reference sequence

The context: I have a DNA-sequence coding for a protein, about 1500 bp in length. Using NGS, a lot of reads of (mutants of) this same sequence were acquired. All of these reads need to be aligned to ...
1
vote
0answers
26 views

Stand alone chaining tool for existing blast hits?

I have a table of blast hits in a database that I'd like to chain. What is a good choice for a 'stand alone' chaining tool? Is it simple to implement in Perl, for example? Example 'hit' data is here. ...
2
votes
1answer
65 views

Minimap2 -ax map-pb doesn't output tlen field

I have used minimap2 to map some pacbio reads to a reference genome. I would like to know the "insert size" (true length of sequence) relative to the reference. More specifically, I want to know the ...
3
votes
0answers
70 views

Glocal\semi-local\Hybride Globale-Local alignment with Python

I was looking for a simple way to do a glocal alignment. The case I have is I have a small sequence which should be find in a bigger one, thus typically a glocal alignment. Also I can not Install ...
0
votes
1answer
81 views

Find a map/correspondence between two versions of a genome

I am working with two versions of the C. elegans genome. I am finding interesting regions (specifically, tRNA genes) in version 1 and then I would like to know if version 2 also has a tRNA gene in ...
2
votes
0answers
25 views

Alignment using secondary and tertiary features of DNA

I m trying to align different sequences of length not greater than 50. Is there any way to incorporate other information such as stacking energy, entropy, bonds etc. If there any way to align ...
2
votes
1answer
205 views

pysam pileup: what reads appear in the pileup?

Cross-posted on Biostars (with no answers currently). I hope that's OK. I thought that iff a read matches the reference at a position, it would appear in the pileup column for that position. But ...
2
votes
1answer
20 views

GRCm37-designed exon target enrichment, which reference to use?

I have exomes from 24 individual mice. The exomes are the product of Roche's SeqCap EZ HyperCap target enrichment kit. I see that the mouse exome design comes from mm9/NCBI37, the previous major ...
5
votes
2answers
1k views

Difference between samtools mark duplicates and samtools remove duplicates?

What is difference between samtools mark duplicates and remove duplicates ? Is it necessary to mark duplicates before removing duplicates with samtools?
3
votes
1answer
245 views

Increase number of threads for GATK 4.0 HaplotypeCaller

I am using GATK version 4.0, I tried to use multiple threads for calling variants using HaplotypeCaller using following command ...
8
votes
1answer
402 views

Better aligner than bowtie2?

Bowtie2 is probably the most widely used aligner because of it's speed. Burrow-wheeler (BW) algorithms (including bwa) tend to ...
1
vote
0answers
26 views

How can I determine a mean sequence divergence for 10k sequences?

I am trying to analyze a number of repetitive sequences and as one step want to calculate a sequence divergence between the elements I found. Now in theory I wanted to generate a MSA of the sequences ...
3
votes
1answer
141 views

Test to determine if two genes/exons share the same evolutionary histories?

In classic phylogenetic inference one is usually given various orthologue sequences of a given gene across various species. Those sequences are then multiple aligned and used to construct a ...
2
votes
1answer
83 views

Interpreting bayesian phylogenetic model using Tracer?

I am doing phylodynamic analysis of virus sequences using bayesian method using BEAUti and BEAST. As a part of my analysis I am using Tracer to check how my model is working. I can interpret some of ...
2
votes
0answers
51 views

How to create a nexus file for BEAST?

I am learning to use BEAUTI and BEAST software. I followed their tutorial and was able to run example data without any problem. Now, I want to do phylodynamic analysis on my data. I have downloaded 70 ...
8
votes
2answers
209 views

Is there a standard definition for “assembly polishing”?

Is there a standard definition for "assembly polishing" in the field? Is there a standard definition for what polishing algorithms do? My understanding of "polishing" is strongly influenced by ...
5
votes
2answers
247 views

Reads mapped to exonic, intronic and intergenic regions

After the alignment step I checked the rnaseq metrics of all the samples. Among 40 samples three samples show high percentage of reads mapped to intronic regions. What could be the reason? ...
7
votes
3answers
348 views

visualisation of genome alignment

I was asked to check the synteny of some genes in two genome assemblies of the same species (PacBio & Illumina). I was given two scaffolds couples (so 4 scaffolds total) in the Illumina genome and ...
1
vote
0answers
54 views

Mugsy error: Can't find species II dna at output/software/mugsy_x86-64-v1r2.3/mugsy line 501

I have installed mugsy in order to create a multiple genome alignment and a phylogenetic tree of several species of nematodes. The following command successfully pulls out the help: ...
5
votes
2answers
231 views

How are Principal Component analyses and Admixture analyses from a genetic alignment different?

How are Principal Component analyses and Admixture analyses from a genetic alignment different? My understanding is that a PCA will take raw genetic differences across the entire alignment and plot ...
2
votes
2answers
55 views

Question regarding the function pairwiseAlignment in R

I'm a bit confused and would appreciate your assistance, please. I'm trying to use the pairwiseAlignment function from the ...
5
votes
2answers
129 views

Is there a Python/R package with the ability to convert an alignment and reference into a CIGAR?

I'm writing a python function from scratch to do this, but I feel like this must exist in some standard bioinformatics library already. In principle, this is a simply regex operation which many must ...