Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	960,653	T→C	M293T (ATG→ACG)	ycaI →	ComEC family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	960,653	0	T	C	100.0%	77.5 / NA	23	M293T (ATG→ACG)	ycaI	ComEC family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/12);  total (11/12)

CCTGGGCGCTGGATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAATGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGCTGCGGTCGCGGCATTG  >  NZ_CP009273/960526‑960786                                                                                                                                |                                                                                                                                      ccTGGGCGCTGGATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCttaagtt                                                                                                                              <  1:4520/137‑1 (MQ=255)      gcgcTGGATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTggg                                                                                                                       <  1:201838/139‑1 (MQ=255)      gcgcTGGATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTggg                                                                                                                       <  2:9735/139‑1 (MQ=255)       cgcTGGATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGc                                                                                                                      <  2:5725/139‑1 (MQ=255)         cTGGATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGa                                                                                                                    <  1:322791/139‑1 (MQ=255)            gATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGAc                                                                                                                   <  1:291668/137‑1 (MQ=255)                 aCTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACagtggag                                                                                                            >  1:392694/1‑139 (MQ=255)                 aCTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACagtggag                                                                                                            >  2:407873/1‑139 (MQ=255)                               ttAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTAt                                                                                              >  2:145566/1‑139 (MQ=255)                                                     ctgctTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAAt                                                                        >  2:140108/1‑139 (MQ=255)                                                        ctTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAAtt                                                                       <  2:392694/137‑1 (MQ=255)                                                                     gTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATc                                                          <  1:145566/137‑1 (MQ=255)                                                                        gACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGt                                                      <  2:30034/138‑1 (MQ=255)                                                                                        cctGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCa                                       >  1:309521/1‑137 (MQ=255)                                                                                                 cGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATcc                                                         <  1:303429/110‑1 (MQ=255)                                                                                                 cGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATcc                                                         >  2:303429/1‑110 (MQ=255)                                                                                                 cGTACCATGGTGGCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGc                            >  2:262027/1‑139 (MQ=255)                                                                                                           tggCGCTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGGGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGc                  >  1:35168/1‑139 (MQ=255)                                                                                                                cTTGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGCTGCgg             >  1:237485/1‑139 (MQ=255)                                                                                                                  tGCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGCTGCGGTc           >  1:93590/1‑139 (MQ=255)                                                                                                                   gCTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGCTGCGGTcg          >  2:358314/1‑139 (MQ=255)                                                                                                                    cTACGTGGGGAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGCTGCGGTcgc         <  2:237485/139‑1 (MQ=255)                                                                                                                           gggAACGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGCTGCGGTCGCGGCATTg  <  2:384353/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                      CCTGGGCGCTGGATCCACTGGCAAATACCATTAATTGGCGGAATCTGCTGTGCTGCTTTTTATGCCTGGTTGACGGGAATGCAACCTCCTGCATTGCGTACCATGGTGGCGCTTGCTACGTGGGGAATGCTTAAGTTAAGTGGGCGACAGTGGAGTGGCTGGGATGTATGGATATGTTGTCTGGCGGCAATTTTGCTGATGGATCCTGTTGCCATTCTCTCGCAAAGTTTATGGCTCTCTGCCGCTGCGGTCGCGGCATTG  >  NZ_CP009273/960526‑960786

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