Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,558,288	T→C	N131S (AAC→AGC)	eutB ←	ethanolamine ammonia‑lyase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,558,288	0	T	C	100.0%	76.3 / NA	23	N131S (AAC→AGC)	eutB	ethanolamine ammonia‑lyase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/11);  total (12/11)

ACTTTGCACGTCGTCACGGGTGTCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGTTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTT  >  NC_000913/2558164‑2558413                                                                                                                             |                                                                                                                              aCTTTGCACGTCGTCACGGGTGTCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCaa                                                                                                                       <  2:403092/133‑1 (MQ=255)            cgtcACGGGTGTCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTAc                                                                                                            >  2:386361/1‑133 (MQ=255)            cgtcACGGGTGTCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTAc                                                                                                            >  2:112644/1‑133 (MQ=255)                     tgtCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGAc                                                                                                   <  1:112644/133‑1 (MQ=255)                      gtCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGAcc                                                                                                  >  2:587318/1‑133 (MQ=255)                      gtCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCCCGCTGGAACAAATCTTCGCTACCGCCGCGAcc                                                                                                  >  2:248594/1‑133 (MQ=255)                              gCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGa                                                                                          >  2:368093/1‑133 (MQ=255)                                               aaGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCggg                                                                         >  1:512078/1‑133 (MQ=255)                                                     cccGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAgg                                                                   >  1:321738/1‑133 (MQ=255)                                                     cccGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACATCCAAGGTCAGCTCTTTAGGGGTAAAgg                                                                   >  2:290584/1‑133 (MQ=255)                                                                   ggtggtATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACg                                                     <  1:368093/133‑1 (MQ=255)                                                                    gtggtATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGc                                                    >  2:494399/1‑133 (MQ=255)                                                                      ggtATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTg                                                  >  1:601870/1‑133 (MQ=255)                                                                            gtccTTTTTATTTACAGGCATTTGTTACCCGCCGTAGATCAAGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCa                                            <  2:308927/131‑1 (MQ=255)                                                                                 ttttGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCt                                       <  2:354163/133‑1 (MQ=255)                                                                                 ttttGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCt                                       <  2:458/133‑1 (MQ=255)                                                                                 ttttGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCt                                       <  2:91397/133‑1 (MQ=255)                                                                                  tttGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAGAGGCAATGTCGTCCACGCTGGTTTCATCGCTc                                      >  2:255970/1‑133 (MQ=255)                                                                                  tttGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTc                                      >  2:327770/1‑133 (MQ=255)                                                                                      aTTACCGGGATTTTCCTCGCGCCGTAGATCGGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGca                                  <  1:248594/133‑1 (MQ=255)                                                                                      aTTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGca                                  <  1:587318/133‑1 (MQ=255)                                                                                        tACCGGCTTTTTCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGcaca                                <  1:386361/133‑1 (MQ=255)                                                                                        tACCGGCAGTTTCTTTGCGCCGGAGATCAGTTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGcaca                                <  1:290584/133‑1 (MQ=255)                                                                                                ttttCTTCGCGCCGTAGATCAGGTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGCACATACTCACg                        <  1:327770/133‑1 (MQ=255)                                                                                                ttttCTTCACGCCGTAGAACAGGTCCGCGCGGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGCACATACTCACg                        <  1:278920/133‑1 (MQ=255)                                                                                                                      gTCCGCGCTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGttt  >  2:57207/1‑133 (MQ=255)                                                                                                                             |                                                                                                                              ACTTTGCACGTCGTCACGGGTGTCATTTGGCTGCAAACGGGCGCTAAAGGTGCCCGGAATACCGATGGTGGTATTGGCCTTTTTGATTACCGGCATTTTCTTCGCGCCGTAGATCAGGTCCGCGTTGGAGCAAATCTTCGCTACCGCCGCGACCACTTCCGAGGTCAGCCCTTTGCGGGTAAAGGCAATGTCGTCCACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTT  >  NC_000913/2558164‑2558413

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