Calculating percentage of reads mapped to a reference genome
After sequencing I received two files, Read1 and Read2 how can I best draw a conclusion about the percentage of paired reads related to my reference?
alignment
• 741 views
•
link
updated
by
shenwei356
•
written
by
trezini •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Quantify the number of sequences in common paired-end sequencing data
written by sil_bioinfo •Hello, I have two fastq files with the forward(Read1) and reverse(Read2) paired reads. How could I count the number of sequences in common between `Read1.fastq` …
-
bowtie2 mapping oritation problem
written by Shaway •Hello everyone, Sorry for post this naive sequencing question. I aligned a paired DNA sample to hg19 by bowtie2 with parameter "-q -N 1 -X …
-
BWA-MEM - Improperly Paired Reads
written by mglasena •After aligning paired-end 100bp reads to a reference genome, I am getting very low properly paired percentage: 369208441 0 total (QC-passed reads + QC-failed reads) …
-
Variant calling and alignment stats
written by nadiabeg.comsats •Hi. I am using samtools flagstats to see the statistics of my alignment file. It looks something like this: I have 10x genomics reads. `575130408 …
-
High G and low A,C,T content in the 1-10bp of Read2 file in a paired-end whole genome bisulfite seq…
written by Floyd •Hi there! I'm new to bisulfite sequencing. Recently received my paired-end WGBS raw file. The Read1 looks OK after fastQC, but the Read2 file has …
-
Hamming distance of paired-end data from target amplicon sequencing
written by jefferson.jss •Hi all, I am analyzing paired-end data of amplicon libraries from a target region of a viral gene. Briefly, I ordered a PCR product where …
-
beginner question BWA program
written by rezahi i want to map my own reads (illumina hiseq) to reference using "bwa mem". i have 2 lane of paired end data (read1.fastq, read2.fastq, …
-
Problem with alignment
written by nayshool •Hello everyone! We have received from BGI china a few samples we sent for exome sequence + Basic analysis. We got the fastQ files and …
-
How To Check That Tophat Works Well
written by M K<p>I used tophat to map my RNAseq reads to human genome hg19 (My data is single-end) and I got some files in the tophat output …
-
Getting Confused With The Flagstat After Pcr Duplicates Removed
written by KJ Lim •<p>Good day. </p> <p>I encountered a situation like below: </p> <p>The flagstat before PCR duplicates removed from paired end mapped reads. </p> <pre><code>:::::::::::::: 0H.flagstat.txt :::::::::::::: …
As you said, mapping reads to the reference genome, with tools like bowtie2 or bwa, and quantifying them with samtools.