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	93.3%	126.7 / 2.0	45	A368A (GCT→GCC)	alaS	alanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (17/25);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAC  >  NZ_CP009273/2814185‑2814350                                                                                   |                                                                                    aaGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACt                                                                              <  2:608360/90‑1 (MQ=255)   gCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGc                                                                            <  1:1378740/90‑1 (MQ=255)     aGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctc                                                                          <  2:103204/90‑1 (MQ=255)      gTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctct                                                                         <  1:1320645/90‑1 (MQ=255)      gTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctct                                                                         <  1:1566868/90‑1 (MQ=255)             cATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcag                                                                  <  2:737827/90‑1 (MQ=255)              aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt                                                                 >  2:843219/1‑90 (MQ=255)              aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt                                                                 <  2:403632/90‑1 (MQ=255)              aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt                                                                 <  1:424101/90‑1 (MQ=255)                   gCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtcttca                                                            >  1:545159/1‑90 (MQ=255)                       aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  1:1228485/90‑1 (MQ=255)                       aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  2:124545/90‑1 (MQ=255)                       aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  1:684936/90‑1 (MQ=255)                       aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  1:327893/90‑1 (MQ=255)                        tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCttcttc                                                                    >  2:600847/1‑77 (MQ=255)                        tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCttcttc                                                                    <  1:600847/77‑1 (MQ=255)                        tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAcc                                                       >  2:943356/1‑90 (MQ=255)                            cAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCt                                                   <  2:375624/90‑1 (MQ=255)                             aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  >  2:1319464/1‑90 (MQ=255)                             aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  >  2:566668/1‑90 (MQ=255)                             aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  <  2:1544249/90‑1 (MQ=255)                             aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  <  2:276061/90‑1 (MQ=255)                             aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTtcttcag                                                                  >  1:1292335/1‑74 (MQ=255)                             aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTtcttcag                                                                  <  2:1292335/74‑1 (MQ=255)                                           gCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCctgct                                    >  2:1300267/1‑90 (MQ=255)                                            ctctTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCctgctg                                   <  1:398287/90‑1 (MQ=255)                                            ctctTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCctgctg                                   <  2:308236/90‑1 (MQ=255)                                                   tCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTc                            <  1:778365/90‑1 (MQ=255)                                                   tCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTc                            <  1:798355/90‑1 (MQ=255)                                                     cAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAg                          >  1:688005/1‑90 (MQ=255)                                                      aGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAgg                         <  2:951247/90‑1 (MQ=255)                                                       gCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGt                        >  1:1367359/1‑90 (MQ=255)                                                          aCGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCtt                     >  1:1056070/1‑90 (MQ=255)                                                               aGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCAccc                >  2:81426/1‑90 (MQ=255)                                                               aGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCAccc                >  1:276076/1‑90 (MQ=255)                                                                    gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCa                                                         >  2:795378/1‑44 (MQ=255)                                                                    gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCa                                                         <  1:795378/44‑1 (MQ=255)                                                                    gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  2:1450185/1‑90 (MQ=255)                                                                    gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  1:206883/1‑90 (MQ=255)                                                                    gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  1:269384/1‑90 (MQ=255)                                                                    gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  2:752289/1‑90 (MQ=255)                                                                    gcgcTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  1:1289922/1‑90 (MQ=255)                                                                          cAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCAt     <  2:1428839/90‑1 (MQ=255)                                                                             agTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAc  <  2:269384/90‑1 (MQ=255)                                                                             agTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAc  <  2:206883/90‑1 (MQ=255)                                                                                   |                                                                                    AAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAC  >  NZ_CP009273/2814185‑2814350

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