Question: pindel output error
gravatar for Bioblazer
11 weeks 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 • 234 views
ADD COMMENTlink modified 10 days ago by Biostar ♦♦ 20 • written 11 weeks 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 10 days ago • written 10 days ago by Pierre Lindenbaum104k
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: 703 users visited in the last hour