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%	182.2 / NA	41	coding (696/1365 nt)	sdaA	L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (23/18);  total (23/18)

AAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCA  >  NZ_CP009273/1891795‑1891969                                                                                          |                                                                                      aaagaGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGGcc                                                                                      >  1:697943/1‑91 (MQ=255)      gaTCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTgc                                                                                >  2:667939/1‑91 (MQ=255)              aGTATTTCGCGCATGTTTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCa                                                                        >  1:508669/1‑91 (MQ=255)                tATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCAcg                                                                      >  2:567535/1‑91 (MQ=255)                   ttCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCAcgtcg                                                                   >  2:68199/1‑91 (MQ=255)                           tGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTg                                                           >  1:109870/1‑91 (MQ=255)                                     aCCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCg                                                 >  1:709382/1‑91 (MQ=255)                                        aTGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGAt                                              >  2:614692/1‑91 (MQ=255)                                        aTGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGAt                                              >  1:190905/1‑91 (MQ=255)                                          gCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCAcgtcgt                                                                  >  2:92506/1‑69 (MQ=255)                                          gCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCAcgtcgt                                                                  <  1:92506/69‑1 (MQ=255)                                          gCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGc                                            <  2:386840/91‑1 (MQ=255)                                               cATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGtt                                       >  1:421917/1‑91 (MQ=255)                                                  gTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTcc                                    >  1:55123/1‑91 (MQ=255)                                                  gTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTcc                                    >  1:151655/1‑91 (MQ=255)                                                  gTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTcc                                    <  1:667939/91‑1 (MQ=255)                                                       gatcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGTGCCGGATGCTGGTTTCCAGCGa                               >  1:392888/1‑91 (MQ=255)                                                       gatcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCg                                                 <  1:363117/73‑1 (MQ=255)                                                       gatcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCg                                                 >  2:363117/1‑73 (MQ=255)                                                       gatcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGAGCCGGATGCTGGTTTCCAGCGa                               >  1:308016/1‑91 (MQ=255)                                                        atcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGAt                              >  2:627982/1‑91 (MQ=255)                                                         tcgCGGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATa                             >  2:819156/1‑91 (MQ=255)                                                           gcgGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATaaa                           >  2:821999/1‑91 (MQ=255)                                                            cgGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAAc                          <  1:720140/91‑1 (MQ=255)                                                            cgGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAAc                          <  1:819156/91‑1 (MQ=255)                                                            cgGGATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAAc                          <  2:339114/91‑1 (MQ=255)                                                              ggATGAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTg                        <  2:536323/91‑1 (MQ=255)                                                                  gAACACCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTa                    >  2:643078/1‑91 (MQ=255)                                                                      acCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGa                <  1:89198/91‑1 (MQ=255)                                                                      acCGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGa                <  1:693548/91‑1 (MQ=255)                                                                        cGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAAc                          >  2:17154/1‑79 (MQ=255)                                                                        cGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAAc                          <  1:17154/79‑1 (MQ=255)                                                                        cGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATc              >  1:699297/1‑91 (MQ=255)                                                                        cGAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATc              <  2:732694/91‑1 (MQ=255)                                                                         gAAGGTGTACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATcc             >  1:751638/1‑91 (MQ=255)                                                                               gtccTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa       <  2:536464/88‑1 (MQ=255)                                                                               gtACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa       <  2:360842/91‑1 (MQ=255)                                                                               gtACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa       <  2:109870/91‑1 (MQ=255)                                                                               gtACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa       <  1:638981/91‑1 (MQ=255)                                                                               gtACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa       <  1:567535/91‑1 (MQ=255)                                                                               gtACTGCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa       <  1:121601/91‑1 (MQ=255)                                                                                    gCCAGG‑CCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCa  >  2:553097/1‑91 (MQ=255)                                                                                          |                                                                                      AAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGATGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCA  >  NZ_CP009273/1891795‑1891969

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