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%	82.1 / NA	24	F417S (TTC→TCC)	pqqL	insulinase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/13);  total (11/13)

AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAC  >  NZ_CP009273/1568111‑1568274                                                                                        |                                                                             aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                            <  1:416707/90‑1 (MQ=255)   ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                          <  1:108448/90‑1 (MQ=255)   ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                          <  1:286322/90‑1 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                          >  1:127527/1‑88 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                          <  2:127527/88‑1 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg                                                                        >  1:883031/1‑90 (MQ=255)      ctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGt                                                                       >  1:616428/1‑90 (MQ=255)       tgcCATTTTTCCGCCAGACTTTGCGCGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGta                                                                      <  2:163539/90‑1 (MQ=255)        gcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGtat                                                                     <  1:380100/90‑1 (MQ=255)        gcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGtat                                                                     <  2:211534/90‑1 (MQ=255)          cATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGtatta                                                                   <  1:102973/90‑1 (MQ=255)                  cGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGa                                                           <  1:688034/90‑1 (MQ=255)                                  gTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGa                                           >  2:908343/1‑90 (MQ=255)                                      tttGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTgg                                       <  1:298650/90‑1 (MQ=255)                                           tgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGc                                  <  1:143877/90‑1 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                            <  1:171460/47‑1 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                            >  2:171460/1‑47 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                                 >  1:330150/1‑90 (MQ=255)                                                    cGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTg                                                       >  1:684793/1‑60 (MQ=255)                                                    cGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTg                                                       <  2:684793/60‑1 (MQ=255)                                                       ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCa                      >  2:45365/1‑90 (MQ=255)                                                                        ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc     >  1:168951/1‑90 (MQ=255)                                                                        ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc     >  1:482650/1‑90 (MQ=255)                                                                           ttCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAc  >  2:841647/1‑90 (MQ=255)                                                                                        |                                                                             AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAC  >  NZ_CP009273/1568111‑1568274

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