Questions tagged [ngs]
use for general applications of high-throughput nucleotide sequencing methods that use short-read technology (e.g. Illumina, IonTorrent). For long-read sequencing, use 'long-read', or a more specific tag if applicable (e.g. nanopore, pacbio).
257
questions
2
votes
0
answers
13
views
Off-target % for whole-exome sequencing panel
My samples have been sequenced with the Twist Exome 2.0 sequencing panel, and I want to assess how efficient the sequencing has been. By efficient, I mean examining metrics such as the uniformity (...
- 171
2
votes
1
answer
50
views
calculate mismatch frequency/rate from a BAM file
I am working with PAR-CLIP data where experimental procedure induces T to C mismatches. The other types of mismatches we can consider coming from artifact such as PCR errors. It could also be an ...
- 31
-1
votes
1
answer
52
views
Variants list to aminoacids
I have the sequencing data of a Monkeypox virus sample. I have the list of variants extracted from Nextclade and I want to find the corresponding aminoacids substitution.
For example I have the ...
0
votes
1
answer
40
views
How do I build the MinusB database for Kraken2? (Taxonomy issues)
I am attempting to build my own custom database for Kraken2. I have two questions:
If I have the MCPyV genome in a file called MCPyV.fasta, how do I build a database with just this?
How do I build ...
0
votes
0
answers
17
views
Kraken2 Standard Database failing to build (unexpected FTP path)
I am attempting to build just the standard Kraken2 Database by using the following command:
kraken2-build --standard --threads 24 --db $DBNAME
I am returned the ...
1
vote
1
answer
126
views
How do I build the "Standard Kraken 2 Database"?
I am at this point in the GitHub tutorial: https://github.com/DerrickWood/kraken2/blob/master/docs/MANUAL.markdown#standard-kraken-2-database
It says to create the standard Kraken database, I use the ...
1
vote
0
answers
15
views
DNASTAR viral-host integration assembly keeps failing
I have two NGS files from an NGS company corresponding to the sequencing data from a tumor sample as follows:
TB_7710391_R1.FASTQ.gz
TB_7710391_R2.FASTQ.gz
I have downloaded the genome for MCPyV as ...
1
vote
1
answer
137
views
Estimate insert size for single-end data with picard CollectInsertSizeMetrics
I have a BAM file generated from single-end data and I want to estimate the insert size using picard's CollectInsertSizeMetrics as follows:
...
- 323
2
votes
1
answer
61
views
Sanger sequencing annotation error
I am a student in a Cancer lab. Working with sanger is new to me. While analyzing a report we found an insertion that has not been reported in any databases so far, we were working on checking if the ...
- 21
1
vote
1
answer
40
views
Problems with BED format for ClinCNV
I have some problems with ClinCNV(https://github.com/imgag/ClinCNV):
I made a BED file from a BAM, then i made a GC count with ngs-bits (link: https://github.com/imgag/ngs-bits )
The main idea of ...
- 15
1
vote
1
answer
288
views
Need to convert .bed to .vcf. Can the reference build (and necessary .fasta) be determined for a .bed file?
I have some .bed files (and .bim and .fam files) containing data for a number of different samples, and I need to convert them to .vcf.
I've found bed2vcf which is ...
2
votes
0
answers
78
views
BAF and LRR calculation and usage
I'm working with a vcf file, in particular I'm trying to identify a gene duplication in a sequenced genome.
I studied bcftools cnv command, but for this tool is needed a vcf file containing a BAF and ...
- 21
1
vote
1
answer
73
views
Is it common to get different number of SNVs+Indels across samples from vcf files generated using GATK and DRAGEN (counts are higher for GATK)?
We have 2 vcf (Whole Exome Sequencing (WES) data; germline samples) files (e.g., vcf_1 and vcf_2). vcf_1 was generated (Ref. genome: hg38) using the GATK pipeline for 250 children and their parents ...
- 11
2
votes
0
answers
77
views
Adapter trimming
I am trying to do adapter trimming, alignment and sorting for a range of large scale paired end fastq files. The code I am using is given below:
...
- 143
2
votes
0
answers
27
views
What data set was used Gene Expression Data of Postmortem Tissues?
I want to run the experiments "Sample Application to Genotype-Tissue Expression (GTEx) Project Gene Expression Data of Postmortem Tissues" mentioned in "Enter the Matrix: Factorization ...
- 487
1
vote
1
answer
58
views
Expression analysis of miRNAs with normal RNA-seq data without small RNA-seq data
I am looking to perform expression analysis of miRNAs with normal RNA-seq data lacking small RNA-seq data?
Which path should I choose for known miRNAs and unknown new miRNAs?
Data set: rna-seq data ...
1
vote
1
answer
51
views
Lower case vs. upper case nucleotids in sequence vs. dots at the end
What is the difference between lower case and upper case nucleotides in a sequence? My other question is what are the dots at the end of the sequence? Some examples are shown below:
GGgG,GGgG,GGgG,...
3
votes
0
answers
32
views
Can multiplexing in Sequel II SMRTcells reduce the coverage?
I would like to have high depth of coverage of the transcriptome, but I need to analyze several tissues. I was suggested to pool all 5 tissue samples in same SMRT 8M cell. Will this result in a ...
- 257
1
vote
1
answer
37
views
How can Iso-Seq reverse transcriptase artifacts be avoided?
My end goal is annotate a de-novo assembled genome. When trying to select the best method for transcriptome assembly I read Iso-seq was the preferred method. However other people suggested Nanopore as ...
- 257
2
votes
2
answers
121
views
What is the most accurate approach for de Novo sequencing?
I'm trying to decide between PacBio HiFi or Illumina sequencing platforms for sequencing the genome of aChrysina scarab. We want to identify the color pattern loci and also perform a phylogenetic ...
- 257
2
votes
1
answer
127
views
Count reads at specific gene features
I have BAM files from an RNA-Seq experiment and for all genes (or a subset) I want to get the number of reads in regions around the TSS (e.g. 2kbp) and the TES (e.g. 2 kbp) and calculate the ratio ...
- 323
1
vote
0
answers
36
views
System specifications for NGS data analysis
I have a 3.7 GB whole genome data of a eukaryote, for which genome assembly, gene prediction, and annotation steps have to be performed. In some time I would also need to analyze the transcriptome ...
- 11
2
votes
3
answers
87
views
Can we do NGS library preparation using UMI with large amount of DNA input?
We know that in next-generation sequencing (NGS), the unique molecular identifier (UMI) can reduce or eliminate sequencing or PCR errors and result in very high accurate data. Therefore, UMI is widely ...
- 21
1
vote
1
answer
149
views
How can I generate ASV file from nanopore sequencing data?
I am converting the fast5 file to fastq by using guppy basecaller after that by using kraken2 classified the sequence.
Now I am trying to generate ASV file.
Is this possible to generate ASV file from ...
- 31
1
vote
1
answer
308
views
Extracting base sequences from ABI/AB1 sanger sequencing chromatogram
I am trying to understand the sanger sequencing ABI/AB1 file format better, and extract base calls from given signal intensities over time.
As I understand, reading in a raw AB1/ABI file into python, ...
- 11
2
votes
2
answers
117
views
BWA-mem and sambamba read group line error
this question has been asked [and answered] on Stack Overflow
This is a two-part question:
help interpreting an error;
help with coding.
I'm trying to run bwa-mem ...
1
vote
2
answers
89
views
How can I get latest .vcf files with annotation data?
I'd like to perform the annotation (1000genomes, COSMIC etc.) of my variants using SnpSift and SnpEff, however so far all I get are vcf files for separate chromosomes:
http://ftp.1000genomes.ebi.ac.uk/...
- 206
4
votes
2
answers
122
views
Does the number of RNA reads per cell obtained from the 10X scRNA experiment depend on amount of mRNA in given cell?
As we know, the amount of RNA reads per cell obtained from 10X scRNA experiment vary between cells. I wonder if this is effect of technical issues or does the number of RNA reads per cell obtained ...
- 43
0
votes
1
answer
137
views
Why must a maximal non-branching path be a contig?
The following is from Bioinformatics Algorithms:
Fortunately, we can derive contigs from the de Bruijn graph. A path in a graph is called non-branching if in(v) = out(v) = 1 for each intermediate ...
- 127
1
vote
1
answer
51
views
Reference variant detected as altered one in bam file
I received (from manufacturer) several .bam files and I used four callers (samtools, freebayes, haplotypecaller, deepvariant) to find some sequence variants. In obtained .vcf files, I took a closer ...
- 206
0
votes
0
answers
73
views
Locate the saved template of minknow
I am using minknow on ubuntu. After saving the template (settings), where is that file locally stored in the system?
Also what is its format/extension?
- 1
3
votes
0
answers
52
views
A question on Homer normalization
When using annotatePeaks.pl script from Homer software to create histograms, the output is normalized per bp per peak (on top of normalization to 10 million tags).
What does it mean to normalize per ...
- 177
3
votes
1
answer
44
views
Metagenomics pipeline recommendations for short-read data
de novo metagenomics on viral NGS data is a hot-topic. On this site alone at least 4 specific algorithms have been used to identify multi-strain/multi-species for a given data set, however these do ...
M__♦
- 10.2k
1
vote
1
answer
49
views
What are Dummy fastqs?
I had to do a file transfer for in our system to run from fastqs and I came across the term "dummy fastq", not sure what exactly it means and what is the purpose of it in the workflow. Can ...
- 19
1
vote
2
answers
293
views
How to retrieve fasta sequence after local blast?
I have created a Blast database using a reference genome. Then, I have performed a local blast search in command line using a gene of interest. I have obtained some hits with the usual Blasting ...
- 75
0
votes
0
answers
28
views
Drug-DNA complex simulation using AMBER
I want to carry out the simulation of drug-dna complex by placing the drug molecule at particular site in my DNA sequence (intercalating site). I have PDB files of both drug and DNA sequence. I will ...
2
votes
1
answer
170
views
When using cutadapt, can I specify an R2 adapter as optional when I have a required R1 adapter?
I'm using cutadapt 3.5 to trim adapters and perform some filtering on paired-end data.
Both R1 and R2 sequences have 3' adapters that might be found depending on the sequence length, but R1 also has ...
- 737
0
votes
1
answer
88
views
Shotgun analysis: Where and how to start?
I am a complete rookie in shotgun analysis, having only performed marker gene analysis before. I have a pipeline that I want to try but I can't wrap my head around how it all works.
Please correct me ...
- 108
1
vote
1
answer
186
views
Dictionary and index of vcf for base recalibration step
I need suggestions in creating index and dictionary with vcf files. For the base recalibration step, I downloaded Homo_sapiens_assembly38.known_indels.vcf.gz from the given link: https://console.cloud....
- 11
0
votes
0
answers
99
views
Why is bcl2fastq2 taking so long to calculate stats?
Our lab has been using bcl2fastq v2.20.0.422 to demultiplex RNA-seq data sequenced on an Illumina Novaseq machine on a beefy EC2 instance and we've run into the strange problem: namely that while ...
- 101
0
votes
0
answers
272
views
BWA-MEM2 alignment-Snakemake
I have started using snakemake 6.5.2 to align fastq files with reference file. I have pasted the error below in this question. How to allocate memory in the snakefile and read the header from samfile, ...
- 11
0
votes
1
answer
156
views
How to Read SCF file in Python?
Is there any way that I can read SCF file in python like in R using sangerseqR, I have tried with Biopython, it seems it does not support this format.
- 151
1
vote
1
answer
79
views
Where to find publicly available Sanger chromatography data?
I am looking for any public available databases where I can download Sanger chromatography data ideally in Ab1 or SCF file format. I need extremely large amounts of this, since I want to use it for ...
- 151
2
votes
1
answer
83
views
Low pass sequencing has been reported to detect common variants. How low can one go and get reliable data? Is 2X pass sequencing analysis possible?
I would like to use low pass sequencing to replace a genotyping chip to be able to detect variants up to 0.1 % allele frequency in available population data. What is the minimum depth I can opt for to ...
- 21
1
vote
1
answer
50
views
Output of allelic association doesn't write the rsID [closed]
I'm using plink (1.9b5) to do allelic association. My problem is that my output does not write the rsID (SNP), it just writes a dot.
Output:
...
- 11
1
vote
0
answers
251
views
Illumina vs PacBio vs Nanopore price
Can someone help me with comparing the cost per run of Illumina, PacBio and Nanopore sequencing technologies?
- 31
0
votes
1
answer
63
views
Which format is the most raw format
I wanted information regarding which format is considered the most raw from an illumina sequencer ?
fasta
fastq
bcl
bam
As per my research it should be bcl but I am not sure.
- 11
0
votes
1
answer
51
views
Which sequencing technologies are considered short read
I couldn't find any information on this online. Which of the following sequencing technologies are considered short read technologies? Illumina, Oxford Nanopore, Ion Torrent, Pacific Biosciences
- 19
1
vote
1
answer
93
views
RNA_Seq data aligned used uniquely or multi mapped reads impact on result interpretation
I have some transcriptomic (Whole) sequencing data that I should analyse. I would like to do raw data alignment to a reference genome taking into account the multi mapped reads and uniquely mapped ...
- 151
2
votes
2
answers
397
views
Construct the Overlap Graph of a Collection of k-mers
I am doing a course on graph algorithms in genome sequencing and the current assignment involves building an overlap graph from a bunch of k-mers. Formally the problem is as follows:
Input: A ...
- 21