Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	551,361	2 bp→AA	intergenic (+29/+165)	ybcF → / ← purK	putative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	551,361	0	G	A	100.0%	171.8 / NA	54	intergenic (+29/+166)	ybcF/purK	putative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/26);  total (28/26)
*	NC_000913	551,362	0	C	A	100.0%	174.1 / NA	54	intergenic (+30/+165)	ybcF/purK	putative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (28/26);  total (28/26)

TGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGG  >  NC_000913/551233‑551481                                                                                                                                 ||                                                                                                                        tGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGa                                                                                                                <  2:578699/139‑1 (MQ=255) tGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGa                                                                                                                <  2:1248112/139‑1 (MQ=255) tGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGa                                                                                                                <  2:906176/139‑1 (MQ=255) tGGTTATGTCAGAAGCCGTGGGAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGa                                                                                                                <  2:870991/139‑1 (MQ=255)  ggTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGAc                                                                                                               >  1:1437137/1‑139 (MQ=255)         tCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCtt                                                                                                        >  2:150474/1‑139 (MQ=255)           aGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTAt                                                                                                      >  1:332574/1‑139 (MQ=255)                             gcgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc                                                                                                     >  1:1309825/1‑122 (MQ=255)                             gcgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc                                                                                                     <  2:1309825/122‑1 (MQ=255)                              cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                                                                   >  2:1470992/1‑139 (MQ=255)                                tGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGc                                                                                 >  1:1211433/1‑139 (MQ=255)                                                         aaGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCa                                                                                                    <  1:1181853/95‑1 (MQ=255)                                                         aaGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCa                                                                                                    >  2:1181853/1‑95 (MQ=39)                                                               cGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc                                                                                                     >  1:1350832/1‑88 (MQ=39)                                                               cGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc                                                                                                     <  2:1350832/88‑1 (MQ=39)                                                               cGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcg                                                  >  2:1293573/1‑139 (MQ=255)                                                                 cTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCaa                                                                        <  2:1020011/115‑1 (MQ=255)                                                                 cTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCaa                                                                        >  1:1020011/1‑115 (MQ=255)                                                                 cTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgat                                                >  2:880341/1‑139 (MQ=255)                                                                  tGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgatg                                               <  1:1293573/139‑1 (MQ=255)                                                                       ggcgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGa                                          <  2:1211433/139‑1 (MQ=255)                                                                            aGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc                                                                                                     <  2:931190/75‑1 (MQ=38)                                                                            aGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc                                                                                                     >  1:931190/1‑75 (MQ=38)                                                                                      tGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                <  2:916686/134‑1 (MQ=255)                                                                                      tGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                >  1:916686/1‑134 (MQ=255)                                                                                        tATTTCGCTGTAGTCGTAGGGATTACACATTTGTGCCAGATAAGACGGTTGCCGCGTATTATCAGGCCTACAACCGGGGCCGCATCCGGCAATTGGTGCACAATGACTGATGCGATGCTTGACGCATCTTATCAGGCCt                         <  1:1470992/139‑1 (MQ=37)                                                                                             cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAAcc                                                                                         <  1:470327/70‑1 (MQ=25)                                                                                             cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAAcc                                                                                         >  2:470327/1‑70 (MQ=25)                                                                                               cTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                <  2:4517/125‑1 (MQ=255)                                                                                               cTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                >  1:4517/1‑125 (MQ=255)                                                                                               cTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATgg                  >  1:1155676/1‑139 (MQ=255)                                                                                                 gTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGt                <  2:1437137/139‑1 (MQ=255)                                                                                                  taatCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTa               >  2:101534/4‑139 (MQ=255)                                                                                                       gTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGa          >  1:68436/1‑139 (MQ=255)                                                                                                       gTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGa          >  1:613112/1‑139 (MQ=255)                                                                                                        tAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACg                                        <  2:1475949/108‑1 (MQ=255)                                                                                                        tAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACg                                        >  1:1475949/1‑108 (MQ=255)                                                                                                          ggCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCg       <  2:332574/139‑1 (MQ=255)                                                                                                          ggCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCg       <  2:613112/139‑1 (MQ=255)                                                                                                               tAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAgg  >  1:356860/1‑139 (MQ=255)                                                                                                                aGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatg                                                    >  1:810872/1‑88 (MQ=255)                                                                                                                aGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatg                                                    <  2:810872/88‑1 (MQ=255)                                                                                                                     tttGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGc                                                                          <  2:1047970/61‑1 (MQ=255)                                                                                                                     tttGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGc                                                                          >  1:1047970/1‑61 (MQ=21)                                                                                                                     tttGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt                    <  2:81834/115‑1 (MQ=255)                                                                                                                     tttGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt                    >  1:81834/1‑115 (MQ=255)                                                                                                                       tgtgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                >  2:491971/1‑101 (MQ=255)                                                                                                                       tgtgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                >  2:339303/1‑101 (MQ=255)                                                                                                                       tgtgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                <  1:491971/101‑1 (MQ=255)                                                                                                                       tgtgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                <  1:339303/101‑1 (MQ=255)                                                                                                                         tgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt                    <  2:539409/111‑1 (MQ=255)                                                                                                                         tgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt                    >  1:539409/1‑111 (MQ=255)                                                                                                                           ccTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                >  1:915077/1‑97 (MQ=255)                                                                                                                           ccTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                                <  2:915077/97‑1 (MQ=255)                                                                                                                                 ||                                                                                                                        TGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGG  >  NC_000913/551233‑551481

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: ‑