Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,634,151	Δ1 bp	coding (164/534 nt)	rrrQ ←	lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,634,150	0	G	.	100.0%	282.6 / NA	63	coding (165/534 nt)	rrrQ	lysozyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (22/41);  total (22/41)

CCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634062‑1634235                                                                                         |                                                                                      ccTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgc                                                                                      <  2:2241709/90‑2 (MQ=255) ccTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgc                                                                                      <  1:621469/90‑2 (MQ=255) ccTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgc                                                                                      <  2:3006855/90‑2 (MQ=255)  cTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgca                                                                                     <  1:827292/90‑3 (MQ=255)   ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    >  1:1818221/1‑87 (MQ=255)   ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    <  2:2683731/90‑4 (MQ=255)      tCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCAcccc                                                                                <  1:654920/90‑1 (MQ=255)         cGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGAc                                                                             <  2:2148830/90‑1 (MQ=255)         cGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGAc                                                                             <  1:2831706/90‑1 (MQ=255)          gatCAATGGCGTTGACTTGGTCGCATTTTTCCTTCGACAGTTTCAGATTGGGAAAAACGGTTTTTCCATCCCCCACTGT‑GCACCCCGACa                                                                            <  1:1595142/88‑1 (MQ=255)          gCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACa                                                                            >  2:2343767/1‑90 (MQ=255)           cTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAg                                                                           >  1:2213212/1‑90 (MQ=255)            tCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGa                                                                          <  2:1686573/90‑1 (MQ=255)            tCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGa                                                                          <  2:1268722/90‑1 (MQ=255)             cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  2:2731833/90‑1 (MQ=255)             cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         >  2:1471151/1‑90 (MQ=255)             cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         >  2:1380063/1‑90 (MQ=255)             cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  1:7980/90‑1 (MQ=255)             cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  2:482669/90‑1 (MQ=255)               aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         >  2:1359359/1‑88 (MQ=255)               aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  1:1359359/88‑1 (MQ=255)               aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATgg                                                                       >  1:2289649/1‑90 (MQ=255)                    gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATg                                                                        <  2:459743/84‑1 (MQ=255)                    gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATg                                                                        >  1:459743/1‑84 (MQ=255)                    gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCat                                                                  <  1:220311/90‑1 (MQ=255)                    gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCat                                                                  <  1:2451417/90‑1 (MQ=255)                    gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCat                                                                  <  2:714388/90‑1 (MQ=255)                    gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCat                                                                  <  1:1471151/90‑1 (MQ=255)                    gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCat                                                                  <  1:2681573/90‑1 (MQ=255)                     ttGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCata                                                                 <  1:870942/90‑1 (MQ=255)                       gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATg                                                               >  2:226707/1‑90 (MQ=255)                       gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATg                                                               >  2:1416653/1‑90 (MQ=255)                       gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATg                                                               >  2:1915913/1‑90 (MQ=255)                            ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          >  2:2702789/1‑90 (MQ=255)                            ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          >  2:1089761/1‑90 (MQ=255)                            ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          >  1:2700941/1‑90 (MQ=255)                                  aTTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATc                                                    <  2:2700941/90‑1 (MQ=255)                                  aTTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATc                                                    <  1:2343767/90‑1 (MQ=255)                                  aTTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATc                                                    <  2:1339263/90‑1 (MQ=255)                                  aTTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATc                                                    <  1:225048/90‑1 (MQ=255)                                    ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  1:2878267/90‑1 (MQ=255)                                    ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  2:2045469/90‑1 (MQ=255)                                      ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  2:1799594/90‑1 (MQ=255)                                      ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  1:1089761/90‑1 (MQ=255)                                      ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  2:49420/90‑1 (MQ=255)                                        ccTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTa                                              <  2:1200959/90‑1 (MQ=255)                                        ccTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTa                                              <  2:2356345/90‑1 (MQ=255)                                         cTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTAt                                             <  2:1058698/90‑1 (MQ=255)                                         cTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTAt                                             <  1:637463/90‑1 (MQ=255)                                         cTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTAt                                             <  2:2213212/90‑1 (MQ=255)                                          ttCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATg                                            >  1:1163820/1‑90 (MQ=255)                                           tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGc                                           >  2:455248/1‑90 (MQ=255)                                                cgTTTCATATTGGGAAAAACGGTTTTTCCATCCACCGCTGT‑GCACCCCGCCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtg                                      <  1:2219082/89‑1 (MQ=255)                                                 gTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtgt                                     >  2:2405003/1‑90 (MQ=255)                                                        atTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTa                              <  1:2077982/90‑1 (MQ=255)                                                        atTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTa                              <  1:619562/90‑1 (MQ=255)                                                           gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGtt                               <  1:2906659/86‑1 (MQ=255)                                                           gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGtt                               >  2:2906659/1‑86 (MQ=255)                                                           gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  1:417947/1‑90 (MQ=255)                                                           gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  1:2226842/1‑90 (MQ=255)                                                           gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  2:292875/1‑90 (MQ=255)                                                                aaaCGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCttt                      >  2:798752/1‑90 (MQ=255)                                                                             atccaCCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc         <  1:798752/90‑1 (MQ=255)                                                                             atccaCCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc         <  1:1951390/90‑1 (MQ=255)                                                                             atccaCCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc         >  2:2190077/1‑90 (MQ=255)                                                                               ccaccaCTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTg       <  1:226707/90‑1 (MQ=255)                                                                                  ccacTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTc    <  2:2061396/90‑1 (MQ=255)                                                                                  ccacTGT‑GCACCCCGACAGATGGACCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTCTTTCGTCCAGAACCTGGTc    <  1:2057978/90‑1 (MQ=255)                                                                                    acTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCaa  <  1:39817/90‑1 (MQ=255)                                                                                         |                                                                                      CCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634062‑1634235

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