Paire-End Alignments And Scripture
0
0
Entering edit mode
12.3 years ago
GPR ▴ 390

Hello, I am slowly getting my hands on Scripture and have stumbled on yet another unclear (at least to me )detail. I have paired-end data and the two pairs are in a sigle TopHat-derived *.bam file. I have walked my way through the Tuxedo protocol and didn't need two files per paired-end alignment in Cufflinks or Cuffdiff. My question regarding Scripture is: How do I do I input my alignment file in the "makePairedFile" step?. So far I sorted by queryname with Picard and removed the headers with "sed '1,2d'. Given that I have a single file already, should I jump directly to sort/index with igvtools and then run Scripture with the flag -pairedEnd instead of alignment?

Help will be appreciated. G.

paired-end • 2.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 1684 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6