Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,403,940	A→G	N214S (AAC→AGC)	hflC →	regulator of FtsH protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,403,940	0	A	G	100.0%	83.0 / NA	25	N214S (AAC→AGC)	hflC	regulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/15);  total (10/15)

CGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAG  >  NC_000913/4403805‑4404075                                                                                                                                        |                                                                                                                                        cggcggCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTg                                                                                                                                     >  1:186133/1‑140 (MQ=255)         gATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGccgaccga                                                                                                                             <  2:242521/140‑1 (MQ=255)           tAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAg                                                                                                                           <  2:252873/140‑1 (MQ=255)                   ttGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGaa                                                                                                                   >  1:51416/1‑140 (MQ=255)                    tGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAGCCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAg                                                                                                                  <  2:340996/140‑1 (MQ=255)                    tGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAg                                                                                                                  <  2:324182/140‑1 (MQ=255)                    tGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAg                                                                                                                  >  1:227241/1‑140 (MQ=255)                    tGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAg                                                                                                                  <  1:77126/140‑1 (MQ=255)                    tGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAg                                                                                                                  <  2:227241/140‑1 (MQ=255)                           gCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTa                                                                                                           <  1:130442/140‑1 (MQ=255)                           gCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTa                                                                                                           <  1:318668/140‑1 (MQ=255)                                                gagaCGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCg                                                                                      <  2:51416/140‑1 (MQ=255)                                                               gTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGcgcg                                                                       >  1:152875/1‑140 (MQ=255)                                                                     gTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGcgt                                                                      <  1:296963/135‑1 (MQ=255)                                                                     gTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCa                                                                 <  2:186133/140‑1 (MQ=255)                                                                     gTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGGAGCGcgt                                                                      >  2:296963/1‑135 (MQ=255)                                                                              ccGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcaca                                                        <  2:152875/140‑1 (MQ=255)                                                                                  aCAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcac                                                         <  2:187782/135‑1 (MQ=255)                                                                                  aCAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcac                                                         >  1:187782/1‑135 (MQ=255)                                                                                                cTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGa                                      >  2:169353/1‑140 (MQ=255)                                                                                                    gTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGaaaaa                                  >  1:280624/1‑140 (MQ=255)                                                                                                                tCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACt                      <  2:280624/140‑1 (MQ=255)                                                                                                                  gATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGc                    <  1:169353/140‑1 (MQ=255)                                                                                                                        cGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTa              >  1:210447/1‑140 (MQ=255)                                                                                                                                    aTCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAg  >  2:231178/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        CGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAG  >  NC_000913/4403805‑4404075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