Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,175	A→G	intergenic (+55/‑3)	creC → / → creD	sensory histidine kinase CreC/putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,175	0	A	G	100.0%	139.0 / NA	41	intergenic (+55/‑3)	creC/creD	sensory histidine kinase CreC/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/18);  total (23/18)

GACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATG  >  NC_000913/4638098‑4638252                                                                              |                                                                              gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTa                                                                            >  2:964959/1‑81 (MQ=255)   cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATg                                                                          >  1:964953/1‑81 (MQ=255)      cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGtt                                                                        >  2:964957/1‑80 (MQ=255)      cACCGTCACTTCACATAGCTTCAAATTCTTCCAACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGt                                                                         <  2:964909/79‑1 (MQ=255)       aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGa                                                                      <  2:964910/81‑1 (MQ=255)        ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGaa                                                                     <  1:964911/81‑1 (MQ=255)         cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGaaa                                                                    <  1:964912/81‑1 (MQ=255)           tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATc                                                                  <  2:964913/81‑1 (MQ=255)            cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATcc                                                                 >  1:964960/1‑81 (MQ=255)              cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATcccc                                                               >  1:964985/1‑81 (MQ=255)                  acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTg                                                           <  2:964914/81‑1 (MQ=255)                   caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGt                                                          >  2:964958/1‑81 (MQ=255)                    aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGtt                                                         >  1:964969/1‑81 (MQ=255)                    aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGtt                                                         >  2:964951/1‑81 (MQ=255)                    aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGtt                                                         <  2:964915/81‑1 (MQ=255)                       gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCt                                                       <  2:964916/80‑1 (MQ=255)                        cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTgg                                                     <  2:964917/81‑1 (MQ=255)                           cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGa                                                    <  2:964918/79‑1 (MQ=255)                           cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGa                                                    <  1:964919/79‑1 (MQ=255)                              aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAAt                                                >  2:964967/1‑80 (MQ=255)                              aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATg                                               <  1:964920/81‑1 (MQ=255)                                 cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGAc                                             >  2:964966/1‑80 (MQ=255)                                   tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAg                                          >  2:964981/1‑81 (MQ=255)                                     ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGGCTAGc                                         <  1:964921/80‑1 (MQ=255)                                       acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTg                                      <  1:964922/81‑1 (MQ=255)                                       acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTg                                      >  1:964961/1‑81 (MQ=255)                                                 cGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                              <  1:964923/79‑1 (MQ=255)                                                   tATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                          >  1:964971/1‑81 (MQ=255)                                                    aTCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                         >  2:964974/1‑81 (MQ=255)                                                       ctgctgCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                      <  1:964924/81‑1 (MQ=255)                                                          ctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                   >  1:964980/1‑81 (MQ=255)                                                          ctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                   >  1:964975/1‑81 (MQ=255)                                                          ctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                   >  2:964968/1‑81 (MQ=255)                                                            gcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                 >  2:964976/1‑81 (MQ=255)                                                             ccATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                >  2:964982/1‑81 (MQ=255)                                                              cATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt               <  1:964925/81‑1 (MQ=255)                                                                 tGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt            >  2:964972/1‑81 (MQ=255)                                                                     aaGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc         >  2:964964/1‑80 (MQ=255)                                                                     aaGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc         >  1:964991/1‑80 (MQ=255)                                                                       ggAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa       >  2:964983/1‑80 (MQ=255)                                                                       ggAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt      <  2:964926/81‑1 (MQ=255)                                                                           aaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg  >  2:964977/1‑81 (MQ=255)                                                                              |                                                                              GACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATG  >  NC_000913/4638098‑4638252

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