Question: Empty .cfg after running
gravatar for shannon.soucy
3.0 years ago by
shannon.soucy0 wrote:

I have been trying to run bam2cfg on my bam files generated with bowtie. The first time I tried I didn't have read groups assigned, but I wasn't comparing SV between samples so it didn't really matter- but I got an empty cfg file. So I re-ran bowtie and added the @RG information, thinking this would help- it didn't I'm getting a different error now than I was before- but still an empty .cfg file!

Here is the bam file header with the @RG:

@HD VN:1.0  SO:coordinate
@SQ SN:DL31-cdc21   LN:2955
@SQ SN:tADL-cdc21   LN:3522
@SQ SN:DL1-cdc21    LN:3492
@SQ SN:Hlac-cdc21   LN:2103
@SQ SN:DL31-gyrB    LN:3231
@SQ SN:DL1-gyrB LN:1914
@SQ SN:Hlac-gyrB    LN:1938
@SQ SN:tADL-gyrB    LN:1950
@SQ SN:DL31-polB    LN:5208
@SQ SN:tADL-polB    LN:2730
@SQ SN:DL1-polB LN:2691
@SQ SN:Hlac-polB    LN:2790
@SQ SN:DL31-pol_II  LN:4296
@SQ SN:Hlac-pol_II  LN:4101
@SQ SN:DL1-pol_II   LN:5031
@SQ SN:tADL-pol_II  LN:5067
@SQ SN:DL1-rir1 LN:3144
@SQ SN:DL31-rir1    LN:4302
@SQ SN:Hlac-rir1    LN:3171
@SQ SN:tADL-rir1    LN:3147
@SQ SN:tADL-top6B   LN:2430
@SQ SN:DL1-top6B    LN:2406
@SQ SN:DL31-top6B   LN:4137
@SQ SN:Hlac-top6B   LN:2457
@RG ID:ID:24m8u PL:illumina
@PG ID:bowtie2  PN:bowtie2  VN:2.2.5    CL:"/export/apps/bowtie2-align-s --wrapper basic-0 -x ./merged_mappings/dl_inteins2 -q 24m-3.0u-sick_trim.fastq --rg-id ID:24m8u --rg PL:illumina -S dl_int-3-bowtie2.sam"

and here is the error that i'm getting when I run

Use of uninitialized value $lib in hash element at line 79, <$bam> line 26.
[Thu Mar  3 15:35:38 2016] Closing BAM file
[Thu Mar  3 15:35:38 2016] Send TERM signal for 8930
[Thu Mar  3 15:35:40 2016] samtools pid process 8930 is still there...
[Thu Mar  3 15:35:40 2016] invoking kill -9 on 8930 ...
[Thu Mar  3 15:35:40 2016] Closing samtools process : 8930
Use of uninitialized value $lib in hash element at line 195.

when I run on the test-data I don't have "Use of uninitialized value" warnings.

Can anyone offer any suggestions as to how I need to reformat my bam files so I can get to generate a *.cfg file??

assembly genome • 1.1k views
ADD COMMENTlink modified 3.0 years ago by Pierre Lindenbaum116k • written 3.0 years ago by shannon.soucy0
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: 2236 users visited in the last hour