Tag: PHRED

Issue with fastq after converting phred 64 to phred 33 quality scores

Hello, I ran seqtk seq -VQ64 read1.fastq.gz > read1_phred33.fastq to convert my 64 based phred score reads to 33 based phred score phred reads. However when I attempted to run them through tophat alignment I got this error: Saw ASCII character 4 but expected 33-based Phred qual. terminate called after…

Continue Reading Issue with fastq after converting phred 64 to phred 33 quality scores

plotting roh from bcftools

plotting roh from bcftools 0 Heys, I am following this small tutorial on how to calculate ROHs from a vcf file using bcftools (samtools.github.io/bcftools/howtos/roh-calling.html) and I am getting this txt file: # This file was produced by: bcftools roh(1.10.2+htslib-1.10.2-3) # The command line was: bcftools roh -G30 –AF-dflt 0.4 my_file.vcf…

Continue Reading plotting roh from bcftools

How can I get PHRED score?

How can I get PHRED score? 1 Hi, all. I am trying to get the assembly stat(Table S1.) according to the following paper about de novo assembly. [www.ncbi.nlm.nih.gov/pmc/articles/PMC7266049/%5D%5B1] In the table, there is an item “Mean read PHRED score after filtering and trimming”. How can I get this? Is there…

Continue Reading How can I get PHRED score?

The sardine run in southeastern Africa is a mass migration into an ecological trap

INTRODUCTION Large-scale annual migrations occur in an extraordinary range of animals, from insects to the great whales. While the driving mechanisms of these migrations are varied and sometimes poorly understood, they often represent a way of optimizing conditions for breeding and adult fitness when these are in conflict. Often, populations…

Continue Reading The sardine run in southeastern Africa is a mass migration into an ecological trap

Trimmomatic error

Trimmomatic error 1 Hi everyone. I’m trying to trim some read data but i’m getting an error message. This is my input: trimmomatic PE -threads 24 -phred 33 /home/tbeckett/lustre/practice/output_data/ Filtered2S1_L3_R1.fastq.gz /home/tbeckett/lustre/practice/output_data/ Filtered2S1_L3_R2.fastq.gz /home/tbeckett/lustre/practice/output_data/trimmed/ TrimmedFiltered2S1_L3_R1_p.fastq /home/tbeckett/lustre/practice/output_data/trimmed/ TrimmedFiltered2S1_L3_R1_un.fastq /home/tbeckett/lustre/practice/output_data/trimmed/ TrimmedFiltered2S1_L3_R2_p.fastq /home/tbeckett/lustre/practice/output_data/trimmed/ TrimmedFiltered2S1_L3_R2_un.fastq ILLUMINACLIP:NexteraPE-PE.fa LEADING:20 TRAILING:20 MINLEN:60 This is the error i’m getting:…

Continue Reading Trimmomatic error

Illumina Q score

Illumina Q score 1 Hi all, I have Illumina sequencing results of a bacterial genome and a quality score of 35.89 is associated with these data. I know that a quality score of 30 is 99.99% of base calling accuracy based on this but what about the meaning of 35.89?…

Continue Reading Illumina Q score

Oncogene Concatenated Enriched Amplicon Nanopore Sequencing for rapid, accurate, and affordable somatic mutation detection | Genome Biology

Stochastic Amplicon Ligation. DNA samples for oncology sequencing are typically extracted from FFPE tissues and can have average lengths of less than 500 nt due to accumulated chemical damage [18]. We developed the Stochastic Amplicon Ligation (SAL) method to enzymatically concatenate many short DNA molecules together to utilize the long-read…

Continue Reading Oncogene Concatenated Enriched Amplicon Nanopore Sequencing for rapid, accurate, and affordable somatic mutation detection | Genome Biology

Rsubread align maximum nthreads

Hi Experts, I am using Rsubread align using following comand- align (index=”my_index”, readfile1 = “SRR123456_1.fastq” ,readfile2= “SRR123456_2.fastq”, type=”rna”,input_format = “FASTQ”, minFragLength=35,maxFragLength=151,useAnnotation=”TRUE”, nthreads=64, annot.ext = “my_annotation.gtf.gz”, isGTF = “TRUE”, sortReadsByCoordinates = “TRUE”, output_format = “BAM”) here i have asigned 64 threads but in console, i see only 40 threads, I dont…

Continue Reading Rsubread align maximum nthreads

Output of samtools view, what does the third column actually represent?

The samtools view outputs information from SAM and BAM files in SAM format. You can find a description of the SAM format here: samtools.github.io/hts-specs/SAMv1.pdf Section 1.4 deals with the meaning of each of the manditory coloumns. It includes the following table: Col Field Type Regexp/Range Brief description |—|——|——-|—————————-|—————————————-| 1 QNAME…

Continue Reading Output of samtools view, what does the third column actually represent?

Convert a VCF-file in a user specific Format

Convert a VCF-file in a user specific Format 0 Hello everyone, I am curious if it is possible to convert a VCF-File (with multiple samples) in a Format whith 5 columns. Column should be Sample ID Column: Position on the chromosome Genotyp Number of reads covering site QUAL phred-scaled quality…

Continue Reading Convert a VCF-file in a user specific Format