Hi all!
I am currently working on a project where I am trying to align assembled scaffolds from the Human Microbiome Project to a gene sequence (agaA from C. striatum). I am having trouble with this since every time I try to make an alignment, I get a result like this:
27528 reads; of these:
27528 (100.00%) were unpaired; of these:
27528 (100.00%) aligned 0 times
0 (0.00%) aligned exactly 1 time
0 (0.00%) aligned >1 times
0.00% overall alignment rate
I'm not sure if this is an accurate result. Is there a better way to approach this problem? The reason I'm doing this is because I'm trying to find mutational hot spots in the agaA gene (which I know for sure is found on the skin microbiome) and so I thought that I could perhaps use the data available already in the HMP.
Thanks!
I am currently working on a project where I am trying to align assembled scaffolds from the Human Microbiome Project to a gene sequence (agaA from C. striatum). I am having trouble with this since every time I try to make an alignment, I get a result like this:
27528 reads; of these:
27528 (100.00%) were unpaired; of these:
27528 (100.00%) aligned 0 times
0 (0.00%) aligned exactly 1 time
0 (0.00%) aligned >1 times
0.00% overall alignment rate
I'm not sure if this is an accurate result. Is there a better way to approach this problem? The reason I'm doing this is because I'm trying to find mutational hot spots in the agaA gene (which I know for sure is found on the skin microbiome) and so I thought that I could perhaps use the data available already in the HMP.
Thanks!
Comment