Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,555,565	G→T	A55E (GCG→GAG)	eutK ←	putative structural protein, ethanolamine utilization microcompartment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,555,565	0	G	T	100.0%	152.6 / NA	43	A55E (GCG→GAG)	eutK	putative structural protein, ethanolamine utilization microcompartment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/21);  total (22/21)

GTCTGGCCGACCGATCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCGCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGACTGAG  >  NC_000913/2555489‑2555635                                                                             |                                                                       gTCTGGCCGACCGATCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCAc                                                                     >  1:465841/1‑80 (MQ=255) gTCTGGCCGACCGATCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCAc                                                                     <  2:465840/80‑1 (MQ=255) gTCTGGCCGACCGATCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACg                                                                    <  2:465839/81‑1 (MQ=255)      gCCGACCGATCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACAcg                                                               <  2:465842/81‑1 (MQ=255)        cgaccgaTCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACAcgcc                                                             <  2:465843/81‑1 (MQ=255)             gaTCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACAcgccgc                                                           >  2:465900/1‑78 (MQ=255)              aTCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACAcgccgc                                                           >  1:465905/1‑77 (MQ=255)               tctcCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACAcgccgc                                                           >  2:465901/1‑76 (MQ=255)                ctcCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACAcgccgc                                                           >  2:465895/1‑75 (MQ=255)                   cTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAg                                                        <  2:465844/75‑1 (MQ=255)                      gCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCt                                                 <  1:465845/79‑1 (MQ=255)                      gCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCt                                                 <  2:465846/79‑1 (MQ=255)                       cGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTc                                               <  2:465847/80‑1 (MQ=255)                        ggCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTcca                                             >  1:465903/1‑81 (MQ=255)                        ggCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTcca                                             >  1:465922/1‑81 (MQ=255)                           tgatgaCACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTccacca                                          >  2:465907/1‑81 (MQ=255)                            gatgaCACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTccaccac                                         >  1:465916/1‑81 (MQ=255)                               gaCACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAg                                      >  1:465909/1‑81 (MQ=255)                                    gCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGtt                                  >  1:465919/1‑80 (MQ=255)                                    gCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTa                                 <  1:465848/81‑1 (MQ=255)                                       cGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGc                              <  2:465849/81‑1 (MQ=255)                                         gTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGa                            >  2:465915/1‑81 (MQ=255)                                            cGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGAc                           >  1:465926/1‑79 (MQ=255)                                            cGCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCa                         <  1:465850/81‑1 (MQ=255)                                             gCATCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAg                        <  1:465851/81‑1 (MQ=255)                                               aTCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATTGCCTTCCACCACCAGCGTTAAGCGACCAggg                      >  2:465908/1‑81 (MQ=255)                                               aTCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAggg                      >  2:465917/1‑81 (MQ=255)                                                tCGCGGCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGt                     <  2:465852/81‑1 (MQ=255)                                                   cgGCAGCACAACCGGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCaa                  <  1:465853/81‑1 (MQ=255)                                                    ggCAGCACAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAg                 >  2:465902/1‑81 (MQ=255)                                                     gcagcaCAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAg                 >  2:465920/1‑80 (MQ=255)                                                       agcaCAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCAc              >  1:465930/1‑81 (MQ=255)                                                       agcaCAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCAc              >  2:465899/1‑81 (MQ=255)                                                       agcaCAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCAc              <  1:465854/81‑1 (MQ=255)                                                        gcaCAACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACt             <  1:465855/81‑1 (MQ=255)                                                          acaACCTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTc           <  2:465856/81‑1 (MQ=255)                                                              ccTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGa       <  1:465857/81‑1 (MQ=255)                                                               cTGCGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGAc      <  1:465858/81‑1 (MQ=255)                                                                  cGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGACTg    <  2:465859/80‑1 (MQ=255)                                                                  cGTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGACTg    <  1:465860/80‑1 (MQ=255)                                                                   gTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGACTGAg  >  1:465913/1‑81 (MQ=255)                                                                   gTCCAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGACTGAg  >  2:465911/1‑81 (MQ=255)                                                                     ccAGCGCCTCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGACTGAg  >  1:465932/1‑79 (MQ=255)                                                                             |                                                                       GTCTGGCCGACCGATCTCCTTGCGGCTGATGACACGCCCGGTACGCATCGCGGCAGCACAACCTGCGTCCAGCGCCGCACGACACGCCGCCAGATCGCCTTCCACCACCAGCGTTAAGCGACCAGGGTCAAGCACTTCGTGACTGAG  >  NC_000913/2555489‑2555635

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