Oct 29, 2018

Some tips to optimize bacterial genome assembly

I just finished a revised genome assembly for a collaborating lab.  We do de novo sequencing of bacterial genomes all the time, sometimes in batches of 50 or 100 different isolates barcoded together on a high-yield HiSeq run. We typically aim for coverage in the 500x to 1000x range with paired-end 150 base reads.

This genome was assembled as part of a "pipeline" script using SPADES, but it did not come out very nicely - over 3000 contigs, most of them small with low coverage.  A FastQC on the raw data looked good (all bases with mean quality  >Q30) except for the last 10 bases show a big drop in quality AND there are Illumina adapters found sometimes as much as 70 bases into the 3' end of quite a lot of the reads.

I also checked some of the largest contigs by a simple BLASTn at the NCBI website, and they matched quite well to one of our favorite bacterial species (e-value = 0.0;  percent identity = 93%). HOWEVER, many of the smaller contigs match the HUMAN genome. This is a big red flag for me, and probably should be checked for every bacterial (or any non-human) de-novo genome assembly.

So here is my improved de novo assembly protocol:

1) To remove the human DNA, I went back to the raw FASTQ files, and aligned to the reference human genome with Bowtie2. The simple way to do this is to use the '--un' flag in the Bowtie2 command, which dumps all unaligned reads to a separate output file.  There are other more sophisticated ways to do this (using SAMtools)  that interrogate if one or both of a set of paired-end reads align concordantly to the genome, but I wanted to be certain to remove any and all human reads. So I aligned each of my paired end data files (R1 and R2) separately and collected the unmatched reads for each.  I found 0.5% human reads in my raw data. This may not seem like a lot, but the assembler will make thousands of contigs from this small amount of contaminant DNA.

This filter will probably remove a few reads (genome regions) that contain simple sequence repeats longer than 25 bases (for example ATATATAT..., or CCCCCCCCC...) that are found in both human and my bacteria, but we know that assembly of repeats is not reliable anyway with 150 base Illumina reads.

2) Using the unmatched files from the human screen, I filtered for quality and removed Illumina adapters with Trimmomatic. This also catches the case where one read matches human and is removed, but the mate-pair is not. These 'unpaired' reads are not included in the 'trimmed' output file. 

3) The trimmed output files are now ready for assembly with SPADES.  I used kmer sizes 21, 33, 55, 77, 99 as recommended by the SPADES manual for 150 bp paired Illumina reads. I also included the --careful flag which is recommended for bacterial genomes with high coverage.

4) I got good results with this assembly - only about 600 contigs, with the largest one containing more than half of the expected genome (N50 > 1 Mb).  This FASTA data file of contigs was small enough to use the NCBI web BLAST server for a single search. I was also able to  Format the output as a text report showing just the top 10 matches for each contig (with no alignments). 

5) In addition to contig length, SPADES also reports the depth of kmer coverage in the header line for each contig in the multi-sequence FASTA file. I observed that all contigs with coverage greater than 100x matched to the same bacterial genus, but contigs with low coverage matched to all sorts of different things - quite a few to E.coli.  I conclude that the low coverage contigs are low abundance contaminants in our sequencing library, and should not be included in the published de novo genome.  [I do not trust pre-filtering by Bowtie against E.coli, we might lose some highly conserved genes]

6) Somebody more diligent than myself (perhaps one of my students in need of a final project in the Intro Bioinformatics course)  can write a Biopython script to filter sequences in a FASTA file by coverage as reported in the SPADES FASTA header.  I did it with some careful work with awk and Excel.  My final genome is 2.2 Mb, with just 39 high coverage contigs which all have highly significant best BLAST matches to the same bacterial genus and over 2000  genes predicted by GeneMarkS.

The takehome message:
1) filter out Human reads from raw sequence data that will be used for de-novo assemblies (I don't know how well this will work for mammals, vertebrates, etc). 
2) filter out low coverage contigs from your final contig assembly FASTA file.


Raw data:




21 comments:

Anonymous said...

