Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,484,916	A→G	Q181Q (CAA→CAG)	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,484,916	0	A	G	100.0%	176.0 / NA	46	Q181Q (CAA→CAG)	evgS	hybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/24);  total (22/24)

AACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTC  >  NC_000913/2484789‑2485049                                                                                                                                |                                                                                                                                      aaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCaaa                                                                                                                            <  1:1007612/139‑1 (MQ=255)  accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCaaaa                                                                                                                           >  1:582700/1‑139 (MQ=255)     aTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGca                                                                                                                        <  1:327330/139‑1 (MQ=255)       tGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGcaac                                                                                                                      >  1:711369/1‑139 (MQ=255)         aTAATTACCTTTCCGGCGCTGGTAGCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGcaacaa                                                                                                                    >  2:839894/1‑139 (MQ=255)             ttACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTAt                                                                                                                <  1:262660/139‑1 (MQ=255)                     ccGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                                                        <  1:293254/139‑1 (MQ=255)                           cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAAtt                                                                                                  >  2:1143147/1‑139 (MQ=255)                                 accaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                                                            <  1:839894/139‑1 (MQ=255)                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAtt                                                                                     >  1:1036283/1‑139 (MQ=255)                                            cGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                                                 <  2:298496/139‑1 (MQ=255)                                            cGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                                                 <  2:582700/139‑1 (MQ=255)                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtctc                                                                         >  2:71624/1‑139 (MQ=255)                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTg                                                                    >  2:81596/1‑139 (MQ=255)                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                                             <  2:446560/111‑1 (MQ=255)                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                                             >  1:446560/1‑111 (MQ=255)                                                                 cTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTAtc                                                                                                               <  1:69847/88‑1 (MQ=255)                                                                 cTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTAtc                                                                                                               >  2:69847/1‑88 (MQ=255)                                                                         aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCCTTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                                                             <  2:818308/98‑1 (MQ=255)                                                                         aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTGCAAATTTatat                                                                                             >  1:818308/1‑98 (MQ=255)                                                                         aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACtt                                                    >  1:710576/1‑139 (MQ=255)                                                                             aGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAATTTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                                                <  2:715515/139‑1 (MQ=255)                                                                                aaaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                                                >  1:844530/1‑136 (MQ=255)                                                                                aaaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                                                <  2:844530/136‑1 (MQ=255)                                                                                    atTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                                                <  1:528312/132‑1 (MQ=255)                                                                                    atTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATAATTACTTTAtt                                                >  2:528312/1‑132 (MQ=255)                                                                                      tGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                                                 <  2:1062171/97‑1 (MQ=255)                                                                                      tGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                                                 >  1:1062171/1‑97 (MQ=255)                                                                                      tGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCt                                                                     >  2:908156/1‑109 (MQ=255)                                                                                      tGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCt                                                                     <  1:908156/109‑1 (MQ=255)                                                                                       gCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAca                                      <  2:1036283/139‑1 (MQ=255)                                                                                        cTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                                                >  2:576351/1‑128 (MQ=255)                                                                                        cTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                                                <  1:576351/128‑1 (MQ=255)                                                                                         tCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTcc                                                                                                                               >  2:342374/1‑48 (MQ=255)                                                                                         tCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTcc                                                                                                                               <  1:342374/48‑1 (MQ=255)                                                                                             gTAGCAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTa                                >  2:180367/1‑139 (MQ=255)                                                                                                 cAAATTACCCCCCAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcag                            <  1:1143147/139‑1 (MQ=255)                                                                                                        ccccccAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGgtagt                                          <  1:572557/118‑1 (MQ=255)                                                                                                        ccccccAGACGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGgtagt                                          >  2:572557/1‑118 (MQ=255)                                                                                                                aCGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                                                            <  1:477765/60‑1 (MQ=255)                                                                                                                aCGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                                                            >  2:477765/1‑60 (MQ=255)                                                                                                                aCGAGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCg             <  2:710576/139‑1 (MQ=255)                                                                                                                   aGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAg                                                                       <  1:777582/78‑1 (MQ=255)                                                                                                                   aGGTAATTCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAg                                                                       >  2:777582/1‑78 (MQ=255)                                                                                                                         ttCATCAGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTtcac    >  1:334474/1‑139 (MQ=255)                                                                                                                           catcaGTCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTtcactc  <  1:180367/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                      AACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTC  >  NC_000913/2484789‑2485049

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