Cum se numără citirile în interval bam per pat cu bedtools

Am instalat recent Ubuntu 16.04 (pentru că încă foloseam 12.04). Dar se pare că scripturile mele pentru instrumente de pat nu mai funcționează corect. Nu pot să-mi dau seama cum să folosesc noile instrumente de pat pentru vechile mele moduri. Ceea ce vreau să fac este să obțin numărul de citiri dintr-un fișier bam, pe interval dintr-un fișier de pat. A fost foarte simplu cu vechea mea versiune de Ubuntu și bedtools:

bedtools coverage -abam file.bam -b All_peaks.bed > file.cov.txt 

Dar acest lucru nu pare să mai funcționeze. Obișnuiam să obțin fișiere mici ca rezultate (~ 3MB), dar acum sunt create fișiere uriașe (~ 4 GB), inclusiv numele citite ??

Acum folosesc bedtools v2.25.

Răspuns

Ordinea -a și -b a trecut la un moment dat. Doriți:

bedtools coverage -a All_peaks.bed -b file.bam > file.cov.txt 

Pentru referință, acesta este sfârșitul ajutorului din versiunea 2.25:

Default Output: After each entry in A, reports: 1) The number of features in B that overlapped the A interval. 2) The number of bases in A that had non-zero coverage. 3) The length of the entry in A. 4) The fraction of bases in A that had non-zero coverage. 

Și acesta este rezultatul echivalent din versiunea 2.19:

Default Output: After each entry in B, reports: 1) The number of features in A that overlapped the B interval. 2) The number of bases in B that had non-zero coverage. 3) The length of the entry in B. 4) The fraction of bases in B that had non-zero coverage. 

Actualizați : Schimbarea comportamentului a avut loc în versiunea 2.24:

Am schimbat comportamentul instrumentului de acoperire astfel încât să fie în concordanță cu în mod specific, acoperirea este acum calculată pentru intervalele din fișierul A pe baza suprapunerii s cu fișierul B, mai degrabă decât invers.

Comentarii

  • Mulțumesc Devon! ' nu știam de comutatorul ' '. Comutarea -a și -b rezolvă într-adevăr problema mea!

Lasă un răspuns

Adresa ta de email nu va fi publicată. Câmpurile obligatorii sunt marcate cu *