POWERFUL LOTTERY SPELL CASTER DR GBOJIE
i am very grateful sharing this great testimonies with you, The best thing that has ever happened in my life is how i win the lottery. I am a woman who believe that one day i will win the lottery.finally my dreams came through when i email Dr gbojie . and tell him i need the lottery numbers. i have come a long way spending money on ticket just to make sure i win. But i never know that winning was so easy until the day i meant the spell caster online which so many people has talked about that he is very great in casting lottery spell, so i decide to give it a try.I contacted this man and he did a spell and he gave me the winning lottery numbers. But believe me when the draws were out i was among winners. i win 1.900.000 million Dollar. Dr. gbojie truly you are the best, with these man you can will millions of money through lottery. i am so very happy to meet these man, i will forever be grateful to you. Email him for your own winning lottery numbers gbojiespiritualspelltemple@gmail.com. OR call him +2349066410185.or check out his website :http://gbojiespiritualtemple.website2.me

Mariela GarcĂ­a said...

I can’t believe this. A great testimony that i must share to all HERPES SIMPLEX VIRUS patient in the world i never believed that their could be any complete cure for Herpes or any cure for herpes,i saw people’s testimony on blog sites of how dr Ero prepare herbal medicine that cure and brought them back to life again. i had to try it too and you can,t believe that in just few weeks i started using it all my pains stop gradually and i had to leave without the pills the doctor gave to me. Right now i can tell you that few months now i have not had any pain,delay in treatment leads to death. Here is his email:(wealthylovespell@gmail.com) whatsapp him with +2348105150446 visit blog http://wealthyspellhome.over-blog.com

Lucy Morgan said...

God is Good!
I promised God that I would share my testimony on this website. I had all the signs of  STD Virus but I was not too sure as to which one. I did a lot of online research and scared myself straight for a whole week before going to see the nurse. She took one look at my genital part and first said that it could just be the anatomy of my body, then she said it looked like genital warts and that I MAY have herpes. I was devastated. She gave me some medicine for the herpes and some cream for the warts. I was also tested for every single STD including herpes. I went home and cried searching the web for all sorts of cures for herpes and awaiting my results. I saw a post whereby Dr. Kham cured Herpes and other diseases, I copied his contacts out and added him on whats app via +2348159922297. The next day my test result was ready and I confirmed  Herpes positive. I told Dr. Kham about my health problems and he assured me of a cure. He prepared his herbal medicine and sent it to me. I took it for 14 days (2 weeks). Before the completion of the 14 days in which I completed the dose, the Blisters and Warts that were on my body were cleared. I went back for check-up and I was told I'm free from the virus. Dr. Kham cures all types of diseases and viruses with the help of his herbal medicine. You can reach Dr. Kham via his email address on (dr.khamcaregiver@gmail.com) or call him on +2348159922297 and whats app him on same number, but You can Also Visit his website to know more about him at > https://drkhamcaregiver.wixsite.com/drkhamcaregiverherba  for more info.  

Anonymous said...

hello dear, are you still suffering from heart break or you are there you are looking for the fruit of the womb or you are there you have been suffering from HIV/AIDS you do not have time to waste crying all day cause that is a total waste of time all you need to do is to kindly contact Dr. Wealthy in his email address: wealthylovespell@gmail.com and explain everything to him, he will bring a solution to your problem.my dear brothers and sisters i want to tell you all that God have been using this man to bless many people in their numerous problems, what is that thing that is eating you up? you do not have to worry Dr. Wealthy is there to solve your problems today. once more contact him through his email address: wealthylovespell@gmail.com or add him up on WhatsApp +2348105150446

Anonymous said...

A GREAT THANKS TO DR.WEALTHY FOR CURING ME FROM GENITAL HERPES I AM SO HAPPY TO GIVE THIS TESTIMONY ARE YOU SUFFERING FROM THIS SICKNESS HERE IS DR WHO CAN HELP
I can’t believe my genital herpes is really cured, oh is by this time last year I start feeling bad about my life, I feel pain everyday of my life am very happy now that am really cured I couldn’t have do this on my own I wish is not God that help me with my helper I was searching the internet about this sickness last 3month when I found about great doctor wealthy, the man that keep his words I write the man email about my problem immediately I get a reply from him asking me to fill a form which I immediately did and send back to him after some mins he reply me that he have work on my cure that I need to provide some materials, which can enable him to work on my cure which I did on the next day of it, after some hours he inform me that he have getting the things needed for the cure and he is about to go on with the curing spell he called me again after 50mins that he is done with the cure that I should check my body and also go for test I cant believe I was negative a big thanks to him am very happy now with my family you can also get your self cured too from this sickness by contacting him through his email: wealthylovespell@gmail.com or whatapp +2348105150446
He also have a herbal cure for 7 other DISEASES;
1.HIV
2.SHINGLES
3.VIRAL HEPATITIS
4.INFLUENZA
5.IMPOTENCE,
6.BARENESS/INFERTILITY,
7.ANTHRAX
8 HPV
Contact him today and you will have a testimony…Good luck!

