Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,574,135	G→A	D57D (GAC→GAT)	eutT ←	cobalamin adenosyltransferase involved in ethanolamine utilization

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,574,135	0	G	A	100.0%	83.7 / NA	27	D57D (GAC→GAT)	eutT	cobalamin adenosyltransferase involved in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/12);  total (15/12)

TGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGTT  >  NC_000913/2574019‑2574269                                                                                                                     |                                                                                                                                       tGAGTCAGCGTATCCGGCTTTGTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGTTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCt                                                                                                                  <  1:239670/139‑1 (MQ=255)   aGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCg                                                                                                                >  1:712051/1‑139 (MQ=255)                                         acaGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAAGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCaa                                                                          <  2:798460/139‑1 (MQ=255)                                              tCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTccc                                                                     >  1:658500/1‑139 (MQ=255)                                                       gCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCaa                                                                          >  2:611106/1‑125 (MQ=255)                                                       gCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCaa                                                                          <  1:611106/125‑1 (MQ=255)                                                       gCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAAATCCCGGGCAGAgg                                                            >  2:613237/1‑139 (MQ=255)                                                           gCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGt                                                        <  2:308439/139‑1 (MQ=255)                                                                           aCTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCa                                        >  1:368018/1‑139 (MQ=255)                                                                                 cAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGgaga                                  >  2:606201/1‑139 (MQ=255)                                                                                   gCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATg                                <  2:712051/139‑1 (MQ=255)                                                                                    cccATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAggcgg                                                                                                                         >  1:33010/1‑49 (MQ=255)                                                                                    cccATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAggcgg                                                                                                                         <  2:33010/49‑1 (MQ=255)                                                                                    cccATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATgg                               >  2:695029/1‑139 (MQ=255)                                                                                       aTGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACGAACAGGCGGcc                                                                                                                       >  1:358546/1‑48 (MQ=255)                                                                                       aTGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACGAACAGGCGGcc                                                                                                                       <  2:358546/48‑1 (MQ=255)                                                                                         gAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGAtctc                          >  1:616072/1‑139 (MQ=255)                                                                                             gggCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTgcg                      <  2:348765/139‑1 (MQ=255)                                                                                              ggctgcggctgcTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGAc                                                <  2:499127/112‑1 (MQ=255)                                                                                              ggctgcggctgcTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGAc                                                >  1:499127/1‑112 (MQ=255)                                                                                              ggctgcggctgcTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTgcgc                     >  2:226929/1‑139 (MQ=255)                                                                                              ggctgcggctgcTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTgcgc                     >  1:603657/1‑139 (MQ=255)                                                                                                   cggctgcTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCg                <  1:606201/139‑1 (MQ=255)                                                                                                     gctgctgctGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCt              <  1:695029/139‑1 (MQ=255)                                                                                                      ctgctgctgTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTg             <  1:756167/139‑1 (MQ=255)                                                                                                         ctgctgTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGc          >  1:230774/1‑139 (MQ=255)                                                                                                         ctgctgTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGc          >  1:112138/1‑139 (MQ=255)                                                                                                                tcatcatcaACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGt   >  1:79656/1‑139 (MQ=255)                                                                                                                 catcatCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGtt  <  2:603657/139‑1 (MQ=255)                                                                                                                 catcatCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGtt  <  2:368018/139‑1 (MQ=255)                                                                                                                 catcatCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGtt  <  2:658500/139‑1 (MQ=255)                                                                                                                 catcatCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGtt  <  2:112138/139‑1 (MQ=255)                                                                                                                     |                                                                                                                                       TGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGTT  >  NC_000913/2574019‑2574269

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