Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	49,884	C→A	P21Q (CCG→CAG)	folA →	type 3 dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	49,884	0	C	A	100.0%	145.9 / NA	40	P21Q (CCG→CAG)	folA	type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (25/15);  total (25/15)

ACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAA  >  NZ_CP009273/49759‑50017                                                                                                                              |                                                                                                                                      aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGGGGAAcctgcct                                                                                                                          >  1:674579/1‑139 (MQ=255)     gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACACCTTGCAGTGGAACCtgactgccc                                                                                                                      >  2:643459/4‑138 (MQ=255)                tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCggg                                                                                                           >  2:275254/1‑139 (MQ=255)                         tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACg                                                                                                   >  1:254453/1‑138 (MQ=255)                         tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACg                                                                                                   >  2:915611/1‑138 (MQ=255)                         tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACg                                                                                                   >  1:390503/1‑138 (MQ=255)                         tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGc                                                                                                  >  2:769608/1‑139 (MQ=255)                         tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGGGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACg                                                                                                   >  1:417788/1‑138 (MQ=255)                           cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                                                <  2:24773/139‑1 (MQ=255)                           cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCa                                                                                                 <  2:663735/138‑1 (MQ=255)                               taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                                                >  2:674405/1‑135 (MQ=255)                               taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                                                <  1:674405/135‑1 (MQ=255)                               taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                                                                            <  1:275254/139‑1 (MQ=255)                                   ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATAGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTta                                                                                         >  2:837122/1‑138 (MQ=255)                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaa                                                                                       >  1:593346/1‑138 (MQ=255)                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaa                                                                                       >  1:479350/1‑138 (MQ=255)                                     cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaa                                                                                       >  2:358770/1‑138 (MQ=255)                                                    tCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATg                                                                       <  1:763645/139‑1 (MQ=255)                                                    tCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATg                                                                       <  2:390503/139‑1 (MQ=255)                                                      gggAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                                     <  1:873457/139‑1 (MQ=255)                                                         aaaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                                                                  <  1:915611/139‑1 (MQ=255)                                                                 aTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                                          <  2:479350/139‑1 (MQ=255)                                                                       aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                    >  1:926644/1‑139 (MQ=255)                                                                       aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                    >  2:144961/1‑139 (MQ=255)                                                                       aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                    >  2:360391/1‑139 (MQ=255)                                                                        gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAGCGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                    >  1:467544/1‑138 (MQ=255)                                                                        gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                    >  1:495281/1‑138 (MQ=255)                                                                               tGCGGCGTTAGCGGTAGATCGCTTTATGGGCATGGAAAACGCCGTCCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                                             <  2:520371/138‑1 (MQ=255)                                                                                gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt                                            <  1:469452/138‑1 (MQ=255)                                                                                   gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGt                                                                              >  2:742145/1‑101 (MQ=255)                                                                                   gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGt                                                                              <  1:742145/101‑1 (MQ=255)                                                                                   gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc                                        >  1:820516/1‑139 (MQ=255)                                                                                   gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc                                        >  1:384850/1‑139 (MQ=255)                                                                                              aGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAAACGGTCGTCCGTTGCCAGGa                             >  1:799724/1‑139 (MQ=255)                                                                                                aTCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                                      >  2:453474/1‑96 (MQ=255)                                                                                                aTCGCGTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                                      <  1:453474/96‑1 (MQ=255)                                                                                                    cgTTATCGGCATGGAAAACGCCATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGGTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                        >  2:796263/1‑138 (MQ=255)                                                                                                                         cATGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCaa  <  1:358770/139‑1 (MQ=255)                                                                                                                           tGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAc                                                                                  <  2:432205/57‑1 (MQ=255)                                                                                                                           tGCAGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAc                                                                                  >  1:432205/1‑57 (MQ=255)                                                                                                                              |                                                                                                                                      ACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAA  >  NZ_CP009273/49759‑50017

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