Combination 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 ArticleWgsim Mutations In Output After Setting Everything To 0
I was just wondering, is there any useful information on wgsim? Tutorial? Anything? I have been stuck with it for the last 2 weeks. I'm really not sure how to use it. I need it for a project of mine....
View ArticlePaired-End Reads Alignment For Variant Calling ?
I'm trying to do variant calling (SNPs, Indels) from exome-sequencing data, and the sequencing was done with paired end reads. I would like to use BWA for mapping/alignment, followed by PiCard and GATK...
View ArticleEstimating Insert Size From Paired End Data.
Hi,I have paired end data from illumina. To estimate the insert size in silico ( from scratch ), I have aligned the reads as single end reads to the genome ( mouse ). Now I have the two alignment files...
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 ArticlePair End Sequencing Problem
Dear all,I have a quick question on pair end sequencing. I used to work with Illumina without the pair end reads and I have dificulties to understand how the pair end reads work.In the "old" system you...
View ArticleHow To Determine If Library Is Strand-Specific
Hello all, this might be a very basic question but I gave it some thought and don't see a satisfactory answer.Let's say I have a FASTQ file from a sequencing experiment. How can one quickly determine...
View ArticleMapping Trimmed Paired End Reads With Mapsplice
Hi! I'm trying to map illumina paired end using MapSplice. I used script named sickle to trim the bad quality tails and the output looks like this:Paired end 1@TUPAC_0006:1:1:3062:1473#0/1...
View ArticleForum: Attempting to understand pair end sequencing, mainly illumina
Hi I am currently attempted to understand how pair end sequencing works. I understand the basics that it sequences from both ends and if there is an overlapping region they can be joined to create...
View ArticleHow To Read Maq Paired-End Alignemnt Data Using Bioconductor Packages
Hi all I have maq paired-end alignment files that I want to read into R. I have tried to browse several packages and they all seem to depend on ShortRead package of bioconductor which does not...
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 ArticleFastq Sort By Sequence
Dear all,I have some questions about manipulating fastq files. Its the first time I do this, so I want advice, how to do it. Its a special case, because I only want to keep reads when the following...
View ArticleGsnap: How To Output Only Alignments With Minimal Number Of Mismatches?
I am aligning paired-end ChIP-Seq data of drosophila transcription factors.I would like to have a .sam file with the following alignments reported: - A maximum of three mismatches, i.e. three wrong...
View ArticleTool Recommendation Wanted For Cleaning Fasta/Fastq Files To Remove Unpaired...
Hi Everyone, I've been digging around the web trying to find a tool that would allow me to clean-up my paired-end Illumina data before mapping. My pipeline thus far has been to:1) FASTQC - my R1 file...
View ArticleHow To Determine If Paired–End Illumina Rnaseq Reads Are Strand–Specific
I've been provided with more than a billion reads of RNAseq data for a poorly annotated nematode species. They appear to be 2x100 paired-end Illumina reads – I currently know frustratingly little about...
View ArticleAlign Paired End Reads Using Blast
Hi all, Has anyone align illumina paired end reads using BLAST, I used gsnap to do the alignment first, then use BLAST to align the reads which were not mapped by gsnap. It seems that BLAST can only...
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 ArticleCollect Read Pairs Where At Least One Read Is Mapped
I might word my initial question like another post, but I really have the opposite meaning, i think:Filtering multiple flags with SAMtoolsI am trying to remove paired-end reads from a .SAM file where...
View ArticleExisting Tools For Post-Processing After Aligning Paired-End Reads With Blat?
Hi,I'm wondering if there are existing tools that can do post-processing on paired-end reads' BLAT output.More specifically, I'd like the tool to "merge" the alignments from the alignments of each...
View ArticleDna Sequencing Using Abyss-Pe Of Paired End Reads
Hello! I've got to run abyss-pe on 2 files I've got and find the best parameters (k and so on) that would create the best contig coverage (sorry if I'm messing things up, I'm a programmer that had a...
View Article