Sorting Fastq Files After Trimming (Orphans And Pe)
I have a bunch of Illumina PE data that has been run through fastx trimmer and clipper. I am ready to map these reads, but am needing to create 2 files for paired end reads (the left and right hand...
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 ArticleTrimming Adapters For Paired-End Sequences
Hi all,I got illumina paired end fastq files. They told me to trim read 2 at the beginning for ~20 to 30 bp due to the WGA adapters. Can we find the adapters by looking in to the quality? Which tool is...
View ArticleWarnings In Bowtie Mapping
Hello, I am trying to use bowtie on small synthetic data for short read mapping. My command is ./bowtie -p 8 -t -S hg19 -1 synthetic_sample1.fq -2 synthetic_sample2.fq > bowtie.sam. The alignment...
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 ArticleHg19 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 ArticleJoining Paired-End Illumina Raw Reads
I recently have amplicons sequenced (Illumina PE250) to investigate microbial community. I want to know what quality detection, and what kind of trims of raw reads by what programs should be performed...
View ArticleTake A Subset Of A Fastq Paired-End Sample
Hi,I have two paired-end fastq compressed files coming from HiSeq RNA-SEq experiment, ie., pair.1.fastq.gz and pair.2.fastq.gz.The files are very large so I wanted to just take a few million/thousand...
View ArticleHow To Check If Illumina Fastq Is Single Or Paired End With Minimal Sequence Id
Hi all, I am trying to check if a FASTQ is single or paired end. From wikipedia I saw that default format has to be like this:@HWUSI-EAS100R:6:73:941:1973#0/1but in my case the sequence id is...
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 ArticleBest Cnv Software?
Hello,I know there are many reviews out there, but I can't seem to find exactly what I like.I found this software called mrCaNaVaR, and it's great since it has it's own aligner which is supposed to be...
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 ArticleWhat Better Way To Get The Paired Reads Aligned Against The Reference Genome?
Hi everybody,I have a group of paired reads sequenced using Solid 4 (50bp each mate). I discovered that reads are contaminated by E.coli. My strategy is to align the reads against the reference genome...
View ArticleNeed help using Shrimp2 on paired end color-space SOLiD data.
Hi, I have SOLiD reads which are paried-end (75bp and 35bp) in .csfasta and .QV.qual format. I would like to use Shrimp2 to align them. So far I have been having trouble using it. I used the following...
View ArticleIn Paired-End Data, If One Read Of The Pair Is Unmapped, Is That Pair...
Hello everyone, I have a simple question about generic paired-end Illumina data. If one read of a pair is unmapped, does this automatically mean that the pair is improper and is (un)flagged in a SAM...
View ArticleTool: Trim Adapters Of Paired-End Reads (Fastq)
Trimming adapter sequences of paired-end experiments is sometimes a problem. If you clip the mates in two steps, it migh happen that you loose one mate, but not the corresponding one, resulting in two...
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 ArticleTophat 2 - Both Pairs Map Concordantly
Hi,I have just run the following command tophat2 --solexa1.3-quals -p 12 -r 80 --max-multihits 1 --no-mixed --no-discordant /home/Turkey/Index/turkeyindex /home/Turkey/WTCHG24920061sequence1.fa...
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 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 Article