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%	160.3 / NA	42	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 (26/16);  total (26/16)

GAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGA  >  NZ_CP009273/49706‑49872                                                                                      |                                                                                  gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGc                                                                               >  2:344655/1‑90 (MQ=255)  agagCGTGGACTCGCCGGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                              >  2:59468/1‑90 (MQ=255)  agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                              >  1:18953/1‑90 (MQ=255)  agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                                              >  1:742584/1‑90 (MQ=255)   gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGa                                                                             >  1:165031/1‑90 (MQ=255)     gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACa                                                                           >  1:516719/1‑90 (MQ=255)        tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAAtt                                                                        >  1:802313/1‑90 (MQ=255)               gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTa                                                                 <  2:742584/90‑1 (MQ=255)                     agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGa                                                           <  2:795910/90‑1 (MQ=255)                        atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAAt                                                        <  1:344655/90‑1 (MQ=255)                        atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAAt                                                        <  1:614607/90‑1 (MQ=255)                             aaaTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACGATTTTTTTTATCGGGAAATCTCaa                                                   >  1:388723/1‑90 (MQ=255)                             aaaTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCaa                                                   >  1:724363/1‑90 (MQ=255)                                ttttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGTTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGa                                                <  2:248343/90‑1 (MQ=255)                                ttttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGa                                                >  1:143261/1‑90 (MQ=255)                                  ttCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                              <  2:862834/90‑1 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCa                                             <  2:654011/90‑1 (MQ=255)                                         aCATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGa                                       >  1:758583/1‑90 (MQ=255)                                             catcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgc                                   >  1:161063/1‑90 (MQ=255)                                               tcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcgg                                 <  2:699437/90‑1 (MQ=255)                                               tcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcgg                                 <  2:516719/90‑1 (MQ=255)                                                   cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                             <  1:5663/90‑1 (MQ=255)                                                   cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                             >  2:361341/1‑90 (MQ=255)                                                   cGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                             >  1:807389/1‑90 (MQ=255)                                                       ccAGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCg                         <  2:777714/90‑1 (MQ=255)                                                         aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGAtt                                     >  2:158235/1‑76 (MQ=255)                                                         aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGAtt                                     <  1:158235/76‑1 (MQ=255)                                                         aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCggc                       >  1:792133/1‑89 (MQ=255)                                                         aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt                       >  2:360749/1‑90 (MQ=255)                                                          gTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTa                      <  2:689775/90‑1 (MQ=255)                                                                     tttACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt           <  2:397187/90‑1 (MQ=255)                                                                      ttACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATc          <  1:546917/90‑1 (MQ=255)                                                                           cTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCAt     <  1:274965/90‑1 (MQ=255)                                                                            tttACGTATAATGGCGACAATTTTTTTTATCGGTAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATg    >  1:261638/1‑90 (MQ=255)                                                                            tttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATg    >  2:209524/1‑90 (MQ=255)                                                                            tttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATg    >  1:503063/1‑90 (MQ=255)                                                                            tttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATg    >  1:583036/1‑90 (MQ=255)                                                                             ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATgg   >  1:129708/1‑90 (MQ=255)                                                                             ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATgg   >  2:114973/1‑90 (MQ=255)                                                                             ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATgg   >  2:728100/1‑90 (MQ=255)                                                                             ttACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATgg   >  2:571251/1‑90 (MQ=255)                                                                              tACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa  >  1:845619/1‑90 (MQ=255)                                                                                      |                                                                                  GAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGA  >  NZ_CP009273/49706‑49872

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: ‑