Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

bam-readcount step stuck #52

Open
StickHu opened this issue Mar 13, 2023 · 2 comments
Open

bam-readcount step stuck #52

StickHu opened this issue Mar 13, 2023 · 2 comments

Comments

@StickHu
Copy link

StickHu commented Mar 13, 2023

Sorry to disturb you,but I was stuck in the bam-readcount step.

The code provided is that
`and run bam-readcount:

while read line
do
mag=$line

echo $mag

cd SplitBam
cd ${mag}

cp ../../Split/${mag}/${mag}_contigs.tsv ${mag}_contigs.tsv
samtools faidx ../../Split/${mag}/${mag}.fa

echo "${mag}_contigs.tsv"
mkdir ReadcountFilter
for bamfile in *_Filter.bam
do
    stub=${bamfile%_Filter.bam}`

Is ${mag}_contig.tsv the same as ${Cluster_1}/Cluster1_contig.tsv?
I will appreciate it if you can reply as soon as possible.

@StickHu
Copy link
Author

StickHu commented Mar 14, 2023

And if I have several samples, i can't figure out what can do to make the protocol done. Can you help me about that?

@StickHu
Copy link
Author

StickHu commented Mar 14, 2023

Besides, I can't download all the data provided,such as the test data or the simulated reads

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant