Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,863,006	T→C	G224G (GGT→GGC)	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,863,006	0	T	C	100.0%	159.4 / NA	42	G224G (GGT→GGC)	ygbK	FliA‑regulated DUF1537 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/23);  total (19/23)

ACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTA  >  NC_000913/2862859‑2863143                                                                                                                                                    |                                                                                                                                          aCAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTtctggcctg                                                                                                                                        <  2:68457/151‑1 (MQ=255)           gATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATtggcc                                                                                                                              <  2:531774/151‑1 (MQ=255)                     gTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCa                                                                                                                    <  1:265015/151‑1 (MQ=255)                       cAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGt                                                                                                                  <  2:282128/151‑1 (MQ=255)                        aaGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTg                                                                                                                 >  2:270975/1‑151 (MQ=255)                             cTGGCCCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGATTGCTGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTATCGGGCGGTTCTGGCATGGCGATTGGCCTGGCCCGGCAGTGGgc                                                                                                              >  1:413215/1‑149 (MQ=255)                                   cGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGaa                                                                                                       <  2:312689/150‑1 (MQ=255)                                      ttACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACg                                                                                                   >  1:157112/1‑151 (MQ=255)                                      ttACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACg                                                                                                   >  1:182542/1‑151 (MQ=255)                                             aaGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAAcc                                                                                             >  1:620267/1‑150 (MQ=255)                                               gagGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAgg                                                                                          <  1:470982/151‑1 (MQ=255)                                                     tACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGc                                                                                     >  2:540794/1‑150 (MQ=255)                                                            aCGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAgg                                                                             >  1:228067/1‑151 (MQ=255)                                                            aCGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAgg                                                                             <  2:34123/151‑1 (MQ=255)                                                                 gTGCTTGATGCGCTGACCGAACACCATATGGAAATTCAGGGAGAAGCCTTGCGCGATGCACCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTc                                                                        >  1:118293/1‑151 (MQ=255)                                                                   gCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTcc                                                                       <  2:434444/150‑1 (MQ=255)                                                                   gCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTcc                                                                       <  2:641073/150‑1 (MQ=255)                                                                      tGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCcgctc                                                                   >  1:435624/1‑151 (MQ=255)                                                                       gATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCcgctcg                                                                  <  2:102297/151‑1 (MQ=255)                                                                            gcTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTggg                                                             <  2:493942/151‑1 (MQ=255)                                                                                 ccGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCgcgg                                                        <  1:584306/151‑1 (MQ=255)                                                                                   gAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCgcggcg                                                      >  1:463775/1‑151 (MQ=255)                                                                                   gAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCgcggcg                                                      >  1:346697/1‑151 (MQ=255)                                                                                    aaCACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGt                                                     >  1:579490/1‑151 (MQ=255)                                                                                      cacCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAg                                                   <  2:144563/151‑1 (MQ=255)                                                                                        ccATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTg                                                 >  1:9983/1‑151 (MQ=255)                                                                                          aTCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGct                                               >  1:526715/1‑151 (MQ=255)                                                                                            cTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGctc                                              >  1:444610/1‑150 (MQ=255)                                                                                             tGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCtctc                                            >  2:79474/1‑151 (MQ=255)                                                                                                 aaTTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCgg                                         >  1:324535/1‑150 (MQ=255)                                                                                                      aGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATg                                   <  2:349889/151‑1 (MQ=255)                                                                                                             aGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCaa                            <  1:421763/151‑1 (MQ=255)                                                                                                                       gATGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAAcc                  <  2:577123/151‑1 (MQ=255)                                                                                                                        aTGCCCCACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCg                 <  2:20775/151‑1 (MQ=255)                                                                                                                           ccccACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCa              <  2:508121/151‑1 (MQ=255)                                                                                                                            cccACTGGTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAg             <  1:409461/151‑1 (MQ=255)                                                                                                                                   gTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGcac      >  1:50045/1‑151 (MQ=255)                                                                                                                                   gTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGcac      <  2:212181/151‑1 (MQ=255)                                                                                                                                   gTAACGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGcaa      <  2:581910/151‑2 (MQ=255)                                                                                                                                      aCGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACAtt   <  1:224453/151‑1 (MQ=255)                                                                                                                                       cGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACAtt   >  2:647786/1‑150 (MQ=255)                                                                                                                                       cGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTa  <  1:111747/151‑1 (MQ=255)                                                                                                                                       cGGGCGGTTCTGGCCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTa  <  2:69880/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                          ACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGCTGACCGAACACCATCTGGAAATTCAGGGAGAAGCCTTGCGCGATGCCCCACTGGTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTA  >  NC_000913/2862859‑2863143

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