Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,557,497	G→A	P395S (CCG→TCG)	eutB ←	ethanolamine ammonia‑lyase, large subunit, heavy chain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,557,497	0	G	A	100.0%	148.6 / NA	46	P395S (CCG→TCG)	eutB	ethanolamine ammonia‑lyase, large subunit, heavy chain
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/24);  total (22/24)

CATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCC  >  NC_000913/2557378‑2557627                                                                                                                        |                                                                                                                                   cATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTg                                                                                                                 <  1:148889/139‑1 (MQ=255)  aTGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTg                                                                                                                 <  1:61396/138‑1 (MQ=255)  aTGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTg                                                                                                                 >  2:61396/1‑138 (MQ=255)  aTGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCA‑CCAGCCGCCTTCCCattgagtaagttgc                                                                                                               >  2:563453/1‑127 (MQ=255)                     cAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAg                                                                                             <  1:604908/139‑1 (MQ=255)                           ccGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGatcat                                                                                       >  2:432341/1‑139 (MQ=255)                                                  aaCTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgt                                                                <  2:697501/139‑1 (MQ=255)                                                  aaCTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgt                                                                <  2:536571/139‑1 (MQ=255)                                                            caGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgtcagcgtcagc                                                      <  2:6928/139‑1 (MQ=255)                                                            caGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgtcagcgtcagc                                                      <  1:432341/139‑1 (MQ=255)                                                              gTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTcagcgtcagcgt                                                    >  1:70223/1‑139 (MQ=255)                                                               tGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGATGTTCTGGTCAGCGTCAGCGtg                                                   >  2:306034/1‑139 (MQ=255)                                                                       tCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGtat                                           >  1:608590/1‑139 (MQ=255)                                                                               tGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcagcag                                   >  2:218989/1‑139 (MQ=255)                                                                                 cGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                 >  2:261967/1‑123 (MQ=255)                                                                                 cGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                 <  1:261967/123‑1 (MQ=255)                                                                                 cGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTc                                 >  1:414210/1‑139 (MQ=255)                                                                                 cGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTc                                 <  1:218989/139‑1 (MQ=255)                                                                                 cGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTc                                 <  1:118967/139‑1 (MQ=255)                                                                                    tggtCTTGTAATTGGGCATGATGGCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcaca                              <  1:563453/139‑1 (MQ=255)                                                                                       tCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTc                                                                    <  2:167741/98‑1 (MQ=255)                                                                                       tCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTc                                                                    >  1:167741/1‑98 (MQ=255)                                                                                       tCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                               >  2:145083/1‑135 (MQ=255)                                                                                       tCTGGTAGTTGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                               <  1:145083/135‑1 (MQ=255)                                                                                               ttGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAgaga                   <  1:306034/139‑1 (MQ=255)                                                                                                tGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgtcagcgtcagc                                                      >  1:226943/1‑103 (MQ=255)                                                                                                tGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgtcagcgtcagc                                                      <  2:226943/103‑1 (MQ=255)                                                                                                tGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgtcagcgtcag                                                       >  1:421265/1‑102 (MQ=255)                                                                                                tGAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGgtcagcgtcag                                                       <  2:421265/102‑1 (MQ=255)                                                                                                 gAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAgaga                   <  2:731243/137‑1 (MQ=255)                                                                                                 gAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAgaga                   >  1:731243/1‑137 (MQ=255)                                                                                                 gAGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATg                 >  1:366796/1‑139 (MQ=255)                                                                                                  aGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATag                      >  2:227931/1‑133 (MQ=255)                                                                                                  aGCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATag                      <  1:227931/133‑1 (MQ=255)                                                                                                   gCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtg                                                   >  1:603030/1‑103 (MQ=255)                                                                                                   gCATGATGTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtg                                                   <  2:603030/103‑1 (MQ=255)                                                                                                     atgatgTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGatcat                                                                                       <  2:489571/65‑1 (MQ=255)                                                                                                     atgatgTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGatcat                                                                                       >  1:489571/1‑65 (MQ=255)                                                                                                         tgTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAg         >  2:279077/1‑139 (MQ=255)                                                                                                         tgTCATCACCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAg         >  1:351496/1‑139 (MQ=255)                                                                                                                aCCCAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGccc  <  2:414210/139‑1 (MQ=255)                                                                                                                   cAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTca                                <  1:838486/106‑1 (MQ=255)                                                                                                                   cAGCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTca                                >  2:838486/1‑106 (MQ=255)                                                                                                                     gCGACATCCCCATGATGTAGTTGCAGCCTGCGGTg                                                                                                     <  2:450621/35‑1 (MQ=255)                                                                                                                     gCGACATCCCCATGATGTAGTTGCAGCCTGCGGTg                                                                                                     >  1:450621/1‑35 (MQ=255)                                                                                                                     gCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                 <  1:726839/87‑1 (MQ=255)                                                                                                                     gCGACATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                 >  2:726839/1‑87 (MQ=255)                                                                                                                        |                                                                                                                                   CATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCC  >  NC_000913/2557378‑2557627

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