Questions tagged [sequencing]

DNA sequencing is the process of determining the precise order of nucleotides within a DNA molecule.

Filter by
Sorted by
Tagged with
0 votes
1 answer
49 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 ...
user avatar
0 votes
0 answers
10 views

Interpreting 'samtools mpileup' output for multiple inputs

I would like to calculate sequencing coverage for a WGS project. Both long and short reads. I've used samtools as following: ...
user avatar
  • 539
0 votes
1 answer
34 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 ...
user avatar
  • 127
0 votes
0 answers
20 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?
user avatar
1 vote
1 answer
36 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 ...
user avatar
  • 19
1 vote
2 answers
56 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 ...
user avatar
0 votes
0 answers
26 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 ...
user avatar
1 vote
1 answer
77 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 ...
user avatar
  • 594
0 votes
1 answer
61 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 ...
user avatar
  • 23
0 votes
1 answer
67 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.
user avatar
  • 151
1 vote
1 answer
48 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 ...
user avatar
  • 151
2 votes
1 answer
76 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 ...
user avatar
  • 21
1 vote
0 answers
159 views

Illumina vs PacBio vs Nanopore price

Can someone help me with comparing the cost per run of Illumina, PacBio and Nanopore sequencing technologies?
user avatar
0 votes
1 answer
40 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.
user avatar
0 votes
1 answer
49 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
user avatar
  • 37
2 votes
2 answers
219 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 ...
user avatar
  • 21
2 votes
1 answer
178 views

What is the right way of calculating a Phred score by hand?

i am trying to calculate mean Phred scores for my sequencing data, but i feel not very comfortable about it. There are actually two ways of calculating. (I just use an existing sample) giving: 3 reads ...
user avatar
  • 35
1 vote
2 answers
86 views

Get frequency of index sequences of fastq file

I recently did some sequencing on a MiniSeq and unexpectedly the Undetermined_S0_L001_R1_001.fastq.gz file contained quite a lot of data (about 40% of total). I ...
user avatar
0 votes
1 answer
30 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 ...
user avatar
-1 votes
2 answers
90 views

Filter rows of VCF file for Match=EXACT?

How do I remove rows in a VCF file on 1kGenome column, where Match=EXACT using bash ?
user avatar
  • 11
0 votes
1 answer
63 views

What is a PacBio "movie file"?

I came across references to "movie files" from PacBio sequencing in this paper: https://www.jimmunol.org/content/204/12/3434 Specifically: Movie files used to generate results presented in ...
user avatar
  • 594
0 votes
1 answer
52 views

How to determine a correct promoter sequence? (Bisulfite sequencing)

I'd like to ask you for help with an university assignment. First let me post the entire question: Design primers for bisulfite sequencing of the promoter of the ACE2 gene in microbat Myotis ...
user avatar
0 votes
1 answer
53 views

Quast duplication ratio and mismatches percent

When analyzing Quast results it seems that it doesn't calculate mismatches and indels in a useful way if the "Duplication ratio" is over 1. For example, that's what I get for an assembly ...
user avatar
6 votes
1 answer
75 views

How to search for high coverage SRA entries

The Question I want to find high coverage SRA entries, e.g., above 100x. I guess the best way is to use https://www.ncbi.nlm.nih.gov/sra with an appropriate search term. I don't mind if the search ...
user avatar
1 vote
1 answer
79 views

Metagenome simulation from a concatenated FASTA file

I am trying to simulate metagenome sequencing. So I will start with a file with a lot of concatenated genomes. From there, I would like to randomly extract 10,000 sequences of length 200bps. I got ...
user avatar
  • 41
1 vote
1 answer
81 views

chipseeker annotation issue

...
user avatar
  • 1,538
0 votes
0 answers
24 views

DNA sequence error annotation

Suppose I generate a DNA sequence of the following pattern. AAGTC And after being passed through a channel with insertion, deletion and substitution errors obtain the following sequence AAAGGC Where ...
user avatar
2 votes
2 answers
147 views

How can I subset WGS data to the level of WES variants?

