Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,095,787	A→G	F142S (TTC→TCC)	trmB ←	tRNA (guanosine(46)‑N7)‑methyltransferase TrmB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,095,787	0	A	G	95.7%	72.9 / ‑3.9	23	F142S (TTC→TCC)	trmB	tRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (9/13);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

GCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGT  >  NZ_CP009273/3095701‑3095871                                                                                       |                                                                                     gCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGgaaa                                                                                   <  2:522834/90‑1 (MQ=255) gCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGgaaa                                                                                   <  1:898259/90‑1 (MQ=255)    gCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAaga                                                                                >  2:228067/1‑90 (MQ=255)      aGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGc                                                                              >  1:850196/1‑90 (MQ=255)      aGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGGTTATTATGGCGCGCTTTGTGCCACGGGGCAGGGGAAAAGAGc                                                                              >  1:1016610/1‑90 (MQ=255)            cTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTTCCACGGGTCAGGGGAAAAGAGCTGCAcc                                                                        <  1:779610/90‑1 (MQ=255)             tttttACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCa                                                                       <  1:1071650/90‑1 (MQ=255)              ttttACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCAt                                                                      <  1:282886/90‑1 (MQ=255)               tttACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATg                                                                     <  2:793741/90‑1 (MQ=255)                         cGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAAt                                                           <  1:76453/90‑1 (MQ=255)                               aCGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATTCGCAATGAATTGTCAg                                                     <  2:68571/90‑1 (MQ=255)                                   cACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCCCGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAAt                                                 <  2:719407/90‑1 (MQ=255)                                             aTACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATgc                                       >  2:384582/1‑90 (MQ=255)                                                           atGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCga                         <  1:472232/90‑2 (MQ=255)                                                             ggCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCAt                       <  1:232097/90‑1 (MQ=255)                                                             ggCGCGCTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCAt                       <  2:84360/90‑1 (MQ=255)                                                                gcgcTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGt                    >  1:646644/1‑90 (MQ=255)                                                                  gcTTTGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGa                  >  1:519136/1‑90 (MQ=255)                                                                     ttGTGCCACGGGTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGacac               >  1:678856/1‑90 (MQ=255)                                                                               ggTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCa     >  2:133676/1‑90 (MQ=255)                                                                               ggTCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCa     <  1:727250/90‑1 (MQ=255)                                                                                 tCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAAcc                           <  1:761642/66‑1 (MQ=255)                                                                                 tCAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAAcc                           >  2:761642/1‑66 (MQ=255)                                                                                  cAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGt  <  1:171937/90‑1 (MQ=255)                                                                                  cAGGGGAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGt  <  2:678856/90‑1 (MQ=255)                                                                                       |                                                                                     GCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGT  >  NZ_CP009273/3095701‑3095871

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