Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,318,019	G→A	R620C (CGT→TGT)	rcsC ←	hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,318,019	0	G	A	100.0%	131.1 / NA	41	R620C (CGT→TGT)	rcsC	hybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/20);  total (21/20)

ATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGC  >  NC_000913/2317883‑2318152                                                                                                                                         |                                                                                                                                      aTCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGa                                                                                                                                     >  1:1184089/1‑139 (MQ=255)   cAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATa                                                                                                                                   <  2:1058758/139‑1 (MQ=255)      tttttCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACACTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATagag                                                                                                                                >  2:1237416/1‑139 (MQ=255)         ttCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAgagaga                                                                                                                             <  1:1457040/139‑1 (MQ=255)                           cagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAAACGCCATCCGCGCGa                                                                                                           <  2:1155512/139‑1 (MQ=255)                                cGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  1:997989/139‑1 (MQ=255)                                                      tACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCg                                                                                >  1:901805/1‑139 (MQ=255)                                                         gCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaa                                                                                                   <  2:369252/117‑1 (MQ=255)                                                         gCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaa                                                                                                   >  1:369252/1‑117 (MQ=255)                                                              ttCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAAttt                                                                        >  2:1687138/1‑139 (MQ=255)                                                                   aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                   <  1:1687138/139‑1 (MQ=255)                                                                        gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  1:1102069/99‑1 (MQ=255)                                                                        gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      >  2:1102069/1‑99 (MQ=255)                                                                        gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTa                                                              <  1:1237416/139‑1 (MQ=255)                                                                          agaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACt                                                            <  2:1184089/139‑1 (MQ=255)                                                                               ggATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACa                                                       >  1:653563/1‑139 (MQ=255)                                                                                  tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAgagaga                                                                                                                             >  1:546129/1‑66 (MQ=255)                                                                                  tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAgagaga                                                                                                                             <  2:546129/66‑1 (MQ=255)                                                                                  tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGt                                                    >  2:79652/1‑139 (MQ=255)                                                                                     aaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTgg                                                 >  1:955216/1‑139 (MQ=255)                                                                                           cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGc                                                                                        <  1:531141/94‑1 (MQ=255)                                                                                           cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGc                                                                                        >  2:531141/1‑94 (MQ=255)                                                                                                  cttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACctgctg                                    >  2:1362572/1‑139 (MQ=255)                                                                                                        tCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  2:319103/67‑1 (MQ=255)                                                                                                        tCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      >  1:319103/1‑67 (MQ=255)                                                                                                        tCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACg                              >  1:1124160/1‑139 (MQ=255)                                                                                                          gcgggTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCa                            <  2:1124160/136‑1 (MQ=255)                                                                                                               tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                   >  1:903056/1‑95 (MQ=255)                                                                                                               tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                   <  2:903056/95‑1 (MQ=255)                                                                                                                      aCGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           >  2:362777/1‑48 (MQ=255)                                                                                                                      aCGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           <  1:362777/48‑1 (MQ=255)                                                                                                                      aCGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCAt                >  1:224428/1‑139 (MQ=255)                                                                                                                        gCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                   <  1:726646/86‑1 (MQ=255)                                                                                                                        gCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                   >  2:726646/1‑86 (MQ=255)                                                                                                                        gCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCAttt              >  1:474594/1‑139 (MQ=255)                                                                                                                           ggTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaa                                                                                                   >  2:1312614/1‑51 (MQ=255)                                                                                                                           ggTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaa                                                                                                   <  1:1312614/51‑1 (MQ=255)                                                                                                                               tCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCa       <  2:474594/139‑1 (MQ=255)                                                                                                                                  cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                   <  2:1649174/76‑1 (MQ=255)                                                                                                                                  cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                   >  1:1649174/1‑76 (MQ=255)                                                                                                                                    aaCACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGc  <  2:955216/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      ATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGC  >  NC_000913/2317883‑2318152

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