Anonymous said...

Hello friends,I am ALISHA ROSS, This is my testimony on how my ex lover come back to me..i want to say a very big thanks and appreciation to ZUMA ZUK solution home for bringing back my lover who left me for almost 2 year,i feel like my life is completely over, so one day as i was searching on the internet for recommend spell caster who will help me to bring my lover back, Finally i met a writing how so many testimony talking about how ZUMA ZUK help to restore relationship back,so i decided to give it a try so i contacted the spell caster called ZUMA ZUK and explain my problems to him,and he cast spell love spell luckily within 24 hours my lover really called me and started apologizing for all he had caused me.and be begging me to accept him back and will are living together and happily married, I am the happiest woman on earth today because ZUMA ZUK has done a wonderful deeds in my life and i will continue to share this testimony, i pray God almighty give you the strength and wisdom to help more people having similar problem like mine, contact him. Email spiritualherbalisthealing@gmail.com Call or Whats App +15068001647

Jessica L. Smith said...

After being in relationship with Wilson for seven years,he broke up with me, I did everything possible to bring him back but all was in vain, I wanted him back so much because of the love I have for him, I begged him with everything, I made promises but he refused. I explained my problem to someone online and she suggested that I should contact a spell caster that could help me cast a spell to bring him back but I am the type that don’t believed in spell, I had no choice than to try it, I meant a spell caster called Dr Zuma zuk and I email him, and he told me there was no problem that everything will be okay before three days, that my ex will return to me before three days, he cast the spell and surprisingly in the second day, it was around 4pm. My ex called me, I was so surprised, I answered the call and all he said was that he was so sorry for everything that happened, that he wanted me to return to him, that he loves me so much. I was so happy and went to him, that was how we started living together happily again. Since then, I have made promise that anybody I know that have a relationship problem, I would be of help to such person by referring him or her to the only real and powerful spell caster who helped me with my own problem and who is different from all the fake ones out there. Anybody could need the help of the spell caster, his email: spiritualherbalisthealing@gmail.com or WhatsApp him +15068001647
you can email him if you need his assistance in your relationship or anything. CONTACT HIM NOW FOR SOLUTION TO ALL YOUR PROBLEMS

Donna said...

I am Doctor Paul I got affected with HIV in the process of attending to my HIV patient I tried all I can to get cured but all to no avail, until I saw a post in a health forum about a herbalist man who prepare herbal medication to cure all kind of diseases including HIV virus, at first I doubted if it was real but decided to give it a try, when I contact this herbalist via his email Blessedlovetemple@gmail.com and he prepared a HIV herbal cure and sent it to me via fed-ex delivery company service, when I received this herbal cure, he gave me step by directions on how to apply it, when I applied it as instructed, I was totally cured of this deadly disease within 5 days of usage, I am now free from the deadly disease called HIV, all thanks to Dr Mark. Contact this great herbal spell caster. Kindly contact him. Blessedlovetemple@gmail.com
He cures all kinds of sickness or diseases such as: 1. HERPES VIRUS 2. LASSA FEVER 3. GONORRHEA 4. HIV/AID 5. EX BACK.
Thanks Dr Mark for saving my life.

Jessica L. Smith said...

