Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,189,312	A→G	T512A (ACC→GCC)	metG →	methionine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,189,312	0	A	G	98.0%	157.3 / ‑3.6	50	T512A (ACC→GCC)	metG	methionine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/32);  total (17/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

ACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGG  >  NZ_CP009273/2189226‑2189392                                                                                       |                                                                                 aCCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCt                                                                               <  1:253286/90‑1 (MQ=255) aCCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCt                                                                               <  1:1140619/90‑1 (MQ=255) aCCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCt                                                                               >  1:185041/1‑90 (MQ=255) aCCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCt                                                                               <  2:1555168/90‑1 (MQ=255)      ttCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGAt                                                                          <  1:1023249/90‑1 (MQ=255)      ttCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGAt                                                                          <  2:263966/90‑1 (MQ=255)          gcgTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTa                                                                      >  2:1405120/1‑90 (MQ=255)               cTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAg                                                                 <  1:1133526/90‑1 (MQ=255)                         aCCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCgctgc                                                       >  2:436628/1‑90 (MQ=255)                         aCCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCgctgc                                                       >  2:566302/1‑90 (MQ=255)                         aCCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCgctgc                                                       >  1:987628/1‑90 (MQ=255)                             gAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTggg                                                   >  1:640463/1‑90 (MQ=255)                                 ccGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCcac                                               <  2:1014397/90‑1 (MQ=255)                                   ggTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACaa                                             <  2:343646/90‑1 (MQ=255)                                   ggTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACaa                                             <  1:901336/90‑1 (MQ=255)                                     tACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAg                                           <  2:385201/90‑1 (MQ=255)                                         gCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGAt                                                                          <  1:952853/55‑1 (MQ=255)                                         gCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGAt                                                                          >  2:952853/1‑55 (MQ=255)                                         gCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGaa                                       >  1:1018204/1‑90 (MQ=255)                                           cGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATc                                     <  1:1148927/90‑1 (MQ=255)                                           cGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATc                                     <  1:1168366/90‑1 (MQ=255)                                              aaCTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGt                                  >  1:893182/1‑90 (MQ=255)                                                  gACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCaa                              >  1:1537302/1‑90 (MQ=255)                                                  gACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCaa                              >  1:1370644/1‑90 (MQ=255)                                                   aCCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAg                             >  1:1283774/1‑90 (MQ=255)                                                     cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                           <  1:436628/90‑1 (MQ=255)                                                     cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                           <  2:185041/90‑1 (MQ=255)                                                     cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                           <  2:281896/90‑1 (MQ=255)                                                     cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                           <  1:1358506/90‑1 (MQ=255)                                                     cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                           <  1:876658/90‑1 (MQ=255)                                                     cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                           <  1:782987/90‑1 (MQ=255)                                                     cGAGCGGGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                           <  1:324120/90‑1 (MQ=255)                                                      gAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcg                          >  1:872057/1‑90 (MQ=255)                                                       aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                         <  2:893182/90‑1 (MQ=255)                                                       aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                         <  2:1512135/90‑1 (MQ=255)                                                       aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                         <  2:1336650/90‑1 (MQ=255)                                                       aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                         <  2:1283774/90‑1 (MQ=255)                                                       aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                         <  2:1041862/90‑1 (MQ=255)                                                           tGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGta                     <  2:875675/90‑1 (MQ=255)                                                                 aGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCg               >  1:159874/1‑90 (MQ=255)                                                                    aTTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCAt            <  2:1210588/90‑1 (MQ=255)                                                                    aTTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCAt            >  1:1141523/1‑90 (MQ=255)                                                                     ttCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATc           <  2:1055764/90‑1 (MQ=255)                                                                     ttCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATc           <  2:64602/90‑1 (MQ=255)                                                                       ccTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTgg                                                    >  2:1524805/1‑47 (MQ=255)                                                                       ccTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTgg                                                    <  1:1524805/47‑1 (MQ=255)                                                                          cAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGatat      <  2:957362/90‑1 (MQ=255)                                                                              aCGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAgg  <  1:254307/90‑1 (MQ=255)                                                                                 gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATc                                     >  2:80970/1‑52 (MQ=255)                                                                                 gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATc                                     <  1:80970/52‑1 (MQ=255)                                                                                       |                                                                                 ACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGG  >  NZ_CP009273/2189226‑2189392

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: ‑