I would like to compare mutational signatures1 in patients from different studies, however some studies are based on exome seq (i.e. ~20,000 coding variants) and some are from whole genome seq (i.e. ~...
user avatar
2 votes
2 answers
163 views

Are sequencing reads written differently when inserts and indexes are sequenced seperately?

I was going over different sequencing library protocols and noticed that the adapter sequence can vary between protocols. In some adapter sequences the index sits between the primer binding site and ...
user avatar
  • 21
1 vote
0 answers
27 views

Fastq dataset after clustering

I'm doing a project for DNA archival storage systems Is there an open-source dataset used as ground truth for clustering? E.g A file containing the ground truth strand sequences, which have been ...
user avatar
3 votes
2 answers
90 views

Virus genome: protein inside another protein

The genome sequence for the GenBank virus JX869059.2 (a human betacoronavirus) contains, among others, two proteins: N protein (range 28566..29807) and Orf8b (range 28762..29100). The interval for ...
user avatar
  • 215
1 vote
2 answers
92 views

Bacterial DNA at the tail of transcriptome reads. What does that mean?

I am assembling a transcriptome obtained from the Internet. The transcriptome was extracted from a human cancer tissue that had been previously grafted into a mouse. I have detected that many ...
user avatar
0 votes
1 answer
101 views

How to anonymise a bam file

I have a bam file of, for instance, RNA-Seq, which contains patient-identifiable data in the form of Single Nucleotide Polymorphisms (SNPs) throughout. I would like a technique to take this aligned ...
user avatar
  • 121
0 votes
0 answers
35 views

How many markers or loci are used in forensic DNA profiling and the effect of the DNA sequencing?

I am reading the Wikipedia page on DNA profiling, specifically the STR analysis. It says the DNA profiling relies on 20 or so loci matching. It also gives an account of finding false identification. ...
user avatar
  • 179
2 votes
4 answers
116 views

Separate multiple sequence into different file, one sequence per file

I have a file contain multiple sequence, and I want to separate them by "gene:" into different file. example: example.fa ...
user avatar
  • 33
1 vote
2 answers
108 views

How to combine two Genome-wide Association Study (GWAS)? [closed]

I did a GWAS analysis in the past for antibiotic resistance of E. Coli and the results were interesting (matching the literature). I did a new GWAS analysis for some new samples, but the results are ...
user avatar
  • 107
1 vote
1 answer
90 views

merge the matrix from computematrix of deeptools

I have several matrix from computematrix of deeptools. I need to merge two of them using "computeMatrixOperations cbind -m input1.mat.gz input2.mat.gz -o output.mat.gz" but I am running to error "/...
user avatar
  • 95
2 votes
1 answer
61 views

Assembly reads having a copy of their beginning in their tail

I am analyzing the reads for the SARS-CoV-2 assembly with id SRR11140748. Apparently these reads were obtained with parallel sequencing by Illumina and Oxford Nanopore Technologies. I have found these ...
user avatar
2 votes
1 answer
135 views

What is the meaning of these misaligned reads in a sequencing run?

I am analyzing some SARS-CoV-2 sequencing runs abd often find read alignments like the one in the image. ...
user avatar
1 vote
2 answers
757 views

Marking optical or PCR duplicates with picard vs. samtools flagstat

I am trying to identify the best library prep method for noninvasive prenatal test samples, to be sequenced with Illumina Novaseq. One metric that I am evaluating is the number/ percentage of PCR or ...
user avatar
  • 183
3 votes
1 answer
79 views

Overlap in Paired-end Reads for Sequencing?

Regarding Sequencing: Do/can the paired-end reads have overlaps?
user avatar
  • 107
1 vote
2 answers
182 views

Can I automate CLUSTALO and output alignment sequence identity?

I've detected homology between targets of ligands in drugbank and proteins in the proteome of a pathogen. I've parsed the output very rudimentary and calculated my query coverage. This exists in an ...
user avatar
  • 21
1 vote
1 answer
53 views

Location of Reads in Sequencing

I have questions regarding sequencing: Are paired-end reads from different stands? What about single reads, are they coming from both strands?
user avatar
  • 107
1 vote
1 answer
48 views

Warning in fastqc

I am checking the quality control of my sequences using the Fastqc tool. For some steps, like "Per base GC content", etc., I received a warning. So, I was wondering whether I Should also take care of ...
user avatar
  • 107
1 vote
1 answer
228 views

Per Base Sequence Content in fastqc

I have a question regarding "Per Base Sequence Content" plot for "fastqc": In the fastqc documentation, it is written: "In a random library you would expect that there would be little to no ...
user avatar
  • 107
1 vote
1 answer
32 views

Cutting the sequence into several sequences with the information of a dataframe

I have a fasta file with several 120-concatenation protein sequences. I also have a data frame with 120 names of proteins in column 1 and their length in column 2. Using this dataframe I want to ...
user avatar
2 votes
1 answer
66 views

Translating a genome sequence into possible frames

I have a sequence below from MYC gene about which I need to translate the sequence in all possible frames AND identify (for each frame) which codons are actually used in MYC ...
user avatar
0 votes
1 answer
299 views

What samples can be used for a Mutect2 Panel of Normals

I'm working on calling somatic variants from solid tissue tumors. I have one normal and one tumor sample for each patient. I plan to use Mutect2 to call somatic variants after preprocessing the data ...
user avatar
2 votes
2 answers
560 views

Importance of Proper Pairs vs Aligned Reads for RNASeq data

I have stranded, paired end RNASeq reads that I have aligned using STAR. I plan do conduct a differential expression analysis with DESeq2. After running quality control checks, a good portion of my ...
user avatar
33 votes
4 answers
9k views

Why does the FASTA sequence for coronavirus look like DNA, not RNA?

I'm looking at a genome sequence for 2019-nCoV on NCBI. The FASTA sequence looks like this: ...
user avatar