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%	185.1 / NA	57	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 (28/29);  total (28/29)

GTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGC  >  NC_000913/2317887‑2318152                                                                                                                                     |                                                                                                                                      gTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATaga                                                                                                                                 <  2:1402420/139‑1 (MQ=255)            attgCCAGACCGAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATaga                                                                                                                                 >  1:1062079/4‑128 (MQ=255)            aTCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATaga                                                                                                                                 <  2:1062079/128‑1 (MQ=255)            aTCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGc                                                                                                                      >  1:1078863/1‑139 (MQ=255)                      ccagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCcgcgcg                                                                                                            <  1:1784809/139‑1 (MQ=255)                       cagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           <  2:99582/139‑1 (MQ=255)                       cagaccGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           <  2:173820/139‑1 (MQ=255)                            cGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  2:1078863/139‑1 (MQ=255)                            cGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  2:829156/139‑1 (MQ=255)                              gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATag                                                                                                                                  <  1:1400791/109‑1 (MQ=255)                              gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATag                                                                                                                                  >  2:1400791/1‑109 (MQ=255)                              gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCa                                                                                                    >  1:1409574/1‑139 (MQ=255)                              gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCa                                                                                                    >  2:1609358/1‑139 (MQ=255)                              gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCa                                                                                                    >  2:1833155/1‑139 (MQ=255)                              gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCa                                                                                                    >  1:319771/1‑139 (MQ=255)                               tCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaa                                                                                                   >  1:875297/1‑139 (MQ=255)                                        aaTTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACa                                                                                          <  2:875297/139‑1 (MQ=255)                                                tGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTAt                                                                                  >  2:683549/1‑139 (MQ=255)                                                      cccGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGa                                                                            <  1:1609358/139‑1 (MQ=255)                                                               aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaa                                                                                                   <  2:1640967/107‑1 (MQ=255)                                                               aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaa                                                                                                   >  1:1640967/1‑107 (MQ=255)                                                                  tGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGt                                                                <  1:683549/139‑1 (MQ=255)                                                                         agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGAt                                                                                                                            <  1:1241665/72‑1 (MQ=255)                                                                         agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGAt                                                                                                                            >  2:1241665/1‑72 (MQ=255)                                                                         agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           <  1:1459011/89‑1 (MQ=255)                                                                         agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           >  2:1459011/1‑89 (MQ=255)                                                                         agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCaa                                                         <  2:494739/139‑1 (MQ=255)                                                                              tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTa                                                              >  1:593745/1‑129 (MQ=255)                                                                              tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTa                                                              <  2:593745/129‑1 (MQ=255)                                                                                      gcgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGc                                                                                        >  1:1327423/1‑95 (MQ=255)                                                                                      gcgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGc                                                                                        <  2:1327423/95‑1 (MQ=255)                                                                                       cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCg                                                                 <  1:405597/117‑1 (MQ=255)                                                                                       cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCg                                                                 >  2:405597/1‑117 (MQ=255)                                                                                           ccacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCTCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAAc                                                        <  1:1194521/122‑1 (MQ=255)                                                                                           ccacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAAc                                                        >  2:1194521/1‑122 (MQ=255)                                                                                            cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTGTATGGCGTTACTCAACAGGTTGGAGATGACctgc                                      >  2:1176970/1‑139 (MQ=255)                                                                                               ttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCcgcg                                                                                                              >  1:7634/1‑64 (MQ=255)                                                                                               ttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCcgcg                                                                                                              <  2:7634/64‑1 (MQ=255)                                                                                                 cttTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           >  1:1575342/1‑65 (MQ=255)                                                                                                 cttTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                           <  2:1575342/65‑1 (MQ=255)                                                                                                    tCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACg                              >  2:149244/1‑139 (MQ=255)                                                                                                     cGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGc                             <  1:149244/139‑1 (MQ=255)                                                                                                          gTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      >  1:1532203/1‑61 (MQ=255)                                                                                                          gTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  2:1532203/61‑1 (MQ=255)                                                                                                           tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  2:1237387/60‑1 (MQ=255)                                                                                                           tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      >  1:1237387/1‑60 (MQ=255)                                                                                                           tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATTCAGCCGGTATCGGTGAATTTTATGGCGTTACTCaa                                                         <  2:1291602/105‑1 (MQ=255)                                                                                                           tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGACTTTTATGGCGTTACTCaa                                                         >  1:1291602/1‑105 (MQ=255)                                                                                                           tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCg                                                                 >  2:1079475/1‑97 (MQ=255)                                                                                                           tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCg                                                                 <  1:1079475/97‑1 (MQ=255)                                                                                                           tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCggg                       >  2:1573571/1‑139 (MQ=255)                                                                                                              ccccACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCg                    <  1:1833155/139‑1 (MQ=255)                                                                                                                        gTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      >  2:62338/1‑47 (MQ=255)                                                                                                                        gTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                      <  1:62338/47‑1 (MQ=255)                                                                                                                              cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGgag                                               <  1:918038/96‑1 (MQ=255)                                                                                                                              cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGgag                                               >  2:918038/1‑96 (MQ=255)                                                                                                                                aaCACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGc  >  1:1617767/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                      GTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGC  >  NC_000913/2317887‑2318152

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: ‑