Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,384,689	Δ1 bp	coding (699/1446 nt)	tldD ←	metalloprotease TldD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,384,688	0	C	.	100.0%	131.0 / NA	31	coding (700/1446 nt)	tldD	metalloprotease TldD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (16/15);  total (16/15)

TCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGA  >  NZ_CP009273/3384600‑3384753                                                                                         |                                                                  tCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAgcgca                                                                  <  2:519038/90‑2 (MQ=255)        aCAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCAc                                                          >  2:415523/1‑90 (MQ=255)         cAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACt                                                         <  1:123088/90‑1 (MQ=255)           gCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGc                                                       <  2:963657/90‑1 (MQ=255)           gCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGc                                                       >  2:231862/1‑90 (MQ=255)                cccGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGcttctt                                                  >  1:233422/1‑90 (MQ=255)                   ggCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTg                                               >  1:132902/1‑90 (MQ=255)                        aCCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCAt                                          <  2:641567/90‑1 (MQ=255)                        aCCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCAt                                          <  2:185137/90‑1 (MQ=255)                          cTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGc                                        <  1:136305/90‑1 (MQ=255)                            gcgcCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCAt                                      >  1:122234/1‑90 (MQ=255)                            gcgcCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCAt                                      >  2:767963/1‑90 (MQ=255)                                  aGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCAc                                <  1:433214/90‑1 (MQ=255)                                  aGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCAc                                <  1:720325/90‑1 (MQ=255)                                      ctacCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCTCGGAc                            <  2:495270/90‑1 (MQ=255)                                        acCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACtt                          >  1:667929/1‑90 (MQ=255)                                           ggCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACtt                          <  2:924917/87‑1 (MQ=255)                                           ggCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACtt                          >  1:924917/1‑87 (MQ=255)                                           ggCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGc                       <  1:767963/90‑1 (MQ=255)                                           ggCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGc                       >  1:104681/1‑90 (MQ=255)                                           ggCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGc                       >  1:662346/1‑90 (MQ=255)                                             cATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCg                     >  2:255563/1‑90 (MQ=255)                                             cATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCg                     >  2:244000/1‑90 (MQ=255)                                                     ccGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATc             <  1:974157/90‑1 (MQ=255)                                                           gTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAg       >  2:258118/1‑90 (MQ=255)                                                           gTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAg       >  1:836177/1‑90 (MQ=255)                                                           gTGCAGCAACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAg       >  2:951978/1‑90 (MQ=255)                                                               agcaACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGgaag   <  1:762709/90‑1 (MQ=255)                                                                gcaACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGgaaga  <  2:132902/90‑1 (MQ=255)                                                                gcaACGGCAGAAAGATTGACCAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGgaaga  <  1:269701/90‑1 (MQ=255)                                                                                   ccAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATc                                     >  2:514683/1‑36 (MQ=38)                                                                                   ccAGCG‑CATACGCACTGCTTCTTTTGCCCATGCATc                                     <  1:514683/36‑1 (MQ=38)                                                                                         |                                                                  TCATGCAACAGCACGCCCGGCCAACCTGCGCCAAGTACTACCGGCATGGTGCCCGCTGGTGCAGCAACGGCAGAAAGATTGACCAGCGCCATACGCACTGCTTCTTTTGCCCATGCATCCGCACGGACTTCGCCGTCGAGATCGGCAAGGAAGA  >  NZ_CP009273/3384600‑3384753

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