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:




24 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

Unknown 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

Unknown 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!

Unknown 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

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

Anonymous 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

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.

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

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

Dawn Drury said...

Good day everybody, This is my testimony on how i won $4,200,201 million I want to use this opportunity to thank Great JOHN for helping me to win the lottery of $4,200,201 million dollars on Mega millions lottery ticket. I have been playing the lottery for the past 6 years now and I have never won. Ever since then I have not been able to win and I was so upset and I needed help to win the lottery. So I decided to go online and search for help, there i saw so many good testimonies about this man called Great JOHN of how he has cast lucky spell lotto for people to win the lottery. I contacted him also and told him I want to win a lottery, he cast a spell for me which I use and I played and won $4,200,201 million dollars. I am so grateful to this man, just in-case you also need him to help you win, you can contact him through his email: Drjohnwinlotterspell@gmail.com, WhatsApp number: +1 (803) 820 2671 ..

Susan Wray said...

The significance of genome cleaning steps shows its positive effect on quality calling exactness and further developed protein-coding potential, which will be helpful to others hoping to work on since a long time ago read gatherings. https://www.theacademicpapers.co.uk/cheap-dissertation-writing-service.php

Anonymous said...

#BuyCovidVaccineCard #covid #pfizer #vaccine #cdc #covidvaccinecertificate #buycovidcertificate #covid-19 #corona #vaccination #certification #cdcvaccinecertificate #covidcertification #corona #coronavaccinecertificate

Do you want to travel abroad an you need a Covid-19 Vaccine Certificate without taking the vaccine?
We are here to help you obtain a valid Covid-19 Vaccine Certificate without taking the vaccine.
Some people do not have time to take the vaccine,
others do not have access to the vaccine.
To travel abroad you need to provide a valid Covid-19 Vaccine Certificate to prove you have taken the vaccine.
Contact Us for your Covid-19 vaccine certificate.

Get you valid card now, buy it here if you do not want to be vaccinated.

immunization record card Honolulu

immunization record card Anchorage

immunization record card Los Angeles ‎

immunization record card Denver

immunization record card Phoenix ‎

immunization record card Chicago

digital vaccine passport

Canadian vaccine passport

vaccine passport for travel

covid-19 vaccination record card nyc

covid 19 vaccination card for sale

get covid vaccine record card with getting vaccinated

covid vaccination record card bc,

covid immunization record card bc,

covid vaccination record card alberta

how to get vaccination records alberta

covid vaccination record card

covid vaccination record card walgreens

covid vaccination record card ontario

covid vaccination record card Toronto

immunization record card california

immunization record card

texas immunization record card

immunization record card texas

missouri immunization record card

kansas immunization record card

washington state immunization record card

immunization record card nc

immunization record card in denver

arizona immunization record card

how to get a new immunization record card

vaccination record card

covid-19 vaccination record card

covid 19 vaccination record card

covid vaccination record card

authentic cdc vaccination record card

vaccination record card for sale

walgreens vaccination record card

travel vaccination record card

vaccination record card cvs

vaccination record card front and back

vaccination record card for travel

vaccine card online

international certificate of vaccination

buy vaccination card

proof of vaccination card

immunization card online

vaccination card for sell

vaccination card sale

Buy COVID vaccine passports in USA

The battle against COVID has been waged for almost two years. With over 2 billion people around the globe now fully vaccinated, some countries have introduced a vaccine passport (certificate) program to allow people with proof of vaccination to travel, return to the office, and participate in public events.

WHATSAPP: +1 (901) 878-9747

Telegram: @Approveddocumente

EMAIL: approvedocs365@gmail.com

BUY CDC COVID-19 VACCINATION RECORD CARD
BUY CANADIAN DIGITAL COVID-19 CERTIFICATE
without taking the vaccine
Buy CDC Vaccination cards online

Erdi Kurniadi said...

Maayong website : Logo
Maayong website : Biografi
Maayong website : lambang
Maayong website : Profil
Maayong website : Terbanyak
Maayong website : One Piece
Maayong website : Terluas
Maayong website : One Piece

Katheryn Veron said...

An all-natural cure that totally eliminates Herpes Virus and Other Stds. This cure follows a clinically proven set of steps that stops, heals and then cures HERPES SIMPLEX VIRUS 2. My African American neighbor used this very same system to get his only daughter cured of Autism which she had suffered for years, and he (the African Father) was kind enough to tell Francis the supplements he used – you know whites hardly withhold information the way we do here. Sad thing! We arranged for a Contact with him Through his email : drrazorherbalhome@gmail.com. I also got to visit his website : https://herbalistrazorherb.wixsite.com/drrazorherbalhome and DOCTOR RAZOR went on to explain to me exactly what supplements to take, and what to do. He was also kind enough to package and ship down the natural supplements I needed to cure the Herpes simplex 2 Needless to say, I followed his directions on how to take the herbal medicine for 18 days and took the supplements as instructed. I am Delighted to announce that my test results came out Negative IN 4 different hospitals. Contact him on WhatsApp Directly +2349065420442

Unknown said...


Dear Sir/Madam

At INDIABULLSL FINANCIAL CONSULT, We offer a competitive rates and fast loan facility, With our group capital fund we can finance your signatory projects and help you to enhance your business plan.

For inquiry and more details, Email me for procedure and terms.
Email: indiabullsfinance366@gmail.com What's App +91-9599047193
Thanks
Mr.Abubakar Abdallah

Scofio morgan said...

Cure herpes forever.. with herbal product…

The best online i have seen…

Thanks to robinsonbuckler@yahoo. com.

urgent loan said...

BUSINESS LOAN PERSONAL LOAN HERE APPLY NOW WhatsApp +918929509036 financialserviceoffer876@gmail.com Dr. James Eric

Physicians Revenue Group Inc said...

Physicians Revenue Group is the best Medical Billing company in Louisiana.

FAST LOAN OFFER said...


QUICK EASY EMERGENCY URGENT LOANS LOAN OFFER EVERYONE APPLY NOW +918929509036 financialserviceoffer876@gmail.com Dr. James Eric