Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,503,481	A→G	T125A (ACC→GCC)	yfeO →	ion channel protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,503,481	0	A	G	100.0%	114.8 / NA	32	T125A (ACC→GCC)	yfeO	ion channel protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

TGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTGCGGCGGCG  >  NZ_CP009273/2503348‑2503609                                                                                                                                      |                                                                                                                                 tGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGt                                                                                                                              >  1:170637/1‑138 (MQ=255)           gACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGccc                                                                                                                   <  2:553723/139‑1 (MQ=255)            aCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCt                                                                                                                  >  1:397028/1‑139 (MQ=255)                      tgAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGc                                                                                                         >  2:82953/1‑138 (MQ=255)                        aaCCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCAt                                                                                                      >  1:44851/1‑139 (MQ=255)                                aTCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGgcgc                                                                                                <  2:550435/137‑1 (MQ=255)                                aTCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTc                                                                                              <  2:161706/139‑1 (MQ=255)                                                cGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGt                                                                              <  2:44851/139‑1 (MQ=255)                                                cGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGt                                                                              <  2:303918/139‑1 (MQ=255)                                                cGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGt                                                                              <  2:170637/139‑1 (MQ=255)                                                  cccTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTc                                                                             >  1:211532/1‑138 (MQ=255)                                                   ccTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCaa                                                                           <  1:499076/139‑1 (MQ=255)                                                      ctGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAAcc                                                                         >  2:484365/1‑138 (MQ=255)                                                       tGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAatat                                                                                                                        <  1:152013/90‑1 (MQ=255)                                                       tGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAatat                                                                                                                        >  2:152013/1‑90 (MQ=255)                                                       tGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCg                                                                        >  1:297835/1‑138 (MQ=255)                                                       tGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCg                                                                        >  1:561117/1‑138 (MQ=255)                                                         cgcTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGa                                                                       <  1:496768/137‑1 (MQ=255)                                                                 ggACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                                             >  1:441913/1‑139 (MQ=255)                                                                  gACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                                             <  2:211532/138‑1 (MQ=255)                                                                  gACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                                             <  2:216985/138‑1 (MQ=255)                                                                       aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                                             <  2:17749/133‑1 (MQ=255)                                                                       aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                                             >  1:17749/1‑133 (MQ=255)                                                                        tCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATg                                                                    >  1:82070/1‑125 (MQ=255)                                                                        tCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATg                                                                    <  2:82070/125‑1 (MQ=255)                                                                                                         tCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTgg                     <  2:441913/139‑1 (MQ=255)                                                                                                                      ggAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTgcg        >  1:467947/1‑139 (MQ=255)                                                                                                                       gAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCct                                                 <  1:372791/97‑1 (MQ=255)                                                                                                                       gAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCct                                                 >  2:372791/1‑97 (MQ=255)                                                                                                                       gAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTgcgg       >  1:331200/1‑139 (MQ=255)                                                                                                                          cATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTgcggcgg    <  1:82953/139‑1 (MQ=255)                                                                                                                            tCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTgcggcggtg  >  2:429380/1‑137 (MQ=255)                                                                                                                                      |                                                                                                                                 TGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGGCACTCCTGTTGCGGCGGCG  >  NZ_CP009273/2503348‑2503609

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