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%	140.8 / NA	38	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 (16/22);  total (16/22)

GACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATG  >  NC_000913/4638098‑4638252                                                                              |                                                                              gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTa                                                                            <  1:873357/81‑1 (MQ=255)  aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTAt                                                                           <  1:873358/81‑1 (MQ=255)       aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTg                                                                       >  1:873397/1‑80 (MQ=255)       aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTACCCTGCTGCCATTGCGAAGGAGAAGGCTATGTTg                                                                       <  2:873359/80‑1 (MQ=255)          gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAAt                                                                   <  2:873360/81‑1 (MQ=255)            cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGTCATTGCAAAGGAGAAGGCTATGTTGAAATcc                                                                 >  2:873402/1‑81 (MQ=255)            cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATcc                                                                 >  2:873400/1‑81 (MQ=255)             aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATccc                                                                >  2:873406/1‑81 (MQ=255)             aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATccc                                                                <  2:873363/81‑1 (MQ=255)             aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATcc                                                                 <  1:873362/80‑1 (MQ=255)             aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATcc                                                                 <  2:873361/80‑1 (MQ=255)                   caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGt                                                          >  1:873398/1‑81 (MQ=255)                     tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTc                                                        <  2:873364/81‑1 (MQ=255)                       gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCt                                                       <  2:873366/80‑1 (MQ=255)                       gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCt                                                       <  2:873365/80‑1 (MQ=255)                        cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTgg                                                     <  2:873367/81‑1 (MQ=255)                              aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATg                                               <  2:873368/81‑1 (MQ=255)                                 cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGAc                                             <  1:873369/80‑1 (MQ=255)                                    cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAg                                          <  2:873370/80‑1 (MQ=255)                                          tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                   <  2:873371/81‑1 (MQ=255)                                           aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                  >  1:873403/1‑81 (MQ=255)                                            gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                 >  1:873408/1‑81 (MQ=255)                                            gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                 >  1:873404/1‑81 (MQ=255)                                              cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                               >  2:873413/1‑81 (MQ=255)                                               ttCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                               >  2:873405/1‑80 (MQ=255)                                                 cGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                             <  2:873372/80‑1 (MQ=255)                                                  gTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                           <  2:873373/81‑1 (MQ=255)                                                   tATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                          >  2:873407/1‑81 (MQ=255)                                                     tCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                        >  2:873410/1‑81 (MQ=255)                                                     tCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                         <  1:873374/80‑1 (MQ=255)                                                      ccTGCTGCCATTGCAAAGGACAAGGGTATGTCGAAATCCCCCCTGTTCTGGAAAACGACTAGCCTGTTTGGCGCAGTATTg                       >  2:873416/1‑81 (MQ=255)                                                        tgctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                     >  1:873418/1‑81 (MQ=255)                                                          ctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                    <  1:873375/80‑1 (MQ=255)                                                             ccATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                >  2:873409/1‑81 (MQ=255)                                                             ccATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                  <  2:873376/79‑1 (MQ=255)                                                                 tGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt            >  2:873412/1‑81 (MQ=255)                                                                   cAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc          <  1:873377/81‑1 (MQ=255)                                                                       ggAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt      <  1:873378/81‑1 (MQ=255)                                                                         agaaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa    >  1:873417/1‑81 (MQ=255)                                                                           aaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg  <  1:873379/81‑1 (MQ=255)                                                                              |                                                                              GACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATG  >  NC_000913/4638098‑4638252

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