Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,854,126	T→C	M121V (ATG→GTG)	ydjK ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,854,126	0	T	C	97.7%	129.7 / ‑4.0	43	M121V (ATG→GTG)	ydjK	MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (17/25);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

CGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGA  >  NZ_CP009273/1854046‑1854209                                                                                 |                                                                                    cGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACg                                                                            >  2:1154041/1‑90 (MQ=255) cGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACg                                                                            >  1:1154263/1‑90 (MQ=255) cGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACg                                                                            >  1:1060717/1‑90 (MQ=255)   ttCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggc                                                                          <  1:335276/90‑1 (MQ=255)     ccACGATTTCTCCCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcag                                                                        <  2:731794/90‑1 (MQ=255)     ccACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcag                                                                        <  1:1329687/90‑1 (MQ=255)     ccACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcag                                                                        <  2:711656/90‑1 (MQ=255)     ccACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcag                                                                        <  2:735386/90‑1 (MQ=255)       aCGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcaggc                                                                      <  2:1123344/90‑1 (MQ=255)       aCGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACggcaggc                                                                      >  1:786187/1‑90 (MQ=255)           tttCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATa                                                                  >  2:1371389/1‑90 (MQ=255)             tCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAg                                                                >  2:19401/1‑90 (MQ=255)               tACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGa                                                              >  2:840286/1‑90 (MQ=255)               tACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGa                                                              >  2:354424/1‑90 (MQ=255)               tACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGa                                                              >  1:973818/1‑90 (MQ=255)                aCCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGaa                                                             <  2:687065/90‑1 (MQ=255)                        tatatTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCatat                                                     >  1:521563/1‑90 (MQ=255)                             tCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCCGCCCAACTCCCACCCCAAAACGGCAGGCAATAAGGAAATCCATATTCGGt                                                <  1:1180254/90‑1 (MQ=255)                              cGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTg                                               >  1:353712/1‑90 (MQ=255)                              cGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTg                                               >  1:993539/1‑90 (MQ=255)                                       cAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCAc                                      <  2:1100160/90‑1 (MQ=255)                                                 cGTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCAt                            <  1:1086423/90‑1 (MQ=255)                                                  gTAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATc                           <  2:563838/90‑1 (MQ=255)                                                   tAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCg                          <  2:1541286/90‑1 (MQ=255)                                                   tAACCAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCg                          <  1:1335813/90‑1 (MQ=255)                                                      ccAGTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAgg                       <  2:1066664/90‑1 (MQ=255)                                                        aGTAAAGCTCCCGGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCa                     >  1:735119/1‑90 (MQ=255)                                                        aGTAAAGCTCCCAGCCCAACTCCCCCCCCAAAACGGAAGGCAATAACGAAATCCATGTTCGGTGAGAAAGCACCGCCAACCATCGAGGCa                     <  1:876781/90‑1 (MQ=255)                                                         gTAAAGCTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCaa                    >  2:541366/1‑90 (MQ=255)                                                               cTCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGaa              >  1:658993/1‑90 (MQ=255)                                                                tCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  2:893107/90‑1 (MQ=255)                                                                tCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  1:810723/90‑1 (MQ=255)                                                                tCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  2:632179/90‑1 (MQ=255)                                                                tCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  2:483258/90‑1 (MQ=255)                                                                tCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  2:353712/90‑1 (MQ=255)                                                                tCCCAGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  1:963271/90‑1 (MQ=255)                                                                tCCCAGCCCAACTCCCACCACAAAACGCCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAAt             <  2:1191593/90‑1 (MQ=255)                                                                  ccAGACCAAATCCCACCAAAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATgg           <  1:523403/90‑1 (MQ=255)                                                                    aGCCCAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGcc         <  2:1546958/90‑1 (MQ=255)                                                                      cccAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTc       <  1:1557470/90‑1 (MQ=255)                                                                        cAACTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCAt     <  2:472513/90‑1 (MQ=255)                                                                          aCTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATag   >  2:1066154/1‑90 (MQ=255)                                                                           cTCCCACCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATaga  >  1:923056/1‑90 (MQ=255)                                                                                 |                                                                                    CGTTCCACGATTTCTACCGGGCATATATTCGGTGAAACCAGCAAACAGCGTAACCAGTAAAGCTCCCAGCCCAACTCCCATCACAAAACGGCAGGCAATAAGGAAATCCATATTCGGTGAGAAAGCACCGACAACCATCGAGGCAATATGAATGGCCTCATAGA  >  NZ_CP009273/1854046‑1854209

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