import os
from query import consistency
from query import getdata
# from query import queryrna
# queryrna()
orig_pdb = int(os.system("ls data/Pdb/ | wc -l"))
rna_pdb = int(os.system("ls data/OnlyNA/ | wc -l"))
inp_pdb = int(os.system("ls data/Inp/ | wc -l"))
rel_date = int(os.system("grep \"REVDAT\" data/Pdb/*.pdb | awk '{ if ($5 ~ /^0/) print $0}' | wc -l"))
consistency(orig_pdb,rna_pdb,inp_pdb)
# Uncomment to check that the release dates match
# consistency(orig_pdb,rna_pdb,rel_date)
getdata()
================================================================================
"""
import subprocess
import query
# from query import analysis
# from query import consistency
query.queryrna()
orig_pdb = subprocess.check_output('ls data/Pdb/ | wc -l', shell=True)
print "total number of pdb files = %d" % int(orig_pdb)
rna_pdb = subprocess.check_output('ls data/OnlyNA/ | wc -l', shell=True)
print "rna pdbs = %d" % int(rna_pdb)
inp_pdb = subprocess.check_output('ls data/Inp/ | wc -l', shell=True)
print "inp pdbs = %d" % int(inp_pdb)
query.consistency(orig_pdb, rna_pdb, inp_pdb)
# query.makedirs()
# query.download()
query.analysis()
query.helices()
query.getdata()
#####################################################################
# -4-
# Plot number of bases vs. year
# This module automatically generates plots after subsetting the