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
33
34
35
36
|
#!/usr/bin/python3
import shasta
import argparse
parser = argparse.ArgumentParser(description=
'Write a FASTA file containing all the reads present in a local read graph.')
parser.add_argument('--readId', type=int, required=True)
parser.add_argument('--strand', type=int, choices=range(2), required=True)
parser.add_argument('--maxDistance', type=int, required=True)
parser.add_argument('--allowChimericReads', action='store_true')
parser.add_argument('--allowCrossStrandEdges', action='store_true')
parser.add_argument('--allowInconsistentAlignmentEdges', action='store_true')
arguments = parser.parse_args()
readId = arguments.readId
strand = arguments.strand
maxDistance = arguments.maxDistance
allowChimericReads = arguments.allowChimericReads
allowCrossStrandEdges = arguments.allowCrossStrandEdges
allowInconsistentAlignmentEdges = arguments.allowInconsistentAlignmentEdges
a = shasta.Assembler()
a.accessAlignmentData()
a.accessReadGraph()
a.accessMarkers()
a.writeLocalReadGraphReads(
readId=readId, strand=strand,
maxDistance=maxDistance,
allowChimericReads=allowChimericReads,
allowCrossStrandEdges=allowCrossStrandEdges,
allowInconsistentAlignmentEdges=allowInconsistentAlignmentEdges)
|
