Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	551,332	G→A	*298* (TAG→TAA)	ybcF →	putative carbonate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	551,332	0	G	A	100.0%	122.0 / NA	39	*298* (TAG→TAA)	ybcF	putative carbonate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/21);  total (18/21)

TGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGG  >  NC_000913/551200‑551465                                                                                                                                     |                                                                                                                                      tGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAg                                                                                                                                 >  1:385416/1‑139 (MQ=255)  ggTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAgg                                                                                                                                <  1:457380/139‑1 (MQ=255)             gCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACAtt                                                                                                                     >  1:29830/1‑139 (MQ=255)                 aaTGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTtgtg                                                                                                                 >  2:518992/1‑139 (MQ=255)                                aGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGc                                                                                                         <  2:215216/132‑1 (MQ=255)                                aGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGc                                                                                                         >  1:215216/1‑132 (MQ=255)                                   ggTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAc                                                                                               <  2:29830/139‑1 (MQ=255)                                   ggTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAc                                                                                               <  1:60511/139‑1 (MQ=255)                                    gTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACg                                                                                                     >  1:977163/1‑132 (MQ=255)                                    gTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACg                                                                                                     <  2:977163/132‑1 (MQ=255)                                     ttATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAcgc                                                                                             >  2:422481/1‑139 (MQ=255)                                        tGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTc                                                                                          >  1:611301/1‑139 (MQ=255)                                          tCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCtt                                                                                        >  2:294953/1‑139 (MQ=255)                                             gAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                     <  2:385416/139‑1 (MQ=255)                                                gCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                  <  2:855307/139‑1 (MQ=255)                                                  cGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTa                                                                                                                           >  2:622218/1‑96 (MQ=255)                                                  cGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTa                                                                                                                           <  1:622218/96‑1 (MQ=255)                                                  cGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                >  2:171300/1‑139 (MQ=255)                                                    tGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                                                <  1:222479/121‑1 (MQ=255)                                                    tGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                                                >  2:222479/1‑121 (MQ=255)                                                    tGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTa                                                                              >  2:1062980/1‑139 (MQ=255)                                                                gTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCg                                                                  <  1:124093/139‑1 (MQ=255)                                                                        gggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGc                                                          <  1:422481/139‑1 (MQ=255)                                                                        gggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGc                                                          <  2:841132/139‑1 (MQ=255)                                                                        gggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGc                                                          <  1:518992/139‑1 (MQ=255)                                                                             ttATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCa                                                                                                                              >  2:274475/1‑66 (MQ=255)                                                                             ttATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCa                                                                                                                              <  1:274475/66‑1 (MQ=255)                                                                             ttATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTg                                                     <  2:611301/139‑1 (MQ=255)                                                                              tATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTa                                                                                                                                  <  2:1023154/61‑1 (MQ=255)                                                                              tATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTa                                                                                                                                  >  1:1023154/1‑61 (MQ=255)                                                                                       ttGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                                                <  2:230285/86‑1 (MQ=255)                                                                                       ttGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                                                >  1:230285/1‑86 (MQ=255)                                                                                                    ggcgggcgAAGCGGGGGCCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTATAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgatgc                              >  1:675444/1‑139 (MQ=255)                                                                                                        ggcgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGTATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGa                          <  1:171300/139‑1 (MQ=255)                                                                                                               cGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCaa                                                        >  2:449745/1‑102 (MQ=255)                                                                                                               cGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCaa                                                        <  1:449745/102‑1 (MQ=255)                                                                                                                            ttCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGTTGCGTTGCTTGACGCATCTTATCAGGCCTACa      <  2:675444/139‑1 (MQ=255)                                                                                                                              cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt    <  1:1062980/139‑1 (MQ=255)                                                                                                                                cTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATgg  >  1:1099540/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                      TGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGG  >  NC_000913/551200‑551465

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