Questions tagged [alignment]
These questions are about sequence alignment.
2
votes
2answers
29 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
25 views
2
votes
1answer
50 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 ...
2
votes
0answers
25 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
126 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
43 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
97 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
73 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
81 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
69 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
12 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
45 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
168 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
48 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
62 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
53 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
44 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
232 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 ...
3
votes
1answer
79 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
81 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
69 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
47 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
63 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
68 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
24 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 ...
0
votes
1answer
64 views
GMAJ error: Exception in thread “main” java.lang.ClassCastException: java.lang.StringBuffer cannot be cast to java.lang.String
I have installed GMAJ using:
...
2
votes
1answer
135 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
17 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
601 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
154 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
262 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
24 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
139 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 ...
1
vote
1answer
52 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
37 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
146 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
186 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?
...
6
votes
3answers
229 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
46 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
151 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
51 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
122 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 ...
5
votes
2answers
175 views
Books on bioinformatics algorithms
I'm looking for a book about bioinformatics algorithms, such as alignment, BLAST search, and variant calling.
I'm hoping reading about this subject will give me a deeper understanding of the ...
5
votes
1answer
78 views
Are there certain alignment methods/tools which perform better with a high density of indels?
I have a set of experiments which should result in WGS reads with a high density of indels.
Question: Are there certain alignment tools/methods which perform respectively "better" with these reads? ...
7
votes
1answer
252 views
Extracting strand specific reads from MinION cDNA-PCR protocol
I recently performed my first MinION run, and now I'm trying to analyze the data. Being pretty new to the bioinformatics field, I was hoping some of you could help me out.
As a bit of background, I'...
1
vote
0answers
67 views
igv_plotter error: Xvfb did not start
I am trying to get screenshots from alignments of reads to a genome in IGV (Integrated Genome Viewer) using igv_plotter.
...
6
votes
1answer
75 views
Does the DNA or RNA of some ethnic groups map better than others' to the human reference genome sequence?
I believe that the human genome reference took DNA samples from different people and that some natural variation is included in extra contigs.
However, the human reference genome comes from a ...
2
votes
0answers
63 views
aligner for 1D^2 oxford nanopore data
I am wondering if anyone has tried aligning 1d^2 RNA-seq data. I am trying that on minimap2. It is giving me very low ~50% of mapped reads. it gives ~95% mapped reads for 1D.
I am using the following ...
5
votes
3answers
193 views
Generating the reconstructed alignment from BAM
I have a (small) BAM file with CIGAR and MD fields.
Question 1: What tools exists in Python and/or R to reconstruct the alignment between the reference and the read in a BAM? Given that this is a ...
