Hg19 Strand Information Of Sam Output.
Hi all,I've got Paired-End Illumina data mapped against the Human Hg19. When viewing the SAM output, how can I check if a pair mapped against the forward Hg19 genome sequence or against the reverse...
View ArticlePicard Matequery Slows Process To A Crawl
I'm looking to iterate through an indexed BAM file using picard and perform various tests on both a read and it's mate. For some I would need the full SAMRecord for the mate so I can't just use the...
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 ArticleUnderstanding 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 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 ArticleHow To Assemble Genome Generated With Bac Clones?
I have 2 fastq files from illimina with reads length 250b. Sequences from one file obtained by sequencing from "right" and from "left" in another. This is paired end sequencing. As it is whole genome...
View ArticleMacs Raises Error: No Such File Or Directory
I am trying to run macs14 with sam files from paired-end data + control. Macs14 returns "No such file":sb7904313:line2 $ macs14 -t /Volumes/Data/G6L2_G6L3/s5/clean_paired_sample.sam -c...
View ArticleIs There Any Advantage Of Paired End Sequencing For Chip-Seq?
Hi,I hope that some of you may have a recommendation regarding the use of paired-end sequencing over single-end for ChIP-seq. In principle we expect only a minor improvement of using paired-end...
View ArticleCoverage For Pair-End Rna-Seq, Extend Reads Or Not?
HI, all I have a question when computing region coverage of pair-end RNASeq data. As showed by the sketch map, when computing region coverage, whether I should use actually mapped reads or extended...
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 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 ArticleAnalyzing Older Illumina Paired-End Data
Hi all, I've got some reads from mid-2010 that I'd like to re-align. I'm not sure how best to proceed. These are from a Illumina Genome Analyzer IIx. They're 40-bp, paired-end. Here's what the text...
View ArticleDownload Large Paired-End Rna-Seq And Microarray Data Of The Same Sample (>...
Hi,I was trying to find if there is any large size paired-end RNA-Seq and microarray data of the same sample, but I wasn't able to find as much data as I wanted.Could somebody please point out where I...
View ArticleCrossbow Final Step Failing On Emr
Hello, I am trying to run Crossbow via EMR command line. I managed to complete all the crossbow steps- Alignment with Bowtie, Calling SNPS and Postprocess. I am getting an error in the final step Get...
View ArticleForum: Mapping Of Ngs Short Reads
This is a simple explanation of how the mapping of short reads works !http://www.youtube.com/watch?v=1ZyoI-4ObSA&feature=related see the first 16 min ! It helped me a lot to understand the basic...
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 ArticleHelp Needed To Run Seal For Genome Mapping.
Hello, I managed to build an index file from the reference genome using Seal. Now I am trying to run Seqal but I am running into errors../seqal /user/hadoop/seal/synthetic_prq1 /user/hadoop/seal_output...
View ArticleBWA MEM mate pair rescue
Greetings,Can someone spell out what this option means? I have several guesses, but would rather ask than guess. -P In the paired-end mode, perform SW to rescue missing hits only but do not try...
View Articlejoin non-overlapping paired-end reads
I have a batch of illumina paired end reads with a shorter than anticipated reverse sequence for a ~550bp amplicon.The joining tools which I have used before (i.e. USEARCH, ea-utils) talk about merging...
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 Article