Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,485,315	(C)6→7	coding (942/3594 nt)	evgS →	hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,485,309	1	.	C	100.0%	92.7 / NA	28	S312S (TCG→TCC)	evgS	hybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (12/16);  total (12/16)

TTTAACAAATGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCG‑CCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCG  >  NC_000913/2485159‑2485454                                                                                                                                                        |                                                                                                                                                  tttAACAAATGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCg                                                                                                                                                    >  1:77917/1‑151 (MQ=255)       aaaTGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCG‑cccccc                                                                                                                                             >  2:197296/1‑151 (MQ=255)        aaTGAAGTTCGCTATGAAGTATCCCAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCG‑ccccccc                                                                                                                                            >  1:605562/1‑150 (MQ=255)         aTGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCG‑ccccccct                                                                                                                                           >  2:210100/1‑149 (MQ=255)           gAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCG‑ccccccctat                                                                                                                                         <  2:612195/151‑5 (MQ=255)            aaGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCG‑ccccccctatt                                                                                                                                        >  1:663945/1‑146 (MQ=255)                     aTGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACg                                                                                                                                 <  1:169129/150‑1 (MQ=255)                                  aaaTTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTc                                                                                                                   <  1:453350/151‑1 (MQ=255)                                     ttGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGt                                                                                                                >  1:359954/1‑151 (MQ=255)                                       ggCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTc                                                                                                              >  1:443240/1‑151 (MQ=255)                                        gCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGtt                                                                                                               <  1:521553/149‑1 (MQ=255)                                            gATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCgggg                                                                                                          <  1:79461/150‑1 (MQ=255)                                            gATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGc                                                                                                         <  2:163338/151‑1 (MQ=255)                                                          ggCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACa                                                                                            <  1:101181/150‑1 (MQ=255)                                                          ggCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACa                                                                                            <  2:595547/150‑1 (MQ=255)                                                            ccTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTc                                                                                         >  1:587080/1‑151 (MQ=255)                                                                   gAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAat                                                                                    <  2:230042/149‑1 (MQ=255)                                                                    aaCAAACCATTAGAACTCACTGAACATTAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCGGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAat                                                                                    <  1:432990/148‑1 (MQ=255)                                                                      cAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAAtattat                                                                               >  1:316572/1‑151 (MQ=255)                                                                             ttAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTg                                                                        >  1:414697/1‑151 (MQ=255)                                                                                gAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCaa                                                                     <  2:319640/151‑1 (MQ=255)                                                                                gAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCaa                                                                     <  2:449903/151‑1 (MQ=255)                                                                                   cTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAAc                                                                   >  2:609169/1‑150 (MQ=255)                                                                                     cACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAg                                                                 <  1:38958/150‑1 (MQ=255)                                                                                           aCATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTaa                                                           >  2:268670/1‑150 (MQ=255)                                                                                               gAAAAACAGTGGATTAAGCAGCAGCCCAATTTAAAGTTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGTCTCGGTTCGGGGCGTTATATGGGACATTCTTAATATTATTACCTTGCAAACAGGGTTAAAttt                                                       <  1:439605/150‑1 (MQ=255)                                                                                                         ggATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATc                                             >  2:507367/1‑150 (MQ=255)                                                                                                         ggATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCa                                            <  2:49512/151‑1 (MQ=255)                                                                                                                 cagcaTCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTc                                     >  2:508606/1‑150 (MQ=255)                                                                                                                          aaTTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATc                           >  2:455835/1‑151 (MQ=255)                                                                                                                                             aTCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCtt         >  2:563907/1‑150 (MQ=255)                                                                                                                                             aTCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCtt         >  2:565480/1‑150 (MQ=255)                                                                                                                                             aTCCTTACTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTa        <  1:150318/151‑1 (MQ=255)                                                                                                                                                   aCTCGCCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCg  <  2:476119/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  TTTAACAAATGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCG‑CCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGACATTCTTAATATTATTACCTTGCAAACAGGTTTAAATTTTTCTCCGATCACCGTTTCACACAATATCCATGCTGGAACACAGCTTAGCCCCG  >  NC_000913/2485159‑2485454

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