Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,568,198	A→G	F417S (TTC→TCC)	pqqL ←	insulinase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,568,198	0	A	G	100.0%	55.4 / NA	18	F417S (TTC→TCC)	pqqL	insulinase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (4/14);  total (4/14)

AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCA  >  NZ_CP009273/1568111‑1568272                                                                                        |                                                                           aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                          <  1:31975/90‑1 (MQ=255) aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                          <  2:422879/90‑1 (MQ=255)   ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                        <  1:567838/90‑1 (MQ=255)          cATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGtatta                                                                 <  1:190446/90‑1 (MQ=255)                     cAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCt                                                      <  1:172712/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                   <  1:582142/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                   <  2:72195/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                   <  2:529321/90‑1 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                               >  2:329218/1‑90 (MQ=255)                                             ccACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCAt                              <  1:277398/90‑1 (MQ=255)                                                     gTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAAt                      <  1:474934/90‑1 (MQ=255)                                                      tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATc                     >  1:418821/1‑90 (MQ=255)                                                      tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATc                     >  2:58691/1‑90 (MQ=255)                                                            aaGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTc               <  2:79267/90‑1 (MQ=255)                                                            aaGTTGACATGCCCGTTCCGGCCACAAGGACGGGGTATTATTTAATGAGCTGGATGCAAGGCGACTGGTCAGCATACGTAACTCACGTTc               <  2:260756/90‑1 (MQ=255)                                                              gTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAg             >  2:129163/1‑90 (MQ=255)                                                                 gATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCtt          <  2:418821/90‑1 (MQ=255)                                                                         tctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCa  <  2:594565/90‑1 (MQ=255)                                                                                        |                                                                           AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCA  >  NZ_CP009273/1568111‑1568272

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