Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,172,602	A→C	E477A (GAG→GCG)	rpoB →	DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,172,602	0	A	C	100.0%	139.7 / NA	40	E477A (GAG→GCG)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/25);  total (15/25)

CATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCG  >  NZ_CP009273/4172468‑4172732                                                                                                                                       |                                                                                                                                   cATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCgt                                                                                                                                >  1:229730/1‑139 (MQ=255)  aTCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGTCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCgt                                                                                                                                <  1:69167/138‑1 (MQ=255)  aTCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCgtc                                                                                                                               >  1:131760/1‑139 (MQ=255)           cGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTgg                                                                                                                      <  2:417178/139‑1 (MQ=255)           cGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTgg                                                                                                                      <  2:281222/139‑1 (MQ=255)              aaCGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGc                                                                                                                    >  2:412664/1‑138 (MQ=255)               aCGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGa                                                                                                                  >  1:124977/1‑139 (MQ=255)               aCGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGa                                                                                                                  <  2:231577/139‑1 (MQ=255)               aCGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGa                                                                                                                  <  2:131760/139‑1 (MQ=255)                cGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGAt                                                                                                                 <  1:319778/139‑1 (MQ=255)                   tAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTg                                                                                                              >  1:340195/1‑139 (MQ=255)                      aGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGAt                                                                                                           <  1:409304/139‑1 (MQ=255)                             gTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                                                                                                    <  1:233080/139‑1 (MQ=255)                                  tgatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCca                                                                                               <  2:64577/139‑1 (MQ=255)                                  tgatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCca                                                                                               <  1:331070/139‑1 (MQ=255)                                    atatCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcac                                                                                              <  1:240169/138‑1 (MQ=255)                                    atatCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcac                                                                                              <  1:150431/138‑1 (MQ=255)                                    atatCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcac                                                                                              <  2:238499/138‑1 (MQ=255)                                    atatCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcac                                                                                              <  2:116019/138‑1 (MQ=255)                                                 cGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAAc                                                                                >  2:262042/1‑139 (MQ=255)                                                                  gttccgttGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCgc                                                               <  1:110776/139‑1 (MQ=255)                                                                    tccgttGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCgcag                                                             >  1:8725/1‑139 (MQ=255)                                                                    tccgttGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCgcag                                                             >  1:344242/1‑139 (MQ=255)                                                                                                         gCCTGGTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCa                        <  2:388699/139‑1 (MQ=255)                                                                                                              gTACGTGTAGAGCTTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTa                   <  2:117376/139‑1 (MQ=255)                                                                                                              gTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTa                   <  2:437727/139‑1 (MQ=255)                                                                                                              gTACGTGTAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTa                   <  2:124977/139‑1 (MQ=255)                                                                                                                   tgtAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGAc              >  1:218533/1‑139 (MQ=255)                                                                                                                   tgtAGAGCGTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGAc              <  1:351012/139‑1 (MQ=255)                                                                                                                           gTGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGaacaa      >  1:47984/1‑139 (MQ=255)                                                                                                                            tGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCgcag                                                             >  2:254884/1‑83 (MQ=255)                                                                                                                            tGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCgcag                                                             <  1:254884/83‑1 (MQ=255)                                                                                                                            tGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTa                   >  1:339932/1‑125 (MQ=255)                                                                                                                            tGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGaacaac     >  2:281588/1‑139 (MQ=255)                                                                                                                            tGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGaacaac     >  1:79294/1‑139 (MQ=255)                                                                                                                            tGCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCGCAGTTTa                   <  2:339932/125‑1 (MQ=255)                                                                                                                             gCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCgc                                                               >  2:192193/1‑80 (MQ=255)                                                                                                                             gCGGTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCgc                                                               <  1:192193/80‑1 (MQ=255)                                                                                                                               ggTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCg  <  1:164642/139‑1 (MQ=255)                                                                                                                               ggTGAAAGCGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCg  <  1:14636/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   CATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCG  >  NZ_CP009273/4172468‑4172732

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