Question: bam-readcount error Expect Library
0
gravatar for 9521ljh
19 months ago by
9521ljh20
9521ljh20 wrote:

After i do

java -jar varscan.jar processSomatic output.basename.snp

Then, i get output.base.snp.Somatic.hc.filter

Then to execute False Positive Filter:

bam-readcount -q 1 -b 20 -f reference.fasta -l output.base.snp.Somatic.hc.filter tumorfile.bam >varscan.variants.readcounts

but the output give me

Expect Library: LB in BAM

what is wrong? , i can't understand what should i input after " -l "

ADD COMMENTlink modified 19 months ago by Pierre Lindenbaum131k • written 19 months ago by 9521ljh20
0
gravatar for Pierre Lindenbaum
19 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum131k wrote:

as far as I understand, there is no read group with a LB attribute. https://software.broadinstitute.org/gatk/documentation/article?id=6472

ADD COMMENTlink written 19 months ago by Pierre Lindenbaum131k
0
gravatar for ATpoint
19 months ago by
ATpoint40k
Germany
ATpoint40k wrote:

Probably the input file for -l is malformatted. It must be a BED file with chr-start-end. The documentation, depending on which website you check (original VarScan page, GitHub, the VarScan2 how-to paper) is misleading at times. Given a VCF file, just do:

awk '$1 ~ /^#/ {next} {print $1, $2-1, $2}' in.vcf | sort -k1,1 -k2,2n -u > out.bed
ADD COMMENTlink modified 19 months ago • written 19 months ago by ATpoint40k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 875 users visited in the last hour