After being in relationship with Wilson for seven years,he broke up with me, I did everything possible to bring him back but all was in vain, I wanted him back so much because of the love I have for him, I begged him with everything, I made promises but he refused. I explained my problem to someone online and she suggested that I should contact a spell caster that could help me cast a spell to bring him back but I am the type that don’t believed in spell, I had no choice than to try it, I meant a spell caster called Dr Zuma zuk and I email him, and he told me there was no problem that everything will be okay before three days, that my ex will return to me before three days, he cast the spell and surprisingly in the second day, it was around 4pm. My ex called me, I was so surprised, I answered the call and all he said was that he was so sorry for everything that happened, that he wanted me to return to him, that he loves me so much. I was so happy and went to him, that was how we started living together happily again. Since then, I have made promise that anybody I know that have a relationship problem, I would be of help to such person by referring him or her to the only real and powerful spell caster who helped me with my own problem and who is different from all the fake ones out there. Anybody could need the help of the spell caster, his email: spiritualherbalisthealing@gmail.com or WhatsApp him +15068001647
you can email him if you need his assistance in your relationship or anything. CONTACT HIM NOW FOR SOLUTION TO ALL YOUR PROBLEMS https://www.youtube.com/channel/UCShqcSowkKNWKCHEktBm81g

Mariam Max said...

I have been HIV positive for 6 years and long for the day to be free of this disease. I would love to be part of any trial that helped find the cure, i have an undetectable viral load and CD4 count of around 1100.I have tried almost everything but I couldn’t find any solution on my disease, despite all these happening to me, i always spend a lot to buy a HIV drugs from hospital and taking some several medications but no relieve, until one day i was just browsing on the internet when i came across a great post of !Mechelle rose ! who truly said that she was been diagnose with HIV and was healed that very week through the help of this great powerful healing spell doctor ,I wonder why he is called the great Dr, Zuma, i never knew it was all because of the great and perfect work that he has been doing that is causing all this. so I quickly contacted him, and he ask me some few questions and so i did all the things he asked me to do, He ask me to buy some herbs and which I did for my cure, only to see that at the very day which he said i will be healed, all the strength that has left me before rush back and i becomes very strong and healthy, this disease almost take my life all because of me, so i went to hospital to give the final test to the disease and the doctor said i am HIV negative, i am very amazed and happy about the healing doctor Zuma gave to me through ancients herbs and roots, you can email him now for your own healing too on his email: spiritualherbalisthealing@gmail.com or WhatsApp him +15068001647

Doris Ebeneezer said...

Hello viewers! Please share this for others to know about the cure for herpes1&2. I'm here to testify on how i was cured from HERPES 1&2 with the use of herbs. I got the herbs from Dr Odion and he gave me instructions on how to use it. I took the herbal product for a week and 5days after which i go for a blood test and the result was Negative. I waited for another two months and retested the result was still negative and my doctor told me I'm completely free from herpes1&2 and this Gladys my heart and i wish to share it for the benefit of those living with herpes1&2 to also get the herbs from Dr Odion and get cured just like me. You can contact Dr Odion on his Email: drodionherbalhome12@gmail.com or call him on his phone number and WhatsApp +2349019421176 also through his Dr Odion also have the cure on the following diseases HERPES, LUPUS, HPV, CANCER, DIABETES, SHINGLES,VAGINAL INFECTION, HSV1&2 and so many more.

Frank Thomas said...

Good luck to anyone reading this true life story of mine, I Was Diagnosed With type 2 Herpes Virus Last year, And I Was Looking For Solution To Be Cured Luckily I Saw Testimonies On How Dr Odion Cure Herpes Virus I Decided To Contact Dr Odion I Contacted Him He Prepared A Herbal Medicine Portion And Sent It To Me, I Started The Herbal Medicine For My Health. He Gave Me Step By Step Instructions On How To Apply It, When I Applied It As Instructed, I Was Cured Of This Deadly Herpes Within 2 weeks, I Am Now Herpes Negative. My Brother And Sister I No That There Are So Many People That Have The Same Herpes Virus Please contact Dr Odion To Help You Too, And Help Me To Thank Dr Odion For Cure Me, I’m Cured By Dr. Odion Herbal Medicine, His Contact Email: [drodionherbalhome12@gmail.com ] Or Cell Whatsapp Number +2349019421176 thank you

Jessica L. Smith said...

