Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,862,827	C→T	R165C (CGC→TGC)	ygbK →	FliA‑regulated DUF1537 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,862,827	0	C	T	100.0%	140.5 / NA	39	R165C (CGC→TGC)	ygbK	FliA‑regulated DUF1537 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

CCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTG  >  NC_000913/2862682‑2862973                                                                                                                                                  |                                                                                                                                                   ccTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGtgctgc                                                                                                                                               <  2:25577/151‑1 (MQ=255)       ggTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGcgtc                                                                                                                                         >  2:392789/1‑151 (MQ=255)           aaCGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGttt                                                                                                                                     >  1:585627/1‑151 (MQ=255)                       gTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGttt                                                                                                                         <  2:452465/151‑1 (MQ=255)                            tCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGa                                                                                                                     >  1:168851/1‑150 (MQ=255)                             cAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAAc                                                                                                                   >  1:277803/1‑151 (MQ=255)                               ggggTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACa                                                                                                                  >  2:740500/1‑150 (MQ=255)                                  gTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGgt                                                                                                              <  2:78019/151‑1 (MQ=255)                                   tATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGgtg                                                                                                             >  2:530585/1‑151 (MQ=255)                                       tGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTgg                                                                                                          >  1:205667/1‑150 (MQ=255)                                          tCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGc                                                                                                      >  1:457853/1‑151 (MQ=255)                                          tCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGc                                                                                                      >  2:166928/1‑151 (MQ=255)                                            gTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCg                                                                                                    >  2:314345/1‑151 (MQ=255)                                             tAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGt                                                                                                   >  2:261037/1‑151 (MQ=255)                                                 gAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGACGCCGTTCGt                                                                                               >  2:688644/1‑151 (MQ=255)                                                  aaTCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTc                                                                                              <  1:87262/151‑1 (MQ=255)                                                      aaCTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCa                                                                                             >  1:672692/1‑148 (MQ=255)                                                      aaCTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCa                                                                                             >  1:27114/1‑148 (MQ=255)                                                               ccGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCg                                                                                 <  2:310885/151‑1 (MQ=255)                                                                     ccGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTAc                                                                            <  1:392740/150‑1 (MQ=255)                                                                      cGGGATGCGCCATCCCCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACa                                                                           <  1:83606/150‑1 (MQ=255)                                                                       gggATGCGCCATCACCTGGTAAATCCCATGACCGACAGCCATCTTCCACGTCTGGTTGAAGCGCAATCCACAGGCTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACa                                                                           >  2:169676/1‑149 (MQ=255)                                                                        ggATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCa                                                                        <  2:700742/151‑1 (MQ=255)                                                                        ggATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCa                                                                        <  2:468968/151‑1 (MQ=255)                                                                        ggATGCGCCATCACCCGGTAAATCCCAGGACCGACAGCCATCTTCCCCGTCGGGTTGAAGCGTAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTTTGGATGCCGGTCGTCAAAAGCTGGCTCGCTTACAGcg                                                                        <  1:609872/151‑2 (MQ=255)                                                                                 aTCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCt                                                                >  2:653610/1‑150 (MQ=255)                                                                                   cACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTAcc                                                             <  1:740426/151‑1 (MQ=255)                                                                                          tAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTAcg                                                       >  2:469619/1‑150 (MQ=255)                                                                                                    gACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGAt                                            <  2:722924/151‑1 (MQ=255)                                                                                                    gACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGAt                                            >  1:558042/1‑151 (MQ=255)                                                                                                         aCAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCt                                       <  2:681039/151‑1 (MQ=255)                                                                                                             ccaTCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGAcc                                   <  1:761034/149‑1 (MQ=255)                                                                                                                           tGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTgg                      >  2:300023/1‑150 (MQ=255)                                                                                                                            ggTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTgg                      <  1:369667/149‑1 (MQ=255)                                                                                                                            ggTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGaa                    <  1:249825/151‑1 (MQ=255)                                                                                                                            ggTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGa                     <  2:644180/150‑1 (MQ=255)                                                                                                                              ttGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAAt                  >  1:591251/1‑151 (MQ=255)                                                                                                                                     gcAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAggg            >  2:269274/1‑150 (MQ=255)                                                                                                                                              aGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTg  <  2:598630/151‑1 (MQ=255)                                                                                                                                               gggTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACGCCATCTGGAAATTCAGGGAGAAGCCtt   <  1:101502/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTG  >  NC_000913/2862682‑2862973

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