We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent fd86820 commit 45f51d1Copy full SHA for 45f51d1
1 file changed
SNPgrabber.py
@@ -27,8 +27,8 @@ def TCStoDict(tcs_file, minqual):
27
28
for lines in tcs:
29
name = re.match('^\w*', lines).group(0) # Contig name
30
- quals = re.split(' *', re.search('\|.{16}\|', lines).
31
- group(0)[2:-2].strip())
+ quals = lines.split("|")[3].split()
+
32
SNP = ''
33
for q, b in zip(quals[:-1], ['A', 'C', 'G', 'T']):
34
try:
0 commit comments