Question: pindel output error
gravatar for Bioblazer
11 months ago by
Bioblazer40 wrote:

i have been running pindel like so:

./pindel -f h3rv.fasta -i bam_config.txt -c ALL --Ploidy ploidy_input.txt -o Result_pindel.out

-- this is part of bam_config.txt

/biobank/seq/muthuk/Mycobacterium_bovis/R100/bamfiles/sorted_bam/SRR5216939_mapped.sort.bam 241.4 
/biobank/seq/muthuk/Mycobacterium_bovis/R100/bamfiles/sorted_bam/SRR5216940_mapped.sort.bam 242.9 

2 Ids path of bam file along with insert size and sample ID.

Here I used 192 samples but the output shown result only for 10 entries in VCF file

#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  .0      .1      .2      .3      .4      .5      .6      .7      .8      .9
NC_000962.3     84      .       CG      C       .       PASS    END=85;HOMLEN=1;HOMSEQ=G;SVLEN=-1;SVTYPE=DEL    GT:AD   0/0:1252,0      0/0:1025,0      0/0:1306,0      0/0:1345,0      0/0:1427,0      0/0:1464,1      0/0:1466,0      0/0:678,0       0/0:2056,0      0/0:1150,0

Please anybody tell me what went it wrong here and suggest the steps to be followed. Thanks in advance.

next-gen alignment genome • 459 views
ADD COMMENTlink modified 9 months ago by Biostar ♦♦ 20 • written 11 months ago by Bioblazer40

run samtools view -H and check the read groups RG in the bam header. Is there 192 distinct SN names ?

ADD REPLYlink modified 9 months ago • written 9 months ago by Pierre Lindenbaum114k
Please log in to add an answer.


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