Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	502,414	T→C	D42G (GAT→GGT)	ybaK ←	Cys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	502,414	0	T	C	96.2%	79.0 / ‑3.3	26	D42G (GAT→GGT)	ybaK	Cys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (13/12);  total (14/12)
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.47e-01

CACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAT  >  NZ_CP009273/502356‑502486                                                           |                                                                         cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                           >  1:165838/1‑90 (MQ=255) cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                           >  2:1031670/1‑90 (MQ=255)            gtTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAt                                <  2:435925/90‑1 (MQ=255)             tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  1:117201/1‑90 (MQ=255)             tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  1:1043555/1‑90 (MQ=255)             tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  2:1100208/1‑90 (MQ=255)             tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                               >  2:888216/1‑90 (MQ=255)                  tatCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGttt                          <  1:930839/90‑1 (MQ=255)                    tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           <  1:460132/87‑1 (MQ=255)                    tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           <  1:373571/87‑1 (MQ=255)                    tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           >  2:373571/1‑87 (MQ=255)                    tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                           >  2:460132/1‑87 (MQ=255)                    tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCa                        >  1:991755/1‑90 (MQ=255)                      aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                      >  2:1081791/1‑90 (MQ=255)                       ccGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGcc                     <  2:73515/90‑1 (MQ=255)                         gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTc                                                                  >  1:823443/1‑43 (MQ=255)                         gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTc                                                                  <  2:823443/43‑1 (MQ=255)                           tCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGAt                 >  1:49152/1‑90 (MQ=255)                             aCTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCg               <  2:1121354/90‑1 (MQ=255)                              cTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGt              >  2:93609/1‑90 (MQ=255)                                    ccaGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGt        <  2:372797/90‑1 (MQ=255)                                     cagcagCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTa       <  2:387775/90‑1 (MQ=255)                                        cagcGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGt    >  1:363812/1‑90 (MQ=255)                                          gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt  <  2:311236/90‑1 (MQ=255)                                          gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt  <  2:85536/90‑1 (MQ=255)                                          gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt  <  2:953874/90‑1 (MQ=255)                                                           |                                                                         CACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAT  >  NZ_CP009273/502356‑502486

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