Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,796,543	G→C	G95G (GGC→GGG)	thrS ←	threonine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,796,543	0	G	C	76.7%	61.8 / 13.6	30	G95G (GGC→GGG)	thrS	threonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base C (7/16);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCT  >  NZ_CP009273/1796483‑1796630                                                             |                                                                                        ttCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAg                                                            >  1:153551/1‑90 (MQ=255)           aaCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATc                                                  >  2:837591/1‑90 (MQ=255)           aaCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATc                                                  >  1:135188/1‑90 (MQ=255)           aaCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATc                                                  >  1:520323/1‑90 (MQ=255)            aCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcg                                                 >  2:146137/1‑90 (MQ=255)            aCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcg                                                 >  1:370069/1‑90 (MQ=255)             cGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcgc                                                <  2:143983/90‑1 (MQ=255)             cGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcgc                                                <  1:582257/90‑1 (MQ=255)              gTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcgcg                                               <  1:520401/90‑1 (MQ=255)               tGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGt                                              <  1:61774/90‑1 (MQ=255)                gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCGAAGTTGTTTAATCGCGTg                                             <  2:252581/90‑1 (MQ=255)                gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  1:478104/90‑1 (MQ=255)                gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  2:667048/90‑1 (MQ=255)                gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  2:627845/90‑1 (MQ=255)                gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  2:436858/90‑1 (MQ=255)                  ggTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGcc                                           >  1:319409/1‑90 (MQ=255)                          tCAACGTCGTAATAAAAACCGTTGTCAATACCCGGGCCGATTGCCGTTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAgg                                   <  2:690106/90‑1 (MQ=255)                             aCGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTgc                                >  1:696988/1‑90 (MQ=255)                             aCGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTg                                 >  1:722067/1‑89 (MQ=255)                             aCGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTg                                 <  2:722067/89‑1 (MQ=255)                               gtcgtAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTgcgc                              <  1:408561/90‑1 (MQ=255)                                   taataaAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAg                          <  2:319409/90‑1 (MQ=255)                                   taataaAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAg                          <  2:422801/90‑1 (MQ=255)                                           aCCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACg                  <  2:54855/90‑1 (MQ=255)                                              gTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAAt               >  1:657803/1‑90 (MQ=255)                                               ttGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATg              <  1:835268/90‑1 (MQ=255)                                               ttGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATg              <  1:640316/90‑1 (MQ=255)                                                tGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGa             <  2:426653/90‑1 (MQ=255)                                                         ccGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGAc    <  2:455673/90‑1 (MQ=255)                                                         ccGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGAc    <  2:716063/90‑1 (MQ=255)                                                           gggCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACct  >  1:246412/1‑90 (MQ=255)                                                             |                                                                                        TTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCT  >  NZ_CP009273/1796483‑1796630

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