Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,558,660	A→G	L7S (TTG→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,558,660	0	A	G	100.0%	174.7 / NA	48	L7S (TTG→TCG)	eutB	ethanolamine ammonia‑lyase, large subunit, heavy chain
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/27);  total (21/27)

TCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCG  >  NC_000913/2558517‑2558802                                                                                                                                                |                                                                                                                                               tattttCCGACAACCCCTGCTTTGCCCCCCCGCGCTCCTGTGAGCTTGCCGCGGCACCGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTgg                                                                                                                                          <  2:591607/147‑1 (MQ=255) tCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGt                                                                                                                                         >  1:452359/1‑151 (MQ=255)    tttCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCtt                                                                                                                                      >  2:287324/1‑151 (MQ=255)    ttcccGACAACACCTGTTTTCCCCCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCt                                                                                                                                       <  1:686235/147‑1 (MQ=255)        cGACACCCCCTGCTTTGCCCCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGt                                                                                                                                  <  2:483610/151‑1 (MQ=255)          acaacaCCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGtt                                                                                                                                 >  2:613803/1‑150 (MQ=255)            acccccTGTTTCGCCGCCGCGCCCTCCTGTGGGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGtt                                                                                                                                 <  1:613763/144‑1 (MQ=255)                       tGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGtt                                                                                                                     >  1:344618/1‑149 (MQ=255)                       tGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTc                                                                                                                    <  2:99143/150‑1 (MQ=255)                         ccgccCCGCGCTCCTGTGAGCTGGCCGCTGCGCCCCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCt                                                                                                                 <  2:428615/151‑1 (MQ=255)                          cgccACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTc                                                                                                                >  2:414253/1‑151 (MQ=255)                               cgcgCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGgaa                                                                                                           <  2:102505/151‑1 (MQ=255)                                   cTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAAt                                                                                                       >  1:446114/1‑151 (MQ=255)                                    tCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATc                                                                                                      >  1:799425/1‑151 (MQ=255)                                     ccTGTGTGCTGGCCGCAGCGACCCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATc                                                                                                      <  2:462863/150‑1 (MQ=255)                                        gtgCGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGgaagaa                                                                                                        <  2:751871/145‑1 (MQ=255)                                        gtgAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCgat                                                                                                   <  1:542545/150‑1 (MQ=255)                                         tgAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCgatg                                                                                                  >  2:496503/1‑150 (MQ=255)                                          gcgttggCCGCGGCGCCGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAGCTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCgatga                                                                                                 <  1:39225/144‑1 (MQ=255)                                           aGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGAtt                                                                                               >  2:856220/1‑151 (MQ=255)                                            gCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGAtt                                                                                               <  2:199369/150‑1 (MQ=255)                                            gCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTg                                                                                              >  2:285687/1‑151 (MQ=255)                                                              cgccAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTc                                                                            >  1:525593/1‑151 (MQ=255)                                                                 cAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTaa                                                                         <  2:656449/151‑1 (MQ=255)                                                                  aGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAAc                                                                        <  1:17069/151‑1 (MQ=255)                                                                       aTCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCa                                                                   <  1:104458/151‑1 (MQ=255)                                                                           cccGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGtt                                                               <  2:143450/151‑1 (MQ=255)                                                                            ccGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGttt                                                              <  2:630671/151‑1 (MQ=255)                                                                                   cAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCa                                                                                           >  1:439753/1‑115 (MQ=255)                                                                                   cAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCa                                                                                           <  2:439787/115‑1 (MQ=255)                                                                                       tCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGcc                                                   >  1:542709/1‑151 (MQ=255)                                                                                         gTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCt                                                 <  1:53344/151‑1 (MQ=255)                                                                                            ggCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGt                                              <  1:464927/151‑1 (MQ=255)                                                                                             gCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTc                                             <  1:757026/151‑1 (MQ=255)                                                                                                ttAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCg                                            <  2:346364/149‑1 (MQ=255)                                                                                                      aGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAggt                                    >  2:140533/1‑151 (MQ=255)                                                                                                       gCACCTCTTTTACATCCTTAACCTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAgg                                     <  2:308005/149‑1 (MQ=255)                                                                                                       gCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAggtg                                   <  1:522613/151‑1 (MQ=255)                                                                                                           ctctTTTACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTg                               >  2:236848/1‑151 (MQ=255)                                                                                                                ttACATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTa                           >  1:75858/1‑150 (MQ=255)                                                                                                                  aCATCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAac                        <  1:392821/151‑1 (MQ=255)                                                                                                                    aTCCTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAacac                      >  1:489419/1‑151 (MQ=255)                                                                                                                       cTTAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGcc                   >  1:812601/1‑151 (MQ=255)                                                                                                                        ttAAACTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGccc                  >  1:550320/1‑151 (MQ=255)                                                                                                                           aaCTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGtt               >  2:328304/1‑151 (MQ=255)                                                                                                                            aCTGATATACATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTa              <  2:62762/151‑1 (MQ=255)                                                                                                                                 tataCATTGCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATtcat         >  1:163279/1‑151 (MQ=255)                                                                                                                                         gCCGAACGATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCg  <  2:371941/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               TCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCG  >  NC_000913/2558517‑2558802

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