Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,113,649	T→G	K330N (AAA→AAC)	glcA ←	glycolate permease GlcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,113,649	0	T	G	86.5%	103.0 / 8.2	37	K330N (AAA→AAC)	glcA	glycolate permease GlcA
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (8/24);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC  >  NZ_CP009273/3113564‑3113731                                                                                      |                                                                                   ccGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaa                                                                                >  1:497081/1‑90 (MQ=255) ccGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaa                                                                                <  2:256846/90‑1 (MQ=255)   gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAAc                                                                              <  1:350103/90‑1 (MQ=255)   gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAAc                                                                              >  2:66179/1‑90 (MQ=255)   gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAAc                                                                              <  2:271663/90‑1 (MQ=255)           cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                      <  1:378347/90‑1 (MQ=255)           cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                      <  2:244528/90‑1 (MQ=255)            acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                     <  1:349793/90‑1 (MQ=255)            acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                     <  1:67278/90‑1 (MQ=255)                gttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTcc                                                                 >  2:359061/1‑90 (MQ=255)                 ttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCa                                                                <  1:432988/90‑1 (MQ=255)                    aTGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGAt                                                             >  1:477079/1‑90 (MQ=255)                     tGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATg                                                            >  2:718171/1‑90 (MQ=255)                      gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                           <  2:6244/90‑1 (MQ=255)                      gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                           <  2:330364/90‑1 (MQ=255)                      gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                           <  2:346213/90‑1 (MQ=255)                      gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                           <  1:630562/90‑1 (MQ=255)                       cAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGAt                                                             <  2:280836/87‑1 (MQ=255)                       cAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGAt                                                             >  1:280836/1‑87 (MQ=255)                           tGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGACGCCCGGCGCCAATAACGCGTTAAACGGCCTCATGGTCCAGCTGGTGAcc                                                      <  2:504968/90‑1 (MQ=255)                              gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                   >  2:502925/1‑90 (MQ=255)                              gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                   <  1:118423/90‑1 (MQ=255)                                gATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCAc                                                 <  2:323113/90‑1 (MQ=255)                                 aTCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCAcc                                                >  2:219216/1‑90 (MQ=255)                                     ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                            <  1:698898/90‑1 (MQ=255)                                          ttAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGAtt                                       >  1:9064/1‑90 (MQ=255)                                             aTCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaaa                                    >  2:721840/1‑90 (MQ=255)                                                  cAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACgg                               <  1:708786/90‑1 (MQ=255)                                                         tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGTCCAGCACCGTTAAGATTAAAAACGGTGACCAc                        <  1:701185/90‑1 (MQ=255)                                                         tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCAc                        <  1:721840/90‑1 (MQ=255)                                                                cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                 >  1:571054/1‑90 (MQ=255)                                                                      ccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAttt           <  2:53311/90‑1 (MQ=255)                                                                          aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc       <  1:73384/90‑1 (MQ=255)                                                                          aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc       <  1:718171/90‑1 (MQ=255)                                                                            cAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCa     <  2:667641/90‑1 (MQ=255)                                                                              aaTAACGCGTTAAACGGCTTCATGGTCCATATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa   <  2:81386/90‑1 (MQ=255)                                                                               aTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAc  <  1:354423/90‑1 (MQ=255)                                                                                      |                                                                                   CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC  >  NZ_CP009273/3113564‑3113731

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