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

AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCAT  >  NZ_CP009273/1568111‑1568279                                                                                        |                                                                                  aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 <  2:2037903/90‑1 (MQ=255) aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 <  1:2799871/90‑1 (MQ=255) aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 <  2:1685216/90‑1 (MQ=255) aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 <  2:307906/90‑1 (MQ=255) aaCTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 <  1:1192695/90‑1 (MQ=255)   ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                               <  1:518064/90‑1 (MQ=255)   ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                               <  1:1837855/90‑1 (MQ=255)   ctgctgCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACgg                                                                               <  2:1662900/90‑1 (MQ=255)    tgctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGc                                                                              >  2:2943406/1‑90 (MQ=255)    tgctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGc                                                                              >  1:1391294/1‑90 (MQ=255)    tgctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGc                                                                              >  1:1518041/1‑90 (MQ=255)    tgctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGc                                                                              >  1:524490/1‑90 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg                                                                             >  1:1149115/1‑90 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg                                                                             >  1:499889/1‑90 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg                                                                             <  2:1149115/90‑1 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg                                                                             >  1:2284192/1‑90 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg                                                                             >  2:2399401/1‑90 (MQ=255)     gctgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCg                                                                             >  2:2994615/1‑90 (MQ=255)       tgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGta                                                                           <  2:2466923/90‑1 (MQ=255)       tgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGta                                                                           <  1:2473643/90‑1 (MQ=255)       tgcCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGta                                                                           <  1:2485728/90‑1 (MQ=255)          cATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGtatta                                                                        <  2:770641/90‑1 (MQ=255)               ttCCGCCAGACTTTTCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTaa                                                                   <  1:2628597/90‑1 (MQ=255)                     cgtaCTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCt                                                             <  1:2299571/87‑1 (MQ=255)                     cAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCt                                                             <  1:475878/90‑1 (MQ=255)                         cTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGAt                                                         <  2:2284192/90‑1 (MQ=255)                         cTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGAt                                                         <  1:2503741/90‑1 (MQ=255)                           ttGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGc                                                       >  1:2342060/1‑90 (MQ=255)                              cACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAg                                                    <  2:2251120/90‑1 (MQ=255)                              cACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAg                                                    <  1:2120962/90‑1 (MQ=255)                                       ttGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGt                                           <  2:2342060/90‑1 (MQ=255)                                       ttGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGt                                           <  2:693859/90‑1 (MQ=255)                                       ttGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGt                                           <  2:2967062/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                          <  2:507480/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                          <  1:1343093/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                          <  2:43751/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                          <  2:524490/90‑1 (MQ=255)                                        tgctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTc                                          <  1:2399401/90‑1 (MQ=255)                                          ctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 >  1:2385033/1‑49 (MQ=255)                                          ctgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 <  2:2385033/49‑1 (MQ=255)                                           tgcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGc                                       >  1:2321126/1‑90 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTg                                                            >  2:259189/1‑68 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTg                                                            <  1:259189/68‑1 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                                      >  1:1234375/1‑90 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                                      >  2:743615/1‑90 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                                      >  2:444827/1‑90 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                                      >  1:1139568/1‑90 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                                      >  1:465346/1‑90 (MQ=255)                                            gcCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCa                                      >  1:260527/1‑90 (MQ=255)                                             ccACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCAt                                     <  2:15776/90‑1 (MQ=255)                                                     gTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAAt                             <  2:1139568/90‑1 (MQ=255)                                                     gTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAAt                             <  2:1458837/90‑1 (MQ=255)                                                      tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGAc                                                                                 <  1:1897591/37‑1 (MQ=255)                                                      tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATc                            <  2:1912719/90‑1 (MQ=255)                                                      tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATc                            >  2:1747291/1‑90 (MQ=255)                                                       ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCa                           >  2:2214861/1‑90 (MQ=255)                                                            aaGTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTc                      <  2:35424/90‑1 (MQ=255)                                                              gTTGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAg                    <  1:2994615/90‑1 (MQ=255)                                                                tGATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCt                  >  1:462367/1‑90 (MQ=255)                                                                 gATATGTCTCTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCtt                 <  1:1516762/90‑1 (MQ=255)                                                                       tctcTTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGAt           >  2:2859470/1‑90 (MQ=255)                                                                        ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc          >  2:2763466/1‑90 (MQ=255)                                                                        ctctTCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATc          >  1:2747277/1‑90 (MQ=255)                                                                           ttCCGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAc       >  1:1221949/1‑90 (MQ=255)                                                                              cGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGc    <  2:2168736/90‑1 (MQ=255)                                                                              cGGCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGc    <  2:1467679/90‑1 (MQ=255)                                                                                gCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCAt  >  1:700068/1‑90 (MQ=255)                                                                                gCGACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCAt  >  1:1646364/1‑90 (MQ=255)                                                                                  gACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAc       <  2:2514445/83‑1 (MQ=255)                                                                                  gACAAGGACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAAc       >  1:2514445/1‑83 (MQ=255)                                                                                        |                                                                                  AACTGCTGCCATTTTTCCGCCAGACTTTGCACGGTAATTTGCTGCCACAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCAT  >  NZ_CP009273/1568111‑1568279

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