Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,485,636	Δ1 bp	coding (1263/3594 nt)	evgS →	sensor histidine kinase EvgS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,485,636	0	A	.	100.0%	134.9 / NA	29	coding (1263/3594 nt)	evgS	sensor histidine kinase EvgS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (16/13);  total (16/13)

TCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCAATTAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTTCGTTCGCTTTTCCTCGCGGAGAACCGGAACTTAAGGATATT  >  NC_000913/2485403‑2485867                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        tctcCGATCCCCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATGAAAAAAAGGAATGAAAGTTGCGATGCCAGATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATcc                                                                                                                                                                                                                        <  1:347908/250‑1 (MQ=255)           ccGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAAt                                                                                                                                                                                                             >  1:467363/1‑251 (MQ=255)             gTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATgg                                                                                                                                                                                                           >  1:318681/1‑251 (MQ=255)                                      cacaGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAACGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTg                                                                                                                                                                                  >  2:206863/1‑251 (MQ=255)                                              aGCCCCGGAGGATGGGATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTcac                                                                                                                                                                          >  1:156740/1‑251 (MQ=255)                                                           gggATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCGTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTaaggaa                                                                                                                                                              <  2:467363/250‑1 (MQ=255)                                                                                         aaGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGc                                                                                                                               >  1:558161/1‑251 (MQ=255)                                                                                              cGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAAt                                                                                                                          >  1:40826/1‑251 (MQ=255)                                                                                                 gaAAATAATGTGTTATTTGCTGAAGCCGTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGATGTTGAATTGATAGAGGTCGATAATGCCAGCGTTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTcg                                                                                                                       <  1:578843/251‑1 (MQ=255)                                                                                                                 ttGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATgatcgatc                                                                                                       >  1:583313/1‑251 (MQ=255)                                                                                                                    cTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCAtt                                                                                                    >  1:510355/1‑251 (MQ=255)                                                                                                                                                  ttGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCtt                                                                       <  1:148709/250‑1 (MQ=255)                                                                                                                                                        tGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTg                                                                >  1:43765/1‑251 (MQ=255)                                                                                                                                                            aaaaGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGt                                                            <  2:139956/251‑1 (MQ=255)                                                                                                                                                            aaaaGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCAGTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGt                                                            <  2:294099/251‑1 (MQ=255)                                                                                                                                                                gcgcCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCg                                                        <  2:422405/251‑1 (MQ=255)                                                                                                                                                                   ccTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTAGCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAAt                                                     <  2:193295/251‑1 (MQ=255)                                                                                                                                                                          gTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGAGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGc                                              >  1:302036/1‑251 (MQ=255)                                                                                                                                                                            gAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCt                                             <  1:202070/250‑1 (MQ=255)                                                                                                                                                                             aacaaacaTTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCttt                                           <  2:214304/251‑1 (MQ=255)                                                                                                                                                                                 aacaTTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTttcgtt                                       >  2:19270/1‑251 (MQ=255)                                                                                                                                                                                  acaTTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTttcgttc                                      >  1:368358/1‑251 (MQ=255)                                                                                                                                                                                   caTTAAAAAAATGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTttcgttcg                                     >  2:358082/1‑251 (MQ=255)                                                                                                                                                                                    aTTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTTCGTTCGc                                    >  1:131062/1‑251 (MQ=255)                                                                                                                                                                                            aaGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTTCGTTCGCTTTTCCTc                            <  2:435673/251‑1 (MQ=255)                                                                                                                                                                                              ggAATGACAGTTGCCATTCCATATTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATCATACCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTCTCCTAATGACCTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTTCGTTCGCTTTTCCTcgc                          >  2:386134/1‑251 (MQ=255)                                                                                                                                                                                                      aGTTGCCATTCCATATTATTAGGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTGCGTTCGCGTTTCCTCGCTTAGAAcc                  <  2:558161/251‑1 (MQ=255)                                                                                                                                                                                                             aTTCCATCTTATTATGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTTCGTTCGCTTTTCCTCGCGGAGAACCGGAACtt           >  2:54412/1‑251 (MQ=255)                                                                                                                                                                                                                       attatGAGCTGCATTCGCA‑TAAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTTCGTTCGCTTTTCCTCGCGGAGAACCGGAACTTAAGGAtatt  <  1:582143/250‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        TCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCGGAGGATGGGATATAATACCTGGCGCTATTTATAGTGAAGATCGAGAAAATAATGTTTTATTTGCTGAAGCCTTCATAACAACGCCTTACGTTTTTGTCATGCAAAAAGCGCCTGACAGTGAACAAACATTAAAAAAAGGAATGAAAGTTGCCATTCCATATTATTATGAGCTGCATTCGCAATTAAAAGAGATGTATCCGGAGGTTGAATGGATACAGGTCGATAATGCCAGCGCTGCATTTCACAAGGTTAAGGAAGGTGAACTTGATGCTCTGGTCGCGACACAGCTAAATTCGCGTTACATGATCGATCATTACTATCCTAATGAACTTTATCATTTTCTTATTCCTGGCGTTCCGAATGCATCGCTTTCGTTCGCTTTTCCTCGCGGAGAACCGGAACTTAAGGATATT  >  NC_000913/2485403‑2485867

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