Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,814,267	A→G	A368A (GCT→GCC)	alaS ←	alanine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,814,267	0	A	G	84.0%	51.1 / 5.9	25	A368A (GCT→GCC)	alaS	alanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (10/11);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-02

CGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGT  >  NZ_CP009273/2814181‑2814352                                                                                       |                                                                                      cGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCa                                                                                    >  1:401062/1‑90 (MQ=255)  ggAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCaa                                                                                   <  1:372211/90‑1 (MQ=255)         aGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctc                                                                            <  1:14428/90‑1 (MQ=255)                 cATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcag                                                                    <  1:3308/90‑1 (MQ=255)                 cATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcag                                                                    <  1:382642/90‑1 (MQ=255)                  aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt                                                                   >  2:57314/1‑90 (MQ=255)                    ccAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtct                                                                 >  1:453002/1‑90 (MQ=255)                       gCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtcttca                                                              >  2:58525/1‑90 (MQ=255)                       gCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtcttca                                                              >  1:700076/1‑90 (MQ=255)                          tATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCa                                                           <  1:623815/90‑1 (MQ=255)                           aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                          <  1:58525/90‑1 (MQ=255)                            tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAcc                                                         >  2:347389/1‑90 (MQ=255)                                 aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                    >  1:103417/1‑90 (MQ=255)                                                  ctTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGc                                   <  2:654739/90‑1 (MQ=255)                                                         cAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAg                            >  1:321303/1‑90 (MQ=255)                                                             aaCGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCt                        <  2:103417/90‑1 (MQ=255)                                                             aaCGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCt                        <  2:423389/90‑1 (MQ=255)                                                                gCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCa                     >  1:362731/1‑90 (MQ=255)                                                                       cgcgCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCag              <  2:378903/90‑1 (MQ=255)                                                                       cgcgCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCag              <  2:453002/90‑1 (MQ=255)                                                                       cgcgCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCag              <  2:641048/90‑1 (MQ=255)                                                                        gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga             >  2:581081/1‑90 (MQ=255)                                                                        gcgcTCCAGAGTACGAGCAAACTGCTCTtcttcagt                                                                   >  1:534528/1‑36 (MQ=255)                                                                        gcgcTCCAGAGTACGAGCAAACTGCTCTtcttcagt                                                                   <  2:534528/36‑1 (MQ=255)                                                                                   tACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGt  >  1:152216/1‑90 (MQ=255)                                                                                       |                                                                                      CGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGT  >  NZ_CP009273/2814181‑2814352

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