Understanding Samtools Flagstat Output
The following is the output of samtools flagstat command on bam file (paired-end) generated after markDuplicate of Picards.7417232 + 0 in total (QC-passed reads + QC-failed reads) 287618 + 0 duplicates...
View ArticleMerging Illumina Paired End Reads
Dear All,I have fastq a dataset containing forward and reverse sequences obtained through paired end module of Illumina platform. I am trying to merge these paired end reads. I have a query which I...
View ArticleMapping Applications On Hadoop Cluster
Hello,I would like to know if there are any new mapping applications for paired end data designed to work on hadoop cluster. The mapping applications that i know are CloudBurst (single ended), Crossbow...
View ArticleHow does the calculated inner-distance between reads change after trimming...
If we had an insert length of 426bp and an adapter size of 65bp, then the length of the central region between two 100bp reads would be 96bp (426bp - 130bp - 100bp - 100bp). So providing the mean...
View ArticleHow To Count Stand-Specific Paired-End Rna-Seq Reads Overlapping Known...
Dear BiostarsDoes any one how to overlap stand-specific paired-end RNA-Seq reads (BAM) with known protein coding genes (BED) ?I tried the following but I think it is not the correct way ? Would...
View ArticleAligning Reads To Specific Chromosome Using Bwa
Hello Everyone,I have whole genome illumina paired end reads and I want to align my reads to specific chromosome (chr 21) using BWA. I was thinking of aligning the entire reads to fasta file of the...
View ArticleIs It Ok To Use One End Of A Set Of Paired-End Reads As A Set Of Single Reads?
Hi,I'm wondering if it's OK to use one end of a set of paired-end reads as a set of single reads?I just want to get a single read data set of a sample, but currently there's only paired-end data...
View ArticleCounting Reads On Paired-End Strand-Specific Rnaseq Data
Hello everybody, I've a strand-specific paired-end library on which I'd like to perform some standard DGE analysis. However I'm quite unclear about how to go about counting reads. Usually when I have a...
View ArticleCombination Of Paired-End And Single-End Samples In Chip-Seq Tf Study
I have 2 batches of chip-seq samples:(A) One biological SE replicateThis batch, actually, consists of 4 SE Chip-seq samples - one treatment, one control, both have IP controls. All are SE sequenced at...
View ArticleUsing Paired End And Orphaned Singles For De Novo Assembly
I have been using FastX to process reads prior to de novo assembly and mapping. What I have discovered and few have pointed out is the FastX will delete reads leaving reads unpaired which changes the...
View ArticleHow To Extract Information From Fastq Pair End Files
dear BioStars users,I would like to extract from my pair-end fastq files information how many times my read is occurring in my fastq file. So output could look - my read (sequence) - how many times I...
View ArticleDoes Every Read In Paired-End Sam File Have The 0X0001 Flag?
Hello everyone, I have a simple, generic, question about SAM format flags in paired-end Illumina data. Does every read in a SAM file from a paired-end sequencing run automatically have the 0x0001 flag?...
View ArticlePaired-End Bam Files
Hi,Having two BAM files from NGS data, how can one check if they are the BAM files (left and right) from a paired end mapping of the same sample? Thanks for the help.
View ArticleClik here to view.
Trinity Question
Hi， I 'd like to use Trinity to analyze my strand-specific pair-end sequencing data (dUTP) and I am very confused about how to choose the argument. I got two reads files, _1.fa and _2.fa. For...
View ArticleOrientation in paired-end sequencing?
I am new to bioinformatics and currently learning how to use Bowtie 2. As written in the manual: A pair that aligns with the expected relative mate orientation and with the expected range of distances...
View ArticleHow To Convert A Soapsnp Output File To Soap, Sam Or Bam Formats
Hello,I would like to know if there is any tool to convert SOAPsnp to SOAP? Once it is converted to SOAP, I can convert it to SAM using Soap2Sam from reseqtools. Hopefully using reseqtools is not that...
View ArticleFilter Paired-End Sam File For Xt:A:U
Dear all,I have a sam file (BWA output, paired-end reads). I would like to retain only reads which are "properly paired". This I would do by:samtools view -f 0x002 file.sam >...
View ArticleResampling Fastq Sequences Without Replacement
Hello, I want to extract a random sample (without replacement) of 7.5 million fastq sequences from illumina sequencing data that contains about 30 million sequences each in of the reads. I want to...
View ArticleIdentifying mutations from Paired-End Sequencing data
Hello! I'm trying to get mutations from paired-end sequenced reads aligned with BWA using SamTools. Coverage is about 16,000. Generally it works fine, but in one fragment (TGGGC) i see that in reads...
View ArticleEfficiently join paired-end read coordinates in the same line?
Dear community, I have a huge paired-end HiC dataset (BAM format) which I want to format like this way: HWI-D00283:117:C5KKJANXX:2:1101:1139:77789 chr6 153338506 153338556 37 + chr6 153338031 153338081...
View Article