File: 00README.txt

package info (click to toggle)
python-pysam 0.15.4+ds-3
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, sid
  • size: 27,992 kB
  • sloc: ansic: 140,738; python: 7,881; sh: 265; makefile: 223; perl: 41
file content (32 lines) | stat: -rw-r--r-- 1,095 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
File ex1.fa contains two sequences cut from the human genome
build36. They were exatracted with command:

  samtools faidx human_b36.fa 2:2043966-2045540 20:67967-69550

Sequence names were changed manually for simplicity. File ex1.sam.gz
contains MAQ alignments exatracted with:

  (samtools view NA18507_maq.bam 2:2044001-2045500;
   samtools view NA18507_maq.bam 20:68001-69500)

and processed with `samtools fixmate' to make it self-consistent as a
standalone alignment.

To try samtools, you may run the following commands:

  samtools faidx ex1.fa                 # index the reference FASTA
  samtools view -bt ex1.fa.fai -o ex1.bam ex1.sam.gz   # SAM->BAM
  samtools index ex1.bam                # index BAM
  samtools tview ex1.bam ex1.fa         # view alignment
  samtools pileup -cf ex1.fa ex1.bam    # pileup and consensus
  samtools pileup -cf ex1.fa -t ex1.fa.fai ex1.sam.gz

In order for the script pysam_test.py to work, you will need pysam
in your PYTHONPATH.

In order for the script example.py to work, you will need pysam
in your PYTHONPATH and run

  make all

beforehand.