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%	125.5 / NA	40	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 (18/22);  total (18/22)

ATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTG  >  NC_000913/2317883‑2318150                                                                                                                                         |                                                                                                                                    aTCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTCGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGa                                                                                                                                   <  1:1015862/139‑1 (MQ=255)   cAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATa                                                                                                                                 >  2:419834/1‑139 (MQ=255)      tttttCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATagag                                                                                                                              >  1:919136/1‑139 (MQ=255)         ttCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAgagaga                                                                                                                           <  2:569785/139‑1 (MQ=255)                                  gTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCa                                                                                                  >  1:1316317/1‑139 (MQ=255)                                        tGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACta                                                                                            <  1:72888/139‑1 (MQ=255)                                                        cGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGt                                                                            >  2:240395/1‑139 (MQ=255)                                                                ccGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    >  2:641029/1‑107 (MQ=255)                                                                ccGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    <  1:641029/107‑1 (MQ=255)                                                                   aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                 <  2:513536/139‑1 (MQ=255)                                                                   aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                 <  1:1349444/139‑1 (MQ=255)                                                                   aCCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATgg                                                                 <  1:240395/139‑1 (MQ=255)                                                                     cTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCg                                                               <  2:919136/139‑1 (MQ=255)                                                                      gggAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGt                                                              <  1:1386831/138‑1 (MQ=255)                                                                        gaagaaGGGATCAACCAAGCGCTCCACTTCTTTCACCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAAATAATCGCCATCCGCGCGAACATACAAAAGTATACAGCCGGTATCGGTGAATTTTATGGCGTTa                                                            <  1:1572065/139‑1 (MQ=255)                                                                        gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTa                                                            >  2:858243/1‑139 (MQ=255)                                                                             agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCaa                                                       <  1:858243/139‑1 (MQ=255)                                                                             agGGATCAAACAAGCGCACCACTTCTTTCGCCGATATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCaa                                                       <  2:1492462/139‑1 (MQ=255)                                                                                 aTCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAgg                                                   >  2:450868/1‑139 (MQ=255)                                                                                   cAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCaa                                                       >  1:175832/1‑133 (MQ=255)                                                                                   cAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCaa                                                       <  2:175832/133‑1 (MQ=255)                                                                                           cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACtat                                                                                           <  1:1408667/89‑1 (MQ=255)                                                                                           cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACtat                                                                                           >  2:1408667/1‑89 (MQ=255)                                                                                           cgcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGa                                         >  1:132682/1‑139 (MQ=255)                                                                                            gcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCACGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGAc                                        <  1:419834/139‑1 (MQ=255)                                                                                               ccacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    <  2:322679/76‑1 (MQ=255)                                                                                               ccacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    >  1:322679/1‑76 (MQ=255)                                                                                                                aTCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    >  2:1474169/1‑59 (MQ=255)                                                                                                                aTCCCCACGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    <  1:1474169/59‑1 (MQ=255)                                                                                                                      aCGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    <  1:445339/53‑1 (MQ=255)                                                                                                                      aCGCCGGTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    >  2:445339/1‑53 (MQ=255)                                                                                                                            gTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCa                                                                                                  >  2:1176130/1‑49 (MQ=255)                                                                                                                            gTATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCa                                                                                                  <  1:1176130/49‑1 (MQ=255)                                                                                                                             tATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTg                                                                           <  2:1104271/71‑1 (MQ=255)                                                                                                                             tATCGCGAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTg                                                                           >  1:1104271/1‑71 (MQ=255)                                                                                                                                  cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    <  1:385695/41‑1 (MQ=255)                                                                                                                                  cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATg                                                                                                    >  2:385695/1‑41 (MQ=255)                                                                                                                                  cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCgg                      <  2:415771/119‑1 (MQ=255)                                                                                                                                  cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCgg                      >  1:415771/1‑119 (MQ=255)                                                                                                                                  cgAACACAGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATCACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTg  >  1:1014975/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                    ATCAGTTTTTCACAAATCGCCAGACCCAGACCGGTCCCCTGGAAATTACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTG  >  NC_000913/2317883‑2318150

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