Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,505,430	A→G	K391E (AAA→GAA)	ydcR →	PLP‑dependent aminotransferase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,505,430	0	A	G	95.5%	68.0 / NA	22	K391E (AAA→GAA)	ydcR	PLP‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/6);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGAC  >  NZ_CP009273/1505343‑1505518                                                                                        |                                                                                         gCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGCAACGTAAACAACGGGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGgaa                                                                                        >  1:503306/1‑90 (MQ=255)    cATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAAtt                                                                                     >  1:400502/1‑90 (MQ=255)    cATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAAtt                                                                                     >  1:436598/1‑90 (MQ=255)     aTCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATtc                                                                                    >  1:281184/1‑90 (MQ=255)                 tGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACa                                                                        <  1:677565/90‑1 (MQ=255)                  gCGTCGCCAGCTTGCGGAACGGAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAg                                                                       >  1:712577/1‑90 (MQ=255)                     tCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTgg                                                                    <  1:680000/90‑1 (MQ=255)                           gCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACtt                                                              <  2:436598/90‑1 (MQ=255)                                              cGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTcccc                                           >  2:384912/1‑90 (MQ=255)                                                tGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGa                                         >  1:585485/1‑90 (MQ=255)                                                    tGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCg                                     >  2:436422/1‑90 (MQ=255)                                                               gctgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                          >  2:554833/1‑90 (MQ=255)                                                               gctgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                          >  2:313454/1‑90 (MQ=255)                                                               gctgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                          >  2:294293/1‑90 (MQ=255)                                                                 tgcTTTATCTGCCTGGGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCg                        <  2:257288/90‑1 (MQ=255)                                                                 tgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCg                        <  2:105782/90‑1 (MQ=255)                                                                       aTCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCTGGCGAATTaa                  >  2:4844/1‑90 (MQ=255)                                                                       aTCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTaa                  >  2:493568/1‑90 (MQ=255)                                                                       aTCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTaa                  >  2:608854/1‑90 (MQ=255)                                                                          tgcctgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGcc               >  2:52269/1‑90 (MQ=255)                                                                              tgcGGAAGTGGAAATTCATCATAATGACAGGGGTTACTTTCTCGGGTTGGAGCCCCCGGAGCCGTTAGATGCCGGCGAATTAAGCCTggc           <  1:750269/90‑1 (MQ=255)                                                                               gcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggcg          >  2:411910/1‑90 (MQ=255)                                                                                       gAAAATTCATCATAATGACAGTGTTTATTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGAc  <  1:384912/90‑1 (MQ=255)                                                                                        |                                                                                         GCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGAC  >  NZ_CP009273/1505343‑1505518

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