Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,403,803	C→T	T168T (ACC→ACT)	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,803	0	C	T	100.0%	67.6 / NA	26	T168T (ACC→ACT)	hflC	regulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/16);  total (10/16)

CGCATCAGCGATTTCAGCCGTTACTACCTGGCAACGGGTGGTGGCGACATTTCGCAAGCGGAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCG  >  NC_000913/4403585‑4404048                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      cgcATCAGCGATTTCAGCCGTTACTACCTGGCAACGGGTGGTGGCGACATTTCGCAAGCGGAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGGGAAAGATATCGTCACCGATTCCCGTGGTCGTTTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGGACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGc                                                                                                                                                                                                                       <  1:268969/251‑1 (MQ=255)               cAGCCGTTACTACCTGGCAACGGGTGGTGGCGACATTTCGCAAGCGGAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACgg                                                                                                                                                                                                         >  1:68046/1‑251 (MQ=255)                                             cGACATTTCGCAAGCGGAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTgag                                                                                                                                                                                                    >  1:103256/1‑226 (MQ=255)                                             cGACATTTCGCAAGCGGAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGAAACGGCTGAGAAGAAGGGCAAAGTTCCGGtcatc                                                                                                                                                                            >  1:293197/1‑250 (MQ=255)                                                cATTTCGCAAGCGGAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAAcc                                                                                                                                                                        <  1:15826/251‑1 (MQ=255)                                                            ggAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACTCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCTAACAGTATng                                                                                                                                                            <  1:253418/251‑1 (MQ=255)                                                                  gCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGATCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATggcggc                                                                                                                                                        <  1:172641/249‑1 (MQ=255)                                                                       tGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTggg                                                                                                                                                  <  1:142391/250‑1 (MQ=255)                                                                                    ctctGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGCAAGATATCGTCACCGATTCCCGTGGTCGTCTGCCCCTCGCAGTACGTGACGCGCTGACCTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCCGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGAGCTGGGTATTGGAGTTGTCg                                                                                                                                    >  1:157343/1‑251 (MQ=255)                                                                                          ccGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGCCCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATggc                                                                                                                                                           <  1:246143/222‑1 (MQ=255)                                                                                                              gTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCACCCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTc                                                                                                                                     <  1:8828/224‑1 (MQ=255)                                                                                                                   cTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGACGTTACTACTCCGGCGGCAGCTAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGAATCAAGCAGATCAACCTGccgac                                                                                                     >  1:216900/1‑251 (MQ=255)                                                                                                                   cTGGACGTGAAAGATATCGTCACCGANTCCCGTNNTCNTCTGACCCTCGAAGTACGTNACNNNCNNNANTCCGGNTCNNCGGGTACAGNAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCNCGTAACGGCNNAGACGAAGGGCAAAGTTCCGGNCATNNANNCNAACNga                                                                                                                                                                >  2:71996/1‑191 (MQ=255)                                                                                                                           gAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCATAGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTCCAGAAGATGAAGTTACTACTCCGGCGGCCGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGCAgg                                                                                                                                                                                             >  1:172297/1‑155 (MQ=255)                                                                                                                                    cGTCACCGATTCCCGGGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCg                                                                                     >  1:107353/1‑250 (MQ=255)                                                                                                                                         ccGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGAGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACg                                                                                                                                                                                                          <  1:271413/128‑1 (MQ=255)                                                                                                                                               cccGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCGGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGTGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTc                                                                                            <  1:53353/232‑1 (MQ=255)                                                                                                                                                     gtcgtcTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTTCTACTCCGGCGGCAGATAACGCCATTGCCGAATCGGCAGAGCGCTTAACGGCTTAGACGGAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGTGTATTGAAGTTGTc                                                                                                                                     <  1:283789/185‑1 (MQ=255)                                                                                                                                                     gtcgtcTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGTCAGATCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGGATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAgn                                                                                      <  1:61195/232‑2 (MQ=255)                                                                                                                                                     gtcgtcTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATgcgn                                                                   <  1:192635/251‑2 (MQ=255)                                                                                                                                                         tcTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCCGGGGGTTGAAGTTGTCGAtggtcggttcaagc                                                                                                                     >  1:143004/1‑185 (MQ=255)                                                                                                                                                                      gTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCgn                                                                                                                                                                                                                             <  1:199939/80‑2 (MQ=255)                                                                                                                                                                      gTACATGACGCGCTGAACTCCGGTTCTGCGGGTCCAGAAGATGAAGTTACTACTCCGGCGGCCGGTAACGCCAATGCCgc                                                                                                                                                                                                                             >  2:199939/1‑79 (MQ=255)                                                                                                                                                                                   tGAACTCCGGTTCTGCGGGTACAAAAGATGAAGTTACTGCTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGCCTGAGACGAAggg                                                                                                                                                                                            <  1:98266/100‑1 (MQ=255)                                                                                                                                                                                                            agaTGAAGTTACTACTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAAGCCGGACAGTAGGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAGCCGTATGCGCGCCGAGCGTGAAGCGGTAgcg                                              <  1:104032/217‑1 (MQ=255)                                                                                                                                                                                                                     tactacTCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTg    <  1:71996/250‑1 (MQ=255)                                                                                                                                                                                                                      actactCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCTGTCATCAACCCGAACAGTATGGCGGCGCGGGGGAGTGAAGGTGTCGATGTGCGTAGCAAGCGGATCAACCTGCCGACCGAAGTGTCTGAAGCGATGTACAACCGTATGCGCGCCGATCGGGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTgcg  <  1:318309/251‑1 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      CGCATCAGCGATTTCAGCCGTTACTACCTGGCAACGGGTGGTGGCGACATTTCGCAAGCGGAAGTGCTGTTGAAACGTAAGTTCTCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCG  >  NC_000913/4403585‑4404048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