Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	412,148	G→A	V87V (GTC→GTT)	sbcD ←	exonuclease subunit SbcD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	412,148	0	G	A	100.0%	164.1 / NA	42	V87V (GTC→GTT)	sbcD	exonuclease subunit SbcD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/22);  total (20/22)

AAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAA  >  NZ_CP009273/412016‑412281                                                                                                                                     |                                                                                                                                      aaaCGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAg                                                                                                                                 <  2:487159/139‑1 (MQ=255)           gggCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTc                                                                                                                       <  2:41242/139‑1 (MQ=255)              cacaGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTg                                                                                                                    <  1:318849/139‑1 (MQ=255)                 aGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGc                                                                                                                   >  2:357265/1‑137 (MQ=255)                    aCTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTa                                                                                                              <  2:492478/139‑1 (MQ=255)                        cgcCTGGCGTCCCGTCGCGACGAGGCAAGCTTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTaccac                                                                                                          >  2:471936/1‑139 (MQ=255)                        cgcCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTaccac                                                                                                          >  1:15367/1‑139 (MQ=255)                        cgcCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTaccac                                                                                                          >  2:206152/1‑139 (MQ=255)                           cTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTac                                                                                                             <  1:257233/133‑1 (MQ=255)                           cTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTac                                                                                                             >  2:257233/1‑133 (MQ=255)                                  cccGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAg                                                                                                >  2:771315/1‑139 (MQ=255)                                        gcgACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGt                                                                                           <  2:15367/138‑1 (MQ=255)                                                          cggcgCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGacaac                                                                         >  1:81959/1‑138 (MQ=255)                                                             cgcATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACaaaa                                                                     <  1:98668/139‑1 (MQ=255)                                                                aTGTCCGGCTCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACgg                                                                  <  1:315437/139‑1 (MQ=255)                                                                       gcgcTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGtata                                                           <  1:771315/139‑1 (MQ=255)                                                                                aCCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCggg                                                  <  1:206152/139‑1 (MQ=255)                                                                                              gAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATACCGTGCGGGCGTAACTGggcgg                                     <  1:423777/138‑1 (MQ=255)                                                                                                 gAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACg                                                                   >  2:841009/1‑105 (MQ=255)                                                                                                 gAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACg                                                                   <  1:841009/105‑1 (MQ=255)                                                                                                 gAACGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAg                                 >  1:188986/1‑139 (MQ=255)                                                                                                  aaCGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGc                                <  2:718958/139‑1 (MQ=255)                                                                                                    cGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACaaaa                                                                     >  2:73421/1‑100 (MQ=255)                                                                                                    cGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACaaaa                                                                     <  1:73421/100‑1 (MQ=255)                                                                                                    cGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACg                                                                   >  1:380048/1‑102 (MQ=255)                                                                                                    cGCCATGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACg                                                                   <  2:380048/102‑1 (MQ=255)                                                                                                         tGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACg                                                                   <  1:128053/97‑1 (MQ=255)                                                                                                         tGATATCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACg                                                                   >  2:128053/1‑97 (MQ=255)                                                                                                            tatCGCGCGATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGaa                      <  2:188986/139‑1 (MQ=255)                                                                                                               cgcgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGt                                                       >  1:265286/1‑103 (MQ=255)                                                                                                               cgcgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGt                                                       <  2:265286/103‑1 (MQ=255)                                                                                                               cgcgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGaaaa                    >  2:453658/1‑138 (MQ=255)                                                                                                               cgcgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGaaaa                    >  2:560765/1‑138 (MQ=255)                                                                                                                 cgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCa                                                                                                                 >  2:348114/1‑43 (MQ=255)                                                                                                                 cgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCa                                                                                                                 <  1:348114/43‑1 (MQ=255)                                                                                                                 cgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGggcgg                                     <  1:653400/119‑1 (MQ=255)                                                                                                                 cgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGggcgg                                     >  2:653400/1‑119 (MQ=255)                                                                                                                 cgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACAt                 >  1:529044/1‑139 (MQ=255)                                                                                                                 cgcgATTCATTCAGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACAt                 >  2:6864/1‑139 (MQ=255)                                                                                                                         attcaGCGTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGc          <  2:130005/138‑1 (MQ=255)                                                                                                                                gTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACaataa  >  2:65048/1‑139 (MQ=255)                                                                                                                                gTGGCAACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACaataa  <  2:714548/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      AAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAA  >  NZ_CP009273/412016‑412281

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