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%	188.9 / NA	58	*298* (TAG→TAA)	ybcF	putative carbonate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/29);  total (29/29)

TGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATG  >  NC_000913/551208‑551464                                                                                                                             |                                                                                                                                     tGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGa                                                                                                                        <  1:1727628/139‑1 (MQ=255)     gCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACAtt                                                                                                                    >  1:489149/1‑139 (MQ=255)        gAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTtgt                                                                                                                 <  2:786606/139‑1 (MQ=255)                     gTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACg                                                                                                    >  2:976874/1‑139 (MQ=255)                          tGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                                               <  2:542060/139‑1 (MQ=255)                            gTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAcg                                                                                             >  1:1939170/1‑139 (MQ=255)                                     gAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    <  2:489149/139‑1 (MQ=255)                                      aaGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    >  1:796005/1‑138 (MQ=255)                                      aaGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    <  2:796005/138‑1 (MQ=255)                                      aaGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCa                                                                                   >  2:733358/1‑139 (MQ=255)                                       aGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                  >  1:1761804/1‑139 (MQ=255)                                           gTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCt                                                                              >  2:408836/1‑139 (MQ=255)                                           gTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGACGCGACGCTTGACGCGTCTTATCAGGCCt                                                                              >  2:1109314/1‑139 (MQ=255)                                                  aCCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTGCAACCg                                                                       >  1:1271547/1‑139 (MQ=255)                                                     cgcgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGACGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTg                                                                    <  1:1109314/139‑1 (MQ=255)                                                      gcgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGc                                                                   >  1:17840/1‑139 (MQ=255)                                                       cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                                                  >  1:2034228/1‑139 (MQ=255)                                                       cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                                                  >  1:697261/1‑139 (MQ=255)                                                       cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTAGCGGGGGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGAGGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                                                  >  2:1375393/1‑139 (MQ=255)                                                             ttGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATcc                                                            >  2:512056/1‑139 (MQ=255)                                                               ggggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCgg                                                          >  1:116524/1‑139 (MQ=255)                                                                     ttATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTg                                                    <  1:883053/139‑1 (MQ=255)                                                                     ttATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTg                                                    >  1:444788/1‑139 (MQ=255)                                                                        tCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTg                                                 <  1:1087646/139‑1 (MQ=255)                                                                                      gaCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatg                                   <  2:1761804/139‑1 (MQ=255)                                                                                       aCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTa                                                                                                                                 >  2:1830357/1‑44 (MQ=255)                                                                                       aCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTa                                                                                                                                 <  1:1830357/44‑1 (MQ=255)                                                                                       aCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTgg                                                   >  1:363745/1‑122 (MQ=255)                                                                                       aCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTgg                                                   <  2:363745/122‑1 (MQ=255)                                                                                          cTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgat                               >  1:1606663/1‑139 (MQ=255)                                                                                            ggcgggcgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAg                                                                                                                         <  1:1673493/47‑1 (MQ=255)                                                                                            ggcgggcgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAg                                                                                                                         >  2:1673493/1‑47 (MQ=255)                                                                                            ggcgggcgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgatgc                             >  1:1628912/1‑139 (MQ=255)                                                                                                ggcgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGa                         <  2:116524/139‑1 (MQ=255)                                                                                                  cgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACg                       <  2:1939170/139‑1 (MQ=255)                                                                                                           gACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGTGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGa                         <  1:1246489/128‑1 (MQ=255)                                                                                                           gACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGa                         >  2:1246489/1‑128 (MQ=255)                                                                                                           gACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATc              <  2:697261/139‑1 (MQ=255)                                                                                                           gACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATc              <  1:976874/139‑1 (MQ=255)                                                                                                           gACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATc              <  2:17840/139‑1 (MQ=255)                                                                                                              cTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                                               >  1:1068133/1‑55 (MQ=255)                                                                                                              cTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGa                                                                                               <  2:1068133/55‑1 (MQ=255)                                                                                                               tGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAt                                                                                     >  1:1569686/1‑64 (MQ=255)                                                                                                               tGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAt                                                                                     <  2:1569686/64‑1 (MQ=255)                                                                                                               tGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGc          >  2:1894404/1‑139 (MQ=255)                                                                                                                 tATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCt        <  1:1375393/139‑1 (MQ=255)                                                                                                                    ttCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACa     <  2:1018284/139‑1 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAg            <  2:308611/130‑1 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAg            >  1:308611/1‑130 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   <  1:512056/139‑1 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   <  2:1628912/139‑1 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   >  2:605421/1‑139 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   <  2:60775/139‑1 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   <  1:267068/139‑1 (MQ=255)                                                                                                                      cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   >  1:29656/1‑139 (MQ=255)                                                                                                                       gCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATg  <  2:444788/139‑1 (MQ=255)                                                                                                                        cTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt               <  2:158082/125‑1 (MQ=255)                                                                                                                        cTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt               >  1:158082/1‑125 (MQ=255)                                                                                                                             |                                                                                                                                     TGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATG  >  NC_000913/551208‑551464

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