Kuinka laskea lukuja bam-vuoteiden välein sängynvälineillä

Asensin äskettäin Ubuntu 16.04: n (koska käytin edelleen 12.04: tä). Mutta tuntuu siltä, että sängyntyökalujen komentosarjat eivät enää toimi kunnolla. En voi selvittää, miten uusia sänkytyökaluja käytetään vanhoilla tavoillani. Haluan tehdä, että saan lukumäärän bam-tiedostosta, jaksoittain sängyn tiedostosta. Vanhan Ubuntun version ja bedtoolien kanssa se oli hyvin yksinkertaista:

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

Mutta tämä ei näytä toimivan enää. Minulla oli tapana saada pieniä tiedostoja tuloksina (~ 3 Mt), mutta nyt luodaan valtavia tiedostoja (~ 4 Gt), mukaan lukien luetut nimet?

Käytän nyt bedtools v2.25 -ohjelmaa.

Vastaa

-a ja -b järjestys vaihdettiin jossain vaiheessa. Haluat:

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

Tässä on ohjearvon loppu versiossa 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. 

Ja tämä on vastaava ulostulo versiosta 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. 

Päivitä : Muutos käyttäytymisessä tapahtui versiossa 2.24:

Olemme muuttaneet peittotyökalun käyttäytymistä siten, että se on yhdenmukainen Tarkemmin sanottuna peitto lasketaan nyt A-tiedoston aikaväleille päällekkäisyyden perusteella s B-tiedostolla pikemminkin päinvastoin.

Kommentit

  • Kiitos Devon! En ollut ' tietoinen kytkimestä ' '. -A-ja -b vaihtaminen ratkaisee todellakin ongelmani!

Vastaa

Sähköpostiosoitettasi ei julkaista. Pakolliset kentät on merkitty *