After being in relationship with Wilson for seven years,he broke up with me, I did everything possible to bring him back but all was in vain, I wanted him back so much because of the love I have for him, I begged him with everything, I made promises but he refused. I explained my problem to someone online and she suggested that I should contact a spell caster that could help me cast a spell to bring him back but I am the type that don’t believed in spell, I had no choice than to try it, I meant a spell caster called Dr Zuma zuk and I email him, and he told me there was no problem that everything will be okay before three days, that my ex will return to me before three days, he cast the spell and surprisingly in the second day, it was around 4pm. My ex called me, I was so surprised, I answered the call and all he said was that he was so sorry for everything that happened, that he wanted me to return to him, that he loves me so much. I was so happy and went to him, that was how we started living together happily again. Since then, I have made promise that anybody I know that have a relationship problem, I would be of help to such person by referring him or her to the only real and powerful spell caster who helped me with my own problem and who is different from all the fake ones out there. Anybody could need the help of the spell caster, his email: spiritualherbalisthealing@gmail.com or WhatsApp him +15068001647 blog: spiritualherbal.blogspot.com you can email him if you need his assistance in your relationship or anything. CONTACT HIM NOW FOR SOLUTION TO ALL YOUR PROBLEMS

juwan james said...

THERE IS CURE FOR HERPES VIRUS. This has been a great opportunity for me to speak about Doctor Ishiaku concerning his HERBAL MEDICINE for the cure of HERPES VIRUS. please visit doctor Ishiaku facebook page  
http://www.facebook.com/Dr-Ishiaku-herbal-home-103704578310470    and his website https://www.ishiakuherbalcure.wordpress.com/    I am a living proof of the cure, it is not a story any more, thank God i gave him a try to cure me, i was confirm 100% negative after the usage of the herbal medicine, it is like a dream to me that Doctor Ishiaku herbal medicine could cure  herpes virus to 100% negative, before i proceed further let me quickly share his contact to the general public so anyone suffering from this virus can also have an opportunity like this to be cured.  WHATSAPP/CALL +2348180828544 EMAIL- ishiakuherbalcure@gmail.com .Please let every Herpes patient give herbal medicine a try. It works perfectly well. don't be deceived any more that there is no cure. This herbal medicine was shipped from Africa to me here in the United State by DHL. I am so happy to have seen my health and life back to normal. Thank you Doctor Ishiaku

Elizabeth Glen said...

“So I’m someone’s mom!” Welcoming our first child, Cecily , at 12:24 p.m. on Tuesday, July 14. We are absolutely head over heels in love with Cecily, and parenthood is already the most insane and beautiful thing in existence, It's made me excited to have a little spitfire of a daughter of my own. I remembered when i found out i was pregnant 3years ago and was about to walk away from the musical. But at eight weeks, i had a miscarriage. i was so unhappy, until i seek help spiritually from Priest Babaka the herbal practitioner, who helps and guide me to get pregnant again, even at the trying times few weeks in April when i battled symptoms of the coronavirus including "a cough that makes it feel like my head is splitting open from the inside out, but luckily, the baby was okay with the Priest help. i am happy to finally be a mother, couples out there that needs help, trying to conceive a baby, contact Priest Babaka on babaka.wolf@gmail.com or Facebook at priest.babaka , you will definitely have a baby to make you a parent. He handle Pregnancy related cases.

Mariam Max said...

I have been HIV positive for 6 years and long for the day to be free of this disease. I would love to be part of any trial that helped find the cure, i have an undetectable viral load and CD4 count of around 1100.I have tried almost everything but I couldn’t find any solution on my disease, despite all these happening to me, i always spend a lot to buy a HIV drugs from hospital and taking some several medications but no relieve, until one day i was just browsing on the internet when i came across a great post of !Mechelle rose ! who truly said that she was been diagnosed with HIV and was healed that very week through the help of this great powerful healing spell doctor ,I wonder why he is called the great Dr, Zuma, i never knew it was all because of the great and perfect work that he has been doing that is causing all this. so I quickly contacted him, and he ask me some few questions and so i did all the things he asked me to do, He ask me to buy some herbs and which I did for my cure, only to see that at the very day which he said i will be healed, all the strength that has left me before rush back and i becomes very strong and healthy, this disease almost take my life all because of me, so i went to hospital to give the final test to the disease and the doctor said i am HIV negative, i am very amazed and happy about the healing doctor Zuma gave to me through ancients herbs and roots, you can email him now for your own healing too on his email: spiritualherbalisthealing@gmail.com or WhatsApp him +15068001647
Blog: https://spiritualherbal.blogspot.com

Unknown said...

