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%	105.0 / NA	33	D57D (GAC→GAT)	eutT	cobalamin adenosyltransferase involved in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/20);  total (13/20)

CCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGTT  >  NC_000913/2573999‑2574269                                                                                                                                         |                                                                                                                                       ccATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCatcatc                                                                                                                                      >  1:650815/1‑139 (MQ=255)  cATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTtcatcatca                                                                                                                                     >  2:1227043/1‑139 (MQ=255)     tttttCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATcaaca                                                                                                                                  <  1:1469234/139‑1 (MQ=255)           cgccgACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAgg                                                                                                                            <  1:1266014/139‑1 (MQ=255)           cgccgACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAgg                                                                                                                            <  1:568049/139‑1 (MQ=255)                     tGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCt                                                                                                                  <  1:1227043/139‑1 (MQ=255)                     tGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCt                                                                                                                  <  1:525110/139‑1 (MQ=255)                                 tCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAAc                                                                                                       >  2:967304/1‑138 (MQ=255)                                 tCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAAc                                                                                                       <  1:967304/138‑1 (MQ=255)                                                          gACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAg                                                                             <  2:650815/139‑1 (MQ=255)                                                               aGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGc                                                                            >  2:794118/1‑135 (MQ=255)                                                               aGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGc                                                                            <  1:794118/135‑1 (MQ=255)                                                                 ttCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTcc                                                                      >  2:1094033/1‑139 (MQ=255)                                                                      agcaCGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATgcg                                                                                               >  1:562334/1‑109 (MQ=255)                                                                      agcaCGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATgcg                                                                                               <  2:562334/109‑1 (MQ=255)                                                                            ccTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAggg                                                           <  1:434410/139‑1 (MQ=255)                                                                                  gATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTc                                                                                                                 <  1:448942/79‑1 (MQ=255)                                                                                  gATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTc                                                                                                                 >  2:448942/1‑79 (MQ=255)                                                                                     gTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCC‑GGCAGAGGGCGTCAGGCGa                                                 >  2:438665/1‑139 (MQ=255)                                                                                         aTCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTg                                              >  2:995029/1‑139 (MQ=255)                                                                                          tCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGAc                                                <  2:1577340/136‑1 (MQ=255)                                                                                          tCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGAc                                                >  1:1577340/1‑136 (MQ=255)                                                                                          tCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGt                                             <  1:1094033/139‑1 (MQ=255)                                                                                           cGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAgg                                                                                                                            <  2:1538751/59‑1 (MQ=255)                                                                                           cGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAgg                                                                                                                            >  1:1538751/1‑59 (MQ=255)                                                                                             cTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCaa                                                                                                            <  1:1085144/73‑1 (MQ=255)                                                                                             cTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCaa                                                                                                            >  2:1085144/1‑73 (MQ=255)                                                                                             cTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAGACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCg                                          <  2:283438/139‑1 (MQ=255)                                                                                                             gAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTcc                                                                      >  2:64363/1‑95 (MQ=255)                                                                                                             gAACGGGCTGCGGCTGCTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTcc                                                                      <  1:64363/95‑1 (MQ=255)                                                                                                                  ggctgcggctgcTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAAATTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTgcgc                     >  1:1067485/1‑139 (MQ=255)                                                                                                                        ggctgcTGCTGTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGc               <  2:692571/139‑1 (MQ=255)                                                                                                                             ctgctgTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAATTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAg                                       >  1:929920/1‑110 (MQ=255)                                                                                                                             ctgctgTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAg                                       <  2:929920/110‑1 (MQ=255)                                                                                                                               gctgTTCATCATCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGt        <  1:995029/139‑1 (MQ=255)                                                                                                                                    tcatcatcaACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAg                                       <  2:1177724/103‑1 (MQ=255)                                                                                                                                    tcatcatcaACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAg                                       >  1:1177724/1‑103 (MQ=255)                                                                                                                                     catcatCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGtt  <  1:1374443/139‑1 (MQ=255)                                                                                                                                     catcatCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGtt  <  2:1067485/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTCAGGCGACTGTCCGCAGGGAGATGGATCTCTGCGCCTTCGCTGAGCGTATGGTT  >  NC_000913/2573999‑2574269

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