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	.	97.6%	168.8 / ‑1.5	42	coding (165/534 nt)	rrrQ	lysozyme
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (15/26);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

TTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634064‑1634235                                                                                       |                                                                                      ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    <  2:795630/90‑4 (MQ=255)       cGCTCAATGGCGTTGACCTGGGCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCggc                                                                             >  1:1270096/1‑88 (MQ=255)         cTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACa                                                                            <  2:1246966/90‑1 (MQ=255)           cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  1:1330015/90‑1 (MQ=255)           cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  1:1505490/90‑1 (MQ=255)           cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  2:692418/90‑1 (MQ=255)             aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATgg                                                                       >  2:301163/1‑90 (MQ=255)                  gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCat                                                                  <  2:1032909/90‑1 (MQ=255)                  gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCat                                                                  >  2:633795/1‑90 (MQ=255)                   ttGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCata                                                                 <  2:400606/90‑1 (MQ=255)                   ttGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCata                                                                 <  2:741924/90‑1 (MQ=255)                          ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          >  1:1470606/1‑90 (MQ=255)                          ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          <  2:155226/90‑1 (MQ=255)                          ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          >  1:379401/1‑90 (MQ=255)                               cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         >  2:878758/1‑70 (MQ=255)                               cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         >  2:1245395/1‑70 (MQ=255)                               cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  1:1245395/70‑1 (MQ=255)                               cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGAt                                                                         <  1:878758/70‑1 (MQ=255)                               cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATgg                                                                       >  1:1498058/1‑72 (MQ=255)                               cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATgg                                                                       <  2:1498058/72‑1 (MQ=255)                               cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCAt                                                     >  1:470085/1‑90 (MQ=255)                                  ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGCCCATATGCCAGAACCATcgc                                                  <  2:174837/90‑1 (MQ=255)                                  tcttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  1:760036/88‑1 (MQ=255)                                    ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  2:579762/90‑1 (MQ=255)                                      ccTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTa                                              <  2:538305/90‑1 (MQ=255)                                      ccTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTa                                              <  2:396412/90‑1 (MQ=255)                                         tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGc                                           <  2:419907/90‑1 (MQ=255)                                                    atatTGGGAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGt                                <  2:1126523/90‑1 (MQ=255)                                                         gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                >  2:577676/1‑69 (MQ=255)                                                         gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  1:577676/69‑1 (MQ=255)                                                         gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtgt                                     >  1:247740/1‑80 (MQ=255)                                                         gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtgt                                     <  2:247740/80‑1 (MQ=255)                                                         gggAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  1:110282/1‑90 (MQ=255)                                                           gAAAAACGGTTTTTCCATCCCCCACTGT‑GCCCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACcttc                         <  1:1133889/90‑1 (MQ=255)                                                           gAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          >  2:1086777/1‑57 (MQ=255)                                                           gAAAAACGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGa                                                          <  1:1086777/57‑1 (MQ=255)                                                                 cGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTc                   <  1:889413/90‑1 (MQ=255)                                                                 cGGTTTTTCCATCCACCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTc                   <  2:978981/90‑1 (MQ=255)                                                                           atccaCCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc         >  2:1076167/1‑90 (MQ=255)                                                                           atccaCCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc         >  2:721614/1‑90 (MQ=255)                                                                            tccaCCACTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACt        <  1:470676/90‑1 (MQ=255)                                                                                  acTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCaa  <  2:1470606/90‑1 (MQ=255)                                                                                  acTGT‑GCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCaa  <  1:1417647/90‑1 (MQ=255)                                                                                       |                                                                                      TTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634064‑1634235

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