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	93.9%	108.7 / ‑1.4	33	K330N (AAA→AAC)	glcA	glycolate permease GlcA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (13/18);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTA  >  NZ_CP009273/3113564‑3113733                                                                                      |                                                                                     ccGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaa                                                                                  >  2:55057/1‑90 (MQ=255)      ttCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGc                                                                             >  2:468489/1‑90 (MQ=255)      ttCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGc                                                                             <  2:205331/90‑1 (MQ=255)      ttCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGc                                                                             <  1:400890/90‑1 (MQ=255)          acacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCa                                                                         <  1:614958/90‑1 (MQ=255)           cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  1:763671/90‑1 (MQ=255)           cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                        <  2:43992/90‑1 (MQ=255)            acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                       <  2:71122/90‑1 (MQ=255)                 ttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                  <  1:286488/90‑1 (MQ=255)                  tgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACgg                                                                              <  1:137334/77‑1 (MQ=255)                  tgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACgg                                                                              >  2:137334/1‑77 (MQ=255)                    aTGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGAt                                                               >  1:69534/1‑90 (MQ=255)                     tGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATg                                                              <  2:650165/90‑1 (MQ=255)                      gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                             <  2:427965/90‑1 (MQ=255)                           tGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGAcc                                                        <  2:554786/90‑1 (MQ=255)                                  tCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCg                                                 >  1:267626/1‑90 (MQ=255)                                   cTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                     <  1:133398/85‑1 (MQ=255)                                   cTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                     >  2:133398/1‑85 (MQ=255)                                                    gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTATACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTg                               >  2:433766/1‑90 (MQ=255)                                                    gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTg                               >  2:45806/1‑90 (MQ=255)                                                    gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTg                               >  2:539776/1‑90 (MQ=255)                                                      gAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGAc                             >  1:505556/1‑90 (MQ=255)                                                       aaTAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGAcc                            <  1:449210/90‑1 (MQ=255)                                                           aaaCGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGc                        >  2:88065/1‑90 (MQ=255)                                                              cgcgCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCg                     >  2:59195/1‑90 (MQ=255)                                                                cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                   >  2:291834/1‑90 (MQ=255)                                                                     cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              <  1:88065/90‑1 (MQ=255)                                                                          aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc         <  1:37976/90‑1 (MQ=255)                                                                              aaTAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa     <  1:468489/90‑1 (MQ=255)                                                                              aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              >  2:153247/1‑81 (MQ=255)                                                                              aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt              <  1:153247/81‑1 (MQ=255)                                                                               aTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAc    <  1:115195/90‑1 (MQ=255)                                                                                 aaCGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACta  <  1:486844/90‑1 (MQ=255)                                                                                      |                                                                                     CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTA  >  NZ_CP009273/3113564‑3113733

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