Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	49,791	G→A	intergenic (+160/‑32)	kefC → / → folA	glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	49,791	0	G	A	100.0%	101.6 / NA	27	intergenic (+160/‑32)	kefC/folA	glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/11);  total (16/11)

AGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAAC  >  NZ_CP009273/49707‑49876                                                                                     |                                                                                      agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                                  >  2:403220/1‑90 (MQ=255) agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                                  >  1:825395/1‑90 (MQ=255) agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                                  >  1:588637/1‑90 (MQ=255) agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                                  >  1:492880/1‑90 (MQ=255)   agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGAc                                                                                >  1:41580/1‑90 (MQ=255)      gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAAt                                                                             >  2:383634/1‑90 (MQ=255)        ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTT‑CGTATTATTGCATCAAtttt                                                                          >  1:757902/1‑90 (MQ=255)                     gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGaa                                                              <  1:617153/90‑1 (MQ=255)                       atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAAt                                                            <  2:116507/90‑1 (MQ=255)                                 ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                                  >  1:565952/1‑90 (MQ=255)                                 ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                                  <  2:591426/90‑1 (MQ=255)                                 ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                                  <  1:294696/90‑1 (MQ=255)                                                  cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                                 <  1:368395/90‑1 (MQ=255)                                                        aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                           >  2:434741/1‑90 (MQ=255)                                                              cgacgGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATcg                     <  1:383634/90‑1 (MQ=255)                                                                acgGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATcgcg                   >  2:247278/1‑90 (MQ=255)                                                                   gTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTa                >  1:351148/1‑90 (MQ=255)                                                                            ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcggcg                                   <  1:401430/62‑1 (MQ=255)                                                                            ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcggcg                                   >  2:401430/1‑62 (MQ=255)                                                                            ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATgg       >  1:482013/1‑90 (MQ=255)                                                                             tACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                                 <  2:802546/63‑1 (MQ=255)                                                                             tACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                                 >  1:180036/1‑63 (MQ=255)                                                                             tACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                                 <  2:180036/63‑1 (MQ=255)                                                                             tACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                                 >  1:802546/1‑63 (MQ=255)                                                                              aCGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaa     <  2:588637/90‑1 (MQ=255)                                                                               cGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaa    <  2:757902/90‑1 (MQ=255)                                                                                 tataATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAAc  >  2:113149/1‑90 (MQ=255)                                                                                     |                                                                                      AGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAAC  >  NZ_CP009273/49707‑49876

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