SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
GeneMark-ES error: parse_ET.pl line 213 Invalid regression data elximo Bioinformatics 2 10-23-2016 07:40 PM
VarScan- invalid format for pileup RBCT Introductions 1 02-27-2014 09:37 AM
RepeatMasker error: line in unrecognized format ashleyrosek Bioinformatics 1 06-12-2013 07:40 PM
Trimmomatic: invalid FASTQ name line alisrpp Bioinformatics 3 12-11-2012 06:21 PM
srma-0.1.8 error - Cannot check readability of null file. allenday Bioinformatics 1 10-11-2010 11:30 PM

Reply
 
Thread Tools
Old 08-14-2015, 05:54 AM   #1
dena.dinesh
Member
 
Location: Europe

Join Date: Feb 2013
Posts: 58
Default Error: Invalid format for pileup at line null in varscan2

Hi,

I am using Varscan2 with readcounts like follows

Code:
samtools mpileup -f /home/tsekara/Transcriptomes/trinity5_frameDPcorrected.fa /home/tsekara/David_Schmidt/Bowtie2_default/Sorted_bams/Sorted_Sample_60_EdU48h1mg.bam | java -jar ./bin/VarScan.v2.3.9.jar readcounts --output-vcf > sample.vcf
but after sometime it gives me error saying that
Code:
Error: Invalid format for pileup at line null
How could I solve this error. Kindly guide me.
dena.dinesh is offline   Reply With Quote
Old 12-14-2015, 07:27 AM   #2
RobUK
Junior Member
 
Location: UK

Join Date: Dec 2015
Posts: 1
Default mpileup format

Hi,

I had the same issue and it turns out that VarScan readcount crashes if the fourth column (number of reads covering the site) is 0.
Try to filter for those manually, e.g.:


HTML Code:
samtools mpileup -f /home/tsekara/Transcriptomes/trinity5_frameDPcorrected.fa /home/tsekara/David_Schmidt/Bowtie2_default/Sorted_bams/Sorted_Sample_60_EdU48h1mg.bam |awk '($4>0)' | java -jar ./bin/VarScan.v2.3.9.jar readcounts --output-vcf > sample.vcf
That worked for me.

Rob

Last edited by RobUK; 12-14-2015 at 07:56 AM. Reason: copy and paste error in the awk command
RobUK is offline   Reply With Quote
Reply

Tags
pileup, rnaseq, variant analysis, varscan2

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 10:57 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO