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%	139.2 / NA	38	T125A (ACC→GCC)	yfeO	ion channel protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/23);  total (15/23)

TCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCT  >  NZ_CP009273/2503393‑2503562                                                                                         |                                                                                  tCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGTCCGGAACATCCGATCATGGcc                                                                                 >  2:227907/1‑91 (MQ=255)   cGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGt                                                                               <  1:215743/91‑1 (MQ=255)    gCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTc                                                                              <  2:440435/91‑1 (MQ=255)         ctGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAatat                                                                         >  2:408049/1‑91 (MQ=255)         ctGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAatat                                                                         >  1:50334/1‑91 (MQ=255)                aCCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCtt                                                                  <  2:389970/91‑1 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGt                                                             <  1:355937/86‑1 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGt                                                             >  2:355937/1‑86 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCa                                                        >  2:626505/1‑91 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCa                                                        >  2:594408/1‑91 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCa                                                        >  2:437471/1‑91 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCa                                                        >  1:383744/1‑91 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCa                                                        >  1:361279/1‑91 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGACCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCa                                                        >  2:171236/1‑91 (MQ=255)                          aTCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGCGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCa                                                        >  1:679001/1‑91 (MQ=255)                                  aTTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCt                                                <  2:601089/91‑1 (MQ=255)                                  aTTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCt                                                <  1:95901/91‑1 (MQ=255)                                  aTTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCt                                                <  1:80454/91‑1 (MQ=255)                                  aTTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCt                                                <  1:22193/91‑1 (MQ=255)                                                  gCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACcgcgcg                                <  2:383744/91‑1 (MQ=255)                                                    tGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATGTCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTc                              <  1:626505/91‑1 (MQ=255)                                                      gcggcgTCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCaa                            <  2:453607/91‑1 (MQ=255)                                                            tCAGCCTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAAt                      <  2:282663/91‑1 (MQ=255)                                                                ccTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAg                  <  1:40901/91‑1 (MQ=255)                                                                 cTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGt                 <  1:227907/91‑1 (MQ=255)                                                                 cTGGGGCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGt                 <  1:408049/91‑1 (MQ=255)                                                                   ggggCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAAc                           <  2:570435/79‑1 (MQ=255)                                                                   ggggCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAAc                           >  1:570435/1‑79 (MQ=255)                                                                   ggggCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTgg               >  1:75461/1‑91 (MQ=255)                                                                     ggCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGAc             <  1:265103/91‑1 (MQ=255)                                                                     ggCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGAc             <  2:50334/91‑1 (MQ=255)                                                                     ggCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGAc             <  1:732650/91‑1 (MQ=255)                                                                     ggCCGGAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGAc             <  1:696923/91‑1 (MQ=255)                                                                         ggAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTAtt         >  1:868133/1‑91 (MQ=255)                                                                         ggAACATCCGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTAtt         >  1:553885/1‑91 (MQ=255)                                                                                ccGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCct  <  2:217371/91‑1 (MQ=255)                                                                                ccGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCct  <  2:455475/91‑1 (MQ=255)                                                                                ccGATCATGGCCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCct  <  1:135251/91‑1 (MQ=255)                                                                                         |                                                                                  TCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCT  >  NZ_CP009273/2503393‑2503562

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