Questions tagged [read-mapping]
When dealing with associating sequencing reads to a matching position in a set of reference sequences (typically a genome).
93
questions
0
votes
1
answer
123
views
Mapping List of Probes/Primers/Short Oligos to a Reference Fasta/q
Could you help me bioinformatics SE people.
So I have been duckduckgo-ing up a storm and even consulted the AI overlords and haven't come up with a solid way to do this, but it is so fundamental to ...
CommunityBot
- 1
2
votes
0
answers
54
views
Converting EntrezGene IDs to Ensemble IDs
I am trying to convert my gene_id accessions given by the RNAseq as EntrezGene gene numbers into Ensemble annotation. I need the Ensemble annotation for further GO analysis. I am using BiomaRT package ...
- 2,464
1
vote
0
answers
26
views
Is it normal to have low mapping to bacterial genome in a total RNA sample of plant root-colonized bacterial cells?
I had isolated the total RNA from a sample of plant roots with colonized bacteria by excising the roots (harboring the bacterial cells) and crushing them into a powder using liquid nitrogen. The RNA ...
- 149
1
vote
2
answers
77
views
How to know if FASTQ/BAM is from reference genome (FASTA)?
I'm new to bioinformatics. I have a problem in which I have a FASTA reference genome and lots of reads in FASTQ files. Some of them could be contaminants, so I'd like to filter them out and get only ...
- 2,464
0
votes
1
answer
264
views
Featurecount output .txt file from bam file
After running this command:
featureCounts -p -s 2 -a $genome -o $dir"/"$specie"/count_table/"$value".txt" $output_loc -T 8
on the BAM output of ...
- 1,510
4
votes
4
answers
656
views
what percentage of the human genome is MAPQ=0?
When doing Illumina 2x150bp sequencing of genomic DNA, and after aligning the reads to GRCh38, what percentage of the non-N fraction of the human genome is MAPQ=0? This is, what part corresponds to ...
- 10.5k
15
votes
2
answers
4k
views
How does the BWA-MEM algorithm assign its mapping qualities?
Is there any resource (paper, blogpost, Github gist, etc.) describing the BWA-MEM algorithm for assigning mapping qualities? I vaguely remember that I have somewhere seen a formula for SE reads, which ...
- 109
1
vote
2
answers
572
views
Parameter optimisation of minimap2
I'm leveraging minimap2 to map select genes from short-read microbial fastq metagenomics zip files.
...
M__♦
- 12.8k
3
votes
2
answers
180
views
Reference genome (bwa -mem)
I want to do mapping with bwa mem. I downloaded a reference genome from ncbi database which is (Danaus plexippus). I got one ncbi_database.zip file which contain 2 ...
- 103
1
vote
0
answers
201
views
Mapping statistics from the bam files
I would like to find the mapping statistics from the sorted bam files.
Samtools flagstat gives the output only for a single file. What is the easiest way to find the mapping statistics for all ...
- 11
4
votes
3
answers
113
views
mapping nucleotide sequences
I am trying to do mapping with multiple sequences against a reference genome. I tried the following code: The outputfile is .fastq.sam. I need only .sam. How do I get it?
...
- 6,575
3
votes
1
answer
56
views
How do I generate a variant list using Illumina reads from a Salmonella genome?
I am planning on performing phylogenetic analysis of Salmonella specimens using WGS data from PulseNet and GenomeTrakr, for the purpose of public health surveillance & to provide context for ...
CommunityBot
- 1
0
votes
1
answer
239
views
Question: Bismark methylation for non-uniquely mapping BS-seq reads
I am looking for the DNA methylation, specifically inside small RNAs; piRNA. I am mapping BS-seq reads on the genome with bismark. From the user guide of bismark, it maps reads uniquely on the genome, ...
- 1
4
votes
1
answer
355
views
Kallisto error: index input file could not be opened!
I am utilizing Kallisto in Anaconda/miniconda for RNA sequencing. I have successfully made ...
- 3,148
2
votes
1
answer
1k
views
Understanding PAF format
Does anyone have any experience with PAF files (minimap2 output)? I don’t fully understand what my output means.
I’ve looked at the minimap2 explanation (https://github.com/lh3/miniasm/blob/master/PAF....
- 14.9k
1
vote
1
answer
341
views
Nextflow with Minimap2 not working
I'm trying to create a Nextflow script that maps assemblies to reference gene sequences. This is what I have so far:
...
- 14.9k
2
votes
1
answer
123
views
Coordinate numbering between forward (+) and reverse (-) strand in Tombo
My question is related to the output of Tombo re-squiggling algorithm on nanopore data.
Given the position/coordinate on one strand (+ or -), I want to find the coordinate of its exact location on the ...
- 309
2
votes
0
answers
32
views
MAGs transcriptomics pipeline question
I have currently a following problem. I have a one sample that I did my metagenomics on (Illumina shotgun + nanopore) and have recovered some high quality MAGs. I also did a metatranscriptomics on the ...
- 118
1
vote
1
answer
155
views
Low mapping reads for RNA sequence
I am new to RNA sequencing analysis. I have RNA samples for bacteria for different treatment. I did contaminant filtration, remove adapters and checked fastQC report: attached below for one sample (...
CommunityBot
- 1
2
votes
1
answer
1k
views
how to retain reads with low mapping quality (MAPQ) scores when using samtools view -q
I have been using the -q option of samtools view to filter out reads whose mapping quality (MAPQ) scores are below a given ...
- 4,212
2
votes
1
answer
253
views
Get number of reads with a single, (almost) exact match to the full length of a reference sequence
I can't find an answer to this in previous questions, so hoping someone can help me now.
We have nanopore sequenced a PCR product, and I've filtered our reads to +/- 10 bp of the expected product ...
- 14.9k
2
votes
1
answer
73
views
What tool I can use to map short-reads sequences to reference genome and get specific mapped size
I have about 300 90-bp sequences which I would like to map to a reference genome to make it longer to a 300bp sequence, wherein the 90bp is at the middle.
Anyone knows what bioinformatics tool I can ...
- 166
9
votes
2
answers
3k
views
How do I generate a variant list (i.e. VCF file) using Illumina reads from a human genome?
This is a problem I have to solve frequently, and I'd be interested in knowing what other methods people use to solve the same problem.
About twice a year, I get asked to determine variants from ...
- 14.9k
3
votes
0
answers
81
views
Patient-sample mapping in GSE72056 dataset
I want to use the single-cell data from the following expression profiling which concerns the RNA-seq of metastatic melanoma:
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE72056
The data was ...
M__♦
- 12.8k
5
votes
1
answer
437
views
Displaying soft-clipped nucleotides in samtools tview
There are nicer genomics visualization tools available, but the samtools tview command is almost always my go-to for a quick first look at read alignments. I just ...
CommunityBot
- 1
1
vote
1
answer
236
views
What is the meaning of split read?
I want to use rna seq data to later perform functional tests on fusion genes.
so before that I need to filter the "best results" (of rnaseq) for deciding which candidates I actually want to ...
CommunityBot
- 1
3
votes
1
answer
51
views
RNASeq read coverage in protein space?
I used samtools depth to find the per base-pair read coverage over a number of isoform contigs from my Trinity assembly.
I have also conducted a multiple sequence ...
- 14.9k
1
vote
0
answers
228
views
using Bowtie to map miRNA-seq data to rfam and also reference genome
I am using Bowtie to remove non-coding RNAs (tRNA, snRNA, rRNA) by rfam and also maping our microRNA-seq data with mirbase using following code:
...
- 4,202
1
vote
2
answers
89
views
Solutions to reference bias issue
I've recently learned about the reference bias issue - inability to properly map NGS reads in certain genomic regions caused by the fact that our standard genomic references are linear and do not ...
- 143
0
votes
1
answer
57
views
Modeling number of reads mapped to a gene
I am looking for a probability distribution of a number of reads mapped to a particular gene in metagenomic sequencing (NGS, shotgun, likely illumina).
Naively one could model it via a binomial (or ...
- 14.9k
3
votes
1
answer
511
views
How to check what software has been used to generate bam file?
I have several bam files from unknown origin. I'd like to know what software was used to obtain those bam files, e.g. bwa-mem. Is it possible to somehow check it?
- 841
2
votes
3
answers
221
views
Quantify low complexity regions in DNA sequences
I have fasta files with multiple sequences. They are reads that mapped to the genome outside of probe-targeted regions.
From a quick perusal, they appear to be repetitive and have low complexity. Is ...
- 957
1
vote
0
answers
328
views
How to rename Read Groups
I periodically run into problems with Read Groups. I am in a desperate need for a tool that would allow me to fix them.
This time, when I was mapping reads, I formatted all the information in the Read ...
- 5,632
0
votes
1
answer
307
views
Reads count in metagenomics
Background: I am developing a pipeline for metagenomic studies of human gut microbiote. In particular, I am mapping the reads data originating from shotgun whole genome sequencing onto a gene ...
- 4,212
1
vote
1
answer
2k
views
What is "unmapped read segments" in the output of samtools idxstats?
samtools idxstats produces a four column output (see here)
...
- 19.8k
1
vote
0
answers
48
views
How to control/normalize for number of reads when calling SNPs using RNA-Seq?
I used the GATK pipeline to call SNPs on males and females using RNA-Seq data. But the males have a higher read count (~43-46M reads) than the females (~40-42M reads). This causes SNP counts to be ...
- 86
3
votes
2
answers
292
views
Contamination on genome assembly
I had a question for the community.
I have a genome of a new species that has been sequenced via 150pb Illumina paired-end.
To verify the quality of the assembly I used the ...
- 5,632
1
vote
1
answer
1k
views
Difference between paired-end, mate-pair and long read
I writing here because I have some questions for you.
I wondered what the essential differences were between paired-end, ...
2
votes
0
answers
227
views
How snippy makes MSA-like aligned fasta output from multiple samples?
From the log file it seems snippy doesn't do assembly. It only does mapping:
fastq --> SAM --> BAM --> VCF --> consensus_seq/snps
But if multiple ...
M__♦
- 12.8k
2
votes
1
answer
57
views
Assemble reads for a specific gene
I have a lot of unassembled sequencing data and I want to build phylogeny for some genes from these data. I can retrieved the sequene by mapping reads to a reference sequence but if the two species ...
- 127
2
votes
4
answers
623
views
Is there a read mapper that allows input in SAM format instead of FASTA/FASTQ?
When we get raw sequenced reads in FASTQ format, we usually need to do some pre-processing (QC, demultiplexing, etc.) prior to sequence alignment. The only way to register the results of this ...
15
votes
2
answers
14k
views
Meaning of BWA-MEM MAPQ scores
Does anyone know what the MAPQ values produced by BWA-MEM mean?
I'm looking for something similar to what Keith Bradnam ...
- 14.9k
1
vote
1
answer
167
views
How do you map nanopore fast5 files?
Is there a best practices for mapping Oxford Nanopore files to a reference?
Is there a tool that can take a tarball of fast5 files and map them directly or do they need to be converted to fastq first?...
- 14.9k
1
vote
0
answers
25
views
How to choose reference to calculate coverage of WGS of organism that has high intra-species variation?
Calculation of coverage involves mapping reads with a reference and depending on the choice of reference we can get very different result. If a species has a genome that has very low size variation ...
- 367
2
votes
2
answers
2k
views
Error with BWA Mem input having multiple fastq files using cat and process substitution
Running BWA Mem on a number of paired end fastq files using process substitution on the inputs results in this error:
...
CommunityBot
- 1
2
votes
1
answer
676
views
What do the symbols mean in minimap2's gap cost equation?
Heng Li's paper on minimap2 is understandable to me right up until the moment that it gives this definition of gap cost:
2.2.1 Alignment with 2-piece affine gap cost
Minimap2 performs DP-based ...
- 155
0
votes
0
answers
51
views
How to obtain the extended unaligned bases of the read in reference based genome assembly?
I did illumina read mapping based on a particular gene sequence (reference sequence) using the following command,
...
- 109
2
votes
1
answer
854
views
Oxford Nanopore mapping Quality and sequencing error
I'm analyzing some Oxford Nanopore sequencing reads, with several articles stating the sequencing error is between 5% and 15%. Since its 1D and not 1D^2 sequencing i would assume the sequencing error ...
- 19.8k
4
votes
1
answer
255
views
Problem: "Pair-end" reads scRNA seq data (Drop-seq)
In case of Drop-seq, we have paired end data.
Read 1: Cell code + UMI (unique molecule identifier)
Read 2: The transcript information
But I have a problem/doubt with the sample I am working on.
...
- 31
2
votes
2
answers
491
views
How can I transform a mapped BAM file into an unmapped BAM file?
In order to use MergeBamAlignment (Picard), I need an unmapped BAM file and a mapped BAM file. I have two mapped BAM files:
one with reads mapped to the reference I want but without metadata such as ...
- 1,971