Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,891,886	Δ1 bp	coding (698/1365 nt)	sdaA →	L‑serine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,891,884	0	C	.	100.0%	152.4 / NA	35	coding (696/1365 nt)	sdaA	L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (18/17);  total (18/17)

GTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGCAGTTAACGAAGAAAACGCCGC  >  NZ_CP009273/1891757‑1892018                                                                                                                                |                                                                                                                                       gTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGt                                                                                                                             <  1:621821/138‑1 (MQ=255)              gaacgaacTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGt                                                                                                               <  2:75097/139‑1 (MQ=255)               aacgaacTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTc                                                                                                              >  2:714109/1‑139 (MQ=255)                           cTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCg                                                                                                  <  2:709125/139‑1 (MQ=255)                                   cAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTgg                                                                                          >  2:61765/1‑139 (MQ=255)                                   cAAGAAAGAGATCGACGAGTATTTCGCGCATGGCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTgg                                                                                          >  1:743108/1‑139 (MQ=255)                                            gaTCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCg                                                                                 <  1:61765/139‑1 (MQ=255)                                             aTCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGa                                                                                >  1:10712/1‑139 (MQ=255)                                             aTCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGa                                                                                >  2:727089/1‑139 (MQ=255)                                                 acgaGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATaaa                                                                            <  2:743108/139‑1 (MQ=255)                                                 acgaGTATTTCGCGCATGTCTGGCAAACCATGCAAGCATGTATCGATCGCGGGATGAACACCGAAGGGGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATaaa                                                                            <  2:325934/139‑1 (MQ=255)                                                     gTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATaaa                                                                            <  2:106367/135‑1 (MQ=255)                                                     gTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATaaa                                                                            >  1:106367/1‑135 (MQ=255)                                                           cgcgCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGCTGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACg                                                                  >  2:241325/1‑139 (MQ=255)                                                                aTGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATcc                                                              >  2:708968/1‑138 (MQ=255)                                                                       gCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGa                                                         <  1:639244/136‑1 (MQ=255)                                                                       gCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATg                                                      >  2:596833/1‑139 (MQ=255)                                                                            ccATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGGGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCAtt                                                 >  2:39388/1‑139 (MQ=255)                                                                              aTGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGa                                               <  2:2104/139‑1 (MQ=255)                                                                              aTGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGa                                               <  1:99114/139‑1 (MQ=255)                                                                              aTGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGa                                               <  1:714109/139‑1 (MQ=255)                                                                               tGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGAc                                              <  2:186810/139‑1 (MQ=255)                                                                                       tGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACt                                                                          >  2:181396/1‑103 (MQ=255)                                                                                       tGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACt                                                                          <  1:181396/103‑1 (MQ=255)                                                                                       tGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTaaa                                      <  1:777112/139‑1 (MQ=255)                                                                                         tATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTg                                                <  1:583970/127‑1 (MQ=255)                                                                                         tATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTg                                                >  2:583970/1‑127 (MQ=255)                                                                                              atcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGttt                               >  2:209396/1‑139 (MQ=255)                                                                                              atcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGttt                               >  2:672959/1‑139 (MQ=255)                                                                                              atcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGttt                               >  2:636573/1‑139 (MQ=255)                                                                                                cgcgGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCAc                                                                                                                        <  2:17263/48‑1 (MQ=255)                                                                                                cgcgGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCAc                                                                                                                        >  1:17263/1‑48 (MQ=255)                                                                                                                      tACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGCAGTTAACGAAGaaaa        <  1:209396/138‑1 (MQ=255)                                                                                                                      tACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGCAGTTAACGAAGAAAAc       >  2:760854/1‑139 (MQ=255)                                                                                                                           ccAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGCAGTTAACGAAGAAAAcgccgc  >  2:344010/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       GTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGCAGTTAACGAAGAAAACGCCGC  >  NZ_CP009273/1891757‑1892018

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