Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,120	G→A	*475* (TAG→TAA)	creC →	sensory histidine kinase in two‑component regulatory system with CreB or PhoB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,120	0	G	A	100.0%	201.6 / NA	62	*475* (TAG→TAA)	creC	sensory histidine kinase in two‑component regulatory system with CreB or PhoB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/34);  total (28/34)

AAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGA  >  NC_000913/4637994‑4638241                                                                                                                               |                                                                                                                          aaGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAttctt                                                                                                               <  2:468234/139‑1 (MQ=255)       cGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                         <  2:1558727/139‑1 (MQ=255)       cGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                         <  1:1590475/139‑1 (MQ=255)       cGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                         <  1:587040/139‑1 (MQ=255)           cTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGt                                                                                                     <  1:1428192/139‑1 (MQ=255)                       ttCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCctg                                                                                         <  2:860786/139‑1 (MQ=255)                                  ggTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTCAACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGGCATTGCa                                                                              <  1:1194595/139‑1 (MQ=255)                                    tCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCaaa                                                                            >  2:1256557/1‑139 (MQ=255)                                      gCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCctgc                                                                                        >  2:581568/1‑125 (MQ=255)                                      gCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCctgc                                                                                        <  1:581568/125‑1 (MQ=255)                                          gtttgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGagaa                                                                      >  2:417549/1‑139 (MQ=255)                                           tttgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAttctt                                                                                                               <  2:1630674/97‑1 (MQ=255)                                           tttgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAttctt                                                                                                               >  1:1630674/1‑97 (MQ=255)                                            ttgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGagaaga                                                                    <  1:525226/139‑1 (MQ=255)                                               tttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTa                                                                 <  1:417549/139‑1 (MQ=255)                                                    cGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCtgctgc                                                                                     >  2:1073050/1‑114 (MQ=255)                                                    cGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCtgctgc                                                                                     <  1:1073050/114‑1 (MQ=255)                                                     ggCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGa                                                           >  1:551145/1‑139 (MQ=255)                                                       cGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                              >  1:996170/1‑118 (MQ=255)                                                       cGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                              <  2:996170/118‑1 (MQ=255)                                                       cGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGaaa                                                         >  2:209312/1‑139 (MQ=255)                                                        gAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                        <  1:870152/139‑1 (MQ=255)                                                             cACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccccc                                                   <  2:551145/139‑1 (MQ=255)                                                             cACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccccc                                                   <  1:1228630/139‑1 (MQ=255)                                                                 cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAtt                                                                                                                  <  2:1666520/72‑1 (MQ=255)                                                                 cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAtt                                                                                                                  >  1:1666520/1‑72 (MQ=255)                                                                 cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                         <  1:969851/81‑1 (MQ=255)                                                                 cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                         >  2:969851/1‑81 (MQ=255)                                                                 cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGt                                               >  1:1632465/1‑139 (MQ=255)                                                                 cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGt                                               >  2:587529/1‑139 (MQ=255)                                                                   gcgcAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAt                                                                                                                   >  1:1315984/1‑69 (MQ=255)                                                                   gcgcAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAt                                                                                                                   <  2:1315984/69‑1 (MQ=255)                                                                    cgcAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                              >  1:1681700/1‑105 (MQ=255)                                                                    cgcAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                              <  2:1681700/105‑1 (MQ=255)                                                                      cAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                              >  1:788688/1‑103 (MQ=255)                                                                      cAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                              <  2:788688/103‑1 (MQ=255)                                                                      cAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTgg                                          >  1:291808/1‑139 (MQ=255)                                                                          gTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                        <  2:1284963/121‑1 (MQ=255)                                                                          gTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                        >  1:1284963/1‑121 (MQ=255)                                                                                        cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGga                                                                         >  2:858101/1‑90 (MQ=255)                                                                                        cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGga                                                                         <  1:858101/90‑1 (MQ=255)                                                                                        cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGttt                        >  1:151273/1‑139 (MQ=255)                                                                                          tGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTa                                                                                              <  2:886852/67‑1 (MQ=255)                                                                                          tGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTa                                                                                              >  1:886852/1‑67 (MQ=255)                                                                                            cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGgag                                                                        <  1:1121505/87‑1 (MQ=255)                                                                                            cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGgag                                                                        >  2:1121505/1‑87 (MQ=255)                                                                                            cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                    <  1:1694327/139‑1 (MQ=255)                                                                                            cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                    >  2:1694327/1‑139 (MQ=255)                                                                                            cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                    <  2:151273/139‑1 (MQ=255)                                                                                            cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                    <  1:587529/139‑1 (MQ=255)                                                                                                  tCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt              <  2:291808/139‑1 (MQ=255)                                                                                                         gACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCg                                                                                                >  2:913122/1‑50 (MQ=255)                                                                                                         gACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCg                                                                                                <  1:913122/50‑1 (MQ=255)                                                                                                              cACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa  >  1:270766/1‑139 (MQ=255)                                                                                                               aCCGTCACTTCACATAACTTCAAATTCTTCCCACATAg                                                                                                      <  2:1711355/38‑1 (MQ=255)                                                                                                               aCCGTCACTTCACATAACTTCAAATTCTTCCCACATAg                                                                                                      >  1:1711355/1‑38 (MQ=255)                                                                                                                     aCTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt     <  1:1247984/129‑1 (MQ=255)                                                                                                                     aCTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt     >  2:1247984/1‑129 (MQ=255)                                                                                                                         cacaTAACTTCAAATTCTTCCCACATAGTCTTCGTa                                                                                              >  1:872717/1‑36 (MQ=255)                                                                                                                         cacaTAACTTCAAATTCTTCCCACATAGTCTTCGTa                                                                                              <  2:872717/36‑1 (MQ=255)                                                                                                                           caTAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccccc                                                   <  1:196268/77‑1 (MQ=255)                                                                                                                           caTAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccccc                                                   >  2:196268/1‑77 (MQ=255)                                                                                                                               |                                                                                                                          AAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGA  >  NC_000913/4637994‑4638241

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