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%	134.9 / NA	42	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 (20/22);  total (20/22)

ACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAG  >  NC_000913/2317894‑2318143                                                                                                                              |                                                                                                                             acaAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATaa                                                                                                                 <  1:292598/139‑1 (MQ=255)     aTCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGc                                                                                                             >  2:487395/1‑139 (MQ=255)     aTCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGc                                                                                                             >  1:1317199/1‑139 (MQ=255)                     cGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                             <  1:487395/139‑1 (MQ=255)                       gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGTTCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATaga                                                                                                                        <  2:811989/110‑1 (MQ=255)                       gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCGCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACa                                                                                               <  2:923674/135‑1 (MQ=255)                       gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACa                                                                                               >  1:923674/1‑135 (MQ=255)                       gTCCCCTGGAAATTACGCTGTACGCCCGTTCCAACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATaga                                                                                                                        >  1:811989/1‑110 (MQ=255)                            cTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACt                                                                                      >  1:633841/1‑139 (MQ=255)                             tGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACta                                                                                     >  1:190346/1‑139 (MQ=255)                                         tGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                  >  2:249359/1‑114 (MQ=255)                                         tGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                  <  1:249359/114‑1 (MQ=255)                                           tACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                             <  1:110363/117‑1 (MQ=255)                                           tACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                             >  2:110363/1‑117 (MQ=255)                                            aCGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                             <  1:1312431/116‑1 (MQ=255)                                            aCGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                             >  2:1312431/1‑116 (MQ=255)                                            aCGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCgg                                                                      >  1:349797/1‑139 (MQ=255)                                             cGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGt                                                                     <  2:349797/139‑1 (MQ=255)                                                     ccGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTa                                                             >  2:1032010/1‑139 (MQ=255)                                                        aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                          <  2:683214/139‑1 (MQ=255)                                                             gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTa                                                     >  2:1085744/1‑139 (MQ=255)                                                               agaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACt                                                   <  1:71069/139‑1 (MQ=255)                                                                  agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCcgcg                                                                                                     >  1:785735/1‑86 (MQ=255)                                                                  agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCcgcg                                                                                                     <  2:785735/86‑1 (MQ=255)                                                                       tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGt                                           <  2:661753/139‑1 (MQ=255)                                                                          aaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                  <  1:647112/81‑1 (MQ=255)                                                                          aaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGa                                                                                                  >  2:647112/1‑81 (MQ=255)                                                                               gcgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATg                                   >  1:162815/1‑139 (MQ=255)                                                                                   accacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACct                               >  2:701061/1‑139 (MQ=255)                                                                                     cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGc                                                                                            <  2:1207882/76‑1 (MQ=255)                                                                                     cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGc                                                                                            >  1:1207882/1‑76 (MQ=255)                                                                                                    tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCcgcg                                                                                                     >  1:522276/1‑52 (MQ=255)                                                                                                    tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCcgcg                                                                                                     <  2:522276/52‑1 (MQ=255)                                                                                                    tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                             >  1:1078661/1‑60 (MQ=255)                                                                                                    tATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                             <  2:1078661/60‑1 (MQ=255)                                                                                                     aTCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGt             <  2:1317199/139‑1 (MQ=255)                                                                                                       ccccACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCg           <  1:701061/139‑1 (MQ=255)                                                                                                       ccccACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCg           <  1:524636/139‑1 (MQ=255)                                                                                                               cGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTaa   <  2:162815/139‑1 (MQ=255)                                                                                                                ggTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGTATCGGGTCGCCATTTAAg  <  2:633841/139‑1 (MQ=255)                                                                                                                         aaCACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaaaa                                                                                        <  1:1245615/44‑1 (MQ=255)                                                                                                                         aaCACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCaaaa                                                                                        >  2:1245615/1‑44 (MQ=255)                                                                                                                              |                                                                                                                             ACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAG  >  NC_000913/2317894‑2318143

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