Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	49,898	G→A	A26T (GCC→ACC)	folA →	type 3 dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	49,898	0	G	A	94.7%	123.2 / ‑1.3	39	A26T (GCC→ACC)	folA	type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (16/20);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

GTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGA  >  NZ_CP009273/49763‑50031                                                                                                                                        |                                                                                                                                      gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCtgcctacc                                                                                                                                     >  1:201635/4‑138 (MQ=255) gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctacc                                                                                                                                     >  2:741565/4‑138 (MQ=255)       gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCg                                                                                                                              >  2:765337/1‑139 (MQ=255)                     tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACg                                                                                                                 >  2:740626/1‑138 (MQ=255)                       cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCa                                                                                                               <  2:844777/138‑1 (MQ=255)                       cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCa                                                                                                               <  1:672893/138‑1 (MQ=255)                       cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCa                                                                                                               <  2:115149/138‑1 (MQ=255)                       cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCa                                                                                                               <  2:832413/138‑1 (MQ=255)                       cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCa                                                                                                               <  2:417170/138‑1 (MQ=255)                               ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACAACTta                                                                                                       >  2:739848/1‑138 (MQ=255)                                  gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTtaaata                                                                                                   <  2:201635/139‑1 (MQ=255)                                      aTTTTTTCTATCGGGAAAACTCAATGATCAGTCTGATTTCGGCGTTAGCGGTAGATCGCGTTATCGGCATGAAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAcc                                                                                               <  1:357349/139‑1 (MQ=255)                                         ttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCg                                                                                             <  1:640464/138‑1 (MQ=255)                                                      aaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                                                                               <  1:740626/139‑1 (MQ=255)                                                            aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                                                             <  1:653492/135‑1 (MQ=255)                                                             aTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                                                        <  2:611470/139‑1 (MQ=255)                                                                    gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCCTGTCGTGGAACCTGCCTCCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                                  >  2:711251/1‑138 (MQ=255)                                                                    gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                                  >  2:168843/1‑138 (MQ=255)                                                                     tCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGttt                                                                                                                      >  2:743738/1‑85 (MQ=255)                                                                     tCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGATGTGGAACCTGCCTACCGATCTCGCCGGGttt                                                                                                                      <  1:743738/85‑1 (MQ=255)                                                                            gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt                                                          <  2:846792/138‑1 (MQ=255)                                                                               gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAccc                                                                                              <  1:324755/99‑1 (MQ=255)                                                                               gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAccc                                                                                              >  2:324755/1‑99 (MQ=255)                                                                               gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGGTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAAACGggcgtc                                                      >  1:796727/1‑139 (MQ=255)                                                                               gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc                                                      >  1:208027/1‑139 (MQ=255)                                                                                           gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGa                                           <  1:765337/138‑1 (MQ=255)                                                                                                cgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      >  1:371766/1‑138 (MQ=255)                                                                                                 gTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaaa                                    <  1:739848/139‑1 (MQ=255)                                                                                                     tCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAata                                 <  2:371766/138‑1 (MQ=255)                                                                                                                   gCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGt                   >  1:727817/1‑138 (MQ=255)                                                                                                                       tGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAAcc              >  1:145023/1‑139 (MQ=255)                                                                                                                       tGCTGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAAcc              >  1:6025/1‑139 (MQ=255)                                                                                                                        gctgTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCg             >  1:109796/4‑139 (MQ=255)                                                                                                                        gCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCaa                                                                                                              >  2:104962/1‑42 (MQ=255)                                                                                                                        gCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCaa                                                                                                              <  1:104962/42‑1 (MQ=255)                                                                                                                           gTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGt          <  1:665682/139‑1 (MQ=255)                                                                                                                            tGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCaa                <  2:154933/132‑1 (MQ=255)                                                                                                                            tGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCaa                >  1:154933/1‑132 (MQ=255)                                                                                                                              gAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTAc        <  1:264522/138‑1 (MQ=255)                                                                                                                              gAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACg       >  1:306576/1‑139 (MQ=255)                                                                                                                                aCCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGa     <  2:145023/139‑1 (MQ=255)                                                                                                                                  ctgcctacCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACg   <  1:361364/139‑1 (MQ=255)                                                                                                                                   tgcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCAAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACgc  >  2:684305/1‑138 (MQ=255)                                                                                                                                        |                                                                                                                                      GTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGA  >  NZ_CP009273/49763‑50031

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