The chicago from joachimwolff

chicago's People

Contributors

Stargazers

Watchers

chicago's Issues

oligo files

Hi,
I want to make oligo files of mm10.
Which script can I make??

Bug in bam2chicago.sh

Hi, found a bug in bam2chicago.sh, I was getting this error from bedtools:

Error: Type checker found wrong number of fields while tokenizing data line.
Perhaps you have extra TAB at the end of your line? Check with "cat -t"

Traced it back to this bit in your code

echo "Flipping all reads that overlap with the bait on to the right-hand side..."
awk 'BEGIN{ OFS="\t" }                                                                                                                                                                                    
     {                                                                                                                                                                                                    
     minRight=$13<$3?$13:$3;                                                                                                                                                                              
     maxLeft=$12>$2?$12:$2;                                                                                                                                                                               
     if($1==$11 && (minRight-maxLeft)/($3-$2)>=0.6){                                                                                                                                                      
              print $4,$5,$6,$1,$2,$3,$7,$8,$10,$9,$11,$12,$13,$14,$15                                                                                                                                    
     }                                                                                                                                                                                                    
     else {                                                                                                                                                                                               
              print $0                                                                                                                                                                                    
     }                                                                                                                                                                                                    
}' ${samplename}/${bamname}_mappedToBaits.bedpe > ${samplename}/${bamname}_mappedToBaits_baitOnRight.bedpe

For some reason, you are trying to output $15 but ${samplename}/${bamname}_mappedToBaits.bedpe only has 14 fields. Remvoing $15 from the line fixed the issue.

Recommend Projects