Bwa Sampe Segmentation Fault
Hi everyone, I'm running bwa in the sampe mode and, after successfully processing >10M reads it fails with a segmentation fault (as follows) on what appears like a set of poorly-alignable reads. Any...
View ArticlePaired End Vs Single End Detection
Hi, I am wondering if anyone has a subroutine or function (hopefully in Perl or pseudocode) to detect whether a fastq file is a shuffled paired-end or not. I think I've seen a few different syntaxes on...
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 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 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 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 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 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 ArticleX And Y Chromsome Crossover Position Alignments
How do I find the genotypes on the X chromosome which match the Y SNPs listed in raw data from 23andMe, Ancestry or FTDNA
View Articleinvalid deflate data (invalid code lengths set)
I am trying to trim paired end reads using Trim-Galore. I have made sure that the files match based on the total reads processed in the output txt file from trim-galore. One of the files trimmed...
View ArticleFastQC Quality per tile and per sequence behaving strange after using Cutadapt
I'm using Cutadapt in my paired-end FASTQ files for removing primers (forward primer and rev. com. of the reverse primer in forward reads, and reverse primer and rev. com. of the forward primer in the...
View Article