Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,848,051	T→G	F257C (TTT→TGT)	emrD →	multidrug efflux MFS transporter EmrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,848,051	0	T	G	100.0%	62.8 / NA	20	F257C (TTT→TGT)	emrD	multidrug efflux MFS transporter EmrD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/6);  total (14/6)

CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCG  >  NZ_CP009273/3847965‑3848112                                                                                       |                                                              cTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGtgttt                                                            >  1:206256/1‑90 (MQ=255)  tGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGtgtttt                                                           <  1:78784/90‑1 (MQ=255)   gCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTg                                                          >  1:244922/1‑90 (MQ=255)   gCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTg                                                          >  2:371544/1‑90 (MQ=255)   gCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTg                                                          >  2:203429/1‑90 (MQ=255)                        cgGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTc                                     <  1:58752/90‑1 (MQ=255)                               aGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATaa                              <  2:535806/90‑1 (MQ=255)                                   cTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGc                          <  1:327360/90‑1 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  2:235839/1‑90 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  2:85448/1‑90 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  2:604498/1‑90 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  2:567467/1‑90 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  2:528240/1‑90 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  1:304846/1‑90 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  1:674697/1‑90 (MQ=255)                                    tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                         >  1:43649/1‑90 (MQ=255)                                                cGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTaa             <  1:351817/90‑1 (MQ=255)                                                    cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg         >  2:389430/1‑90 (MQ=255)                                                    cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg         >  2:600862/1‑90 (MQ=255)                                                        aTTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAg     <  1:567467/90‑1 (MQ=255)                                                          tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc   <  1:389430/90‑1 (MQ=255)                                                           ttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCg  >  1:317147/1‑90 (MQ=255)                                                                                       |                                                              CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCG  >  NZ_CP009273/3847965‑3848112

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