{"action":"created","issue":{"url":"https://api.github.com/repos/chiulab/surpi/issues/15","labels_url":"https://api.github.com/repos/chiulab/surpi/issues/15/labels{/name}","comments_url":"https://api.github.com/repos/chiulab/surpi/issues/15/comments","events_url":"https://api.github.com/repos/chiulab/surpi/issues/15/events","html_url":"https://github.com/chiulab/surpi/issues/15","id":100040546,"number":15,"title":"AMI run returns no results for 200k reads containing viruses","user":{"login":"bede","id":1007616,"avatar_url":"https://avatars.githubusercontent.com/u/1007616?v=3","gravatar_id":"","url":"https://api.github.com/users/bede","html_url":"https://github.com/bede","followers_url":"https://api.github.com/users/bede/followers","following_url":"https://api.github.com/users/bede/following{/other_user}","gists_url":"https://api.github.com/users/bede/gists{/gist_id}","starred_url":"https://api.github.com/users/bede/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/bede/subscriptions","organizations_url":"https://api.github.com/users/bede/orgs","repos_url":"https://api.github.com/users/bede/repos","events_url":"https://api.github.com/users/bede/events{/privacy}","received_events_url":"https://api.github.com/users/bede/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":3,"created_at":"2015-08-10T11:25:46Z","updated_at":"2015-09-23T00:29:38Z","closed_at":null,"body":"I ran SURPI to test it on a small sample comprising 200k already prefiltered Illumina reads (93MB fastq). I know there are viruses in there... Kraken spots perfect 31mer matches for a handful of RefSeq viruses. I was hoping to get some taxonomic assignments, but SURPI doesn't seem to have generated any useful information in its output, with most of the output files being empty.\r\n\r\nI deployed the AMI to a m4.2xlarge instance and rsynced the reads across and then ran SURPI.sh in the default comprehensive mode. The fastq is 60MB after cutadapt, which seems reasonable. I get a SAM file from SNAP, but don't get any taxonomic information, nor a de novo assembly.\r\n\r\nCan anyone suggest what might be going wrong here? I'd love to be able to use and cite SURPI.\r\n\r\nHere is a screenshot of the output file listing:\r\nhttps://www.dropbox.com/s/1a80v3raf03nem9/Screen%20Shot%202015-08-10%20at%2012.03.44.png?dl=0"},"comment":{"url":"https://api.github.com/repos/chiulab/surpi/issues/comments/142461844","html_url":"https://github.com/chiulab/surpi/issues/15#issuecomment-142461844","issue_url":"https://api.github.com/repos/chiulab/surpi/issues/15","id":142461844,"user":{"login":"sfederman","id":4042632,"avatar_url":"https://avatars.githubusercontent.com/u/4042632?v=3","gravatar_id":"","url":"https://api.github.com/users/sfederman","html_url":"https://github.com/sfederman","followers_url":"https://api.github.com/users/sfederman/followers","following_url":"https://api.github.com/users/sfederman/following{/other_user}","gists_url":"https://api.github.com/users/sfederman/gists{/gist_id}","starred_url":"https://api.github.com/users/sfederman/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/sfederman/subscriptions","organizations_url":"https://api.github.com/users/sfederman/orgs","repos_url":"https://api.github.com/users/sfederman/repos","events_url":"https://api.github.com/users/sfederman/events{/privacy}","received_events_url":"https://api.github.com/users/sfederman/received_events","type":"User","site_admin":false},"created_at":"2015-09-23T00:29:38Z","updated_at":"2015-09-23T00:29:38Z","body":"@voutcn \r\n\r\nYou're also likely getting empty coverage plots due to your installed version of matplotlib, which looks to be at least v1.3.1. This issue should be able to be fixed by using the latest version of coveragePlot.py from the repository. (matplotlib v1.3.1 deprecated the load command).\r\n\r\nAs for the empty countable, it's hard to say what is happening. It looks like SURPI ran properly for you at least through the alignment to NT.\r\n\r\nAlso - would you mind starting a new thread with this topic? The original thread is regarding issues on the AMI, which you do not appear to be using here. I'll move my response if you can do so.\r\n\r\nThanks,\r\n\r\n-Scot"}}