File: WriteRead.py

package info (click to toggle)
shasta 0.14.0-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 29,636 kB
  • sloc: cpp: 82,262; python: 2,348; makefile: 222; sh: 143
file content (19 lines) | stat: -rwxr-xr-x 338 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
#!/usr/bin/python3

import shasta
import sys

helpMessage = """
Invoke with one argument, the id of the read to be written out.
"""

if not len(sys.argv) == 2:
    print(helpMessage)
    exit(1)

readId = int(sys.argv[1])
fileName = str(readId) + '.fasta'

a = shasta.Assembler()
a.getReads().writeRead(readId=readId, fileName=fileName)