Hello everyone out there, I'm Patricia and I'm delighted to share this with you as one of the beneficiaries of DR. EKPEN Herbs. After trying all possible means to Get Cured of Hepatitis B Virus, I got to know that this virus is totally curable using Natural herbs. It was like a mere joke or something I never thought could happen because I was in distress and uncomfortable with my body system and symptoms. I want to use this means to encourage those who have given up hope of recovering from whatever Illness you're suffering from, that it's never the end of life you can as well be free and get your healing after using DR. EKPEN herbal products and following his instructions, I put it to the whole world that there is a reliable, recommendable and affordable Hepatitis and other virus and Diseases cure. I will say to everyone reading this article, that I am a living testimony of DR.EKPEN HERBAL PRODUCTS. If you are interested in getting cured of any Virus infections and Diseases,, HIV, HVS 1&2CANCEL etc you can Email DR. EKPEN on (doctorekpen222@gmail.com or whats-app him on +2348102454875
also have cure to all kinds of diseases , Herpes cure , HIV / aid,Cancer , etc

Unknown said...

DR. JOROMI IS TRULY A POWERFUL MAN. Contacting Dr Joromi is the best thing that has ever happened to me. My husband left me for another woman for the reason which i know nothing about. I tried every means to get him back but to no avail, all my efforts proved abortive. I was even contemplating suicide. On a lucky day when i was browsing the internet, i came across some amazing testimonies about Dr. Joromi on how he has helped a lot of people with his powerful love spells. I contacted him and explained all my problems to him and he assured me that he will help me. He told me what to do and assured me that my man will be back in 12 hours. I obeyed him and followed the instructions and to my Greatest surprise, my man called me in 12 hours time just as Dr. Joromi said. Thank you Dr Joromi. May the good God continue to bless you. If you are in need of help please contact Dr Joromi now on his website: https://www.joromispells.com and you can also reach him with his email: Joromispells@gmail.com and here is his WhatsApp number which you can also use to reach him +2348138695387 Thank you Dr Joromi. Thank you once again

Ric Clayton said...

I really want to thank Dr Emu for saving my marriage. My wife really treated me badly and left home for almost 3 month this got me sick and confused. Then I told my friend about how my wife has changed towards me. Then she told me to contact Dr Emu that he will help me bring back my wife and change her back to a good woman. I never believed in all this but I gave it a try. Dr Emu casted a spell of return of love on her, and my wife came back home for forgiveness and today we are happy again. If you are going through any relationship stress or you want back your Ex or Divorce husband you can contact his whats app +2347012841542 or email emutemple@gmail.com website: Https://emutemple.wordpress.com/

ushaservicesonline@gmail.com said...

Do you need a Loan?
Are you looking for Finance?
Are you looking for a Loan to enlarge your business?
I think you have come to the right place.
We offer Loans at low interest rate.
Interested people should please contact us on
For immediate response to your application, Kindly
reply to this emails below only.
Whats app +918256953815
ushaservicesonline@gmail.com
Please, do provide us with the Following information if interested.

1) Full Name:………
2) Gender:………
3) Loan Amount Needed:………
4) Loan Duration:………
5) Country:………
6) Home Address:………
7) Mobile Number:………
8)Monthly Income:…………………
9)Occupation:………………………
)Which site did you here about us…………………

Thanks and Best Regards.
ushaservicesonline@gmail.com
+918256953815

Aaron Jeremy said...

CHIEF DR LUCKY HERBAL THERAPY
Thank you very much for curing my Human Papillomavirus (HPV) totally. You are really a Godsent, I promise you that I will always recommend you. and to share your testimony to people. Contact him now thanks all.
CHIEF DR LUCKY HERBAL THERAPY CAN AS WELL CURE THE FOLLOWING VIRUSES :-
1: HIV/AIDS
2: HSV
3: LUPUS
4: ALS
5: HEPATITIS B
6: DIABETES
7: HUMAN PAPILLOMAVIRUS (HPV)
8: CANCER
9: PAD
10: PENIS ENLARGEMENT
11: THYROID ETC..

WhatsApp him on +2348132777335 or email (chiefdrlucky@gmail.com).

Website: https://chiefdrluckyherbaltherapy.wordpress.com/

Http://facebook.com/chiefdrlucky