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%	124.8 / NA	35	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 (14/21);  total (14/21)

GACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATG  >  NC_000913/4638098‑4638252                                                                              |                                                                              gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTa                                                                            >  2:940445/1‑81 (MQ=255) gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTa                                                                            >  1:940478/1‑81 (MQ=255)      cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGtt                                                                        <  1:940428/80‑1 (MQ=255)      cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGtt                                                                        <  2:940427/80‑1 (MQ=255)        ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGa                                                                      >  2:940484/1‑80 (MQ=255)                     tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTc                                                        >  2:940481/1‑81 (MQ=255)                       gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCt                                                       <  1:940430/80‑1 (MQ=255)                       gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCt                                                       <  2:940429/80‑1 (MQ=255)                        cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTg                                                      >  2:940486/1‑80 (MQ=255)                                tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGAc                                             >  1:940503/1‑81 (MQ=255)                                    cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACt                                            <  1:940431/78‑1 (MQ=255)                                    cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAg                                          <  1:940432/80‑1 (MQ=255)                                     ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGcc                                        <  2:940433/81‑1 (MQ=255)                                     ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGcc                                        >  1:940474/1‑81 (MQ=255)                                     ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGc                                         <  2:940434/80‑1 (MQ=255)                                      cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCt                                       <  1:940435/81‑1 (MQ=255)                                      cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCt                                       >  2:940482/1‑81 (MQ=255)                                         aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                    >  2:940488/1‑81 (MQ=255)                                             tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                <  2:940436/81‑1 (MQ=255)                                                 cGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                              <  2:940437/79‑1 (MQ=255)                                                  gTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                           <  1:940438/81‑1 (MQ=255)                                                  gTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                           <  1:940439/81‑1 (MQ=255)                                                   tATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                          >  2:940491/1‑81 (MQ=255)                                                     tCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                        <  2:940440/81‑1 (MQ=255)                                                     tCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                         <  1:940441/80‑1 (MQ=255)                                                     tCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCCGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                        >  2:940480/1‑81 (MQ=255)                                                         gctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                    >  2:940485/1‑81 (MQ=255)                                                          ctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                   <  1:940443/81‑1 (MQ=255)                                                          ctgcCATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                     <  1:940442/79‑1 (MQ=255)                                                             ccATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                 <  1:940445/80‑1 (MQ=255)                                                             ccATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                  <  1:940444/79‑1 (MQ=255)                                                              cATTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt               <  2:940446/81‑1 (MQ=255)                                                               aTTGCAAAGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt               >  2:940499/1‑80 (MQ=255)                                                                    aaaGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc         <  2:940447/81‑1 (MQ=255)                                                                     aaGGAGAAGGCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc         >  2:940505/1‑80 (MQ=255)                                                                         agaaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa    >  2:940487/1‑81 (MQ=255)                                                                         agaaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa    <  2:940448/81‑1 (MQ=255)                                                                         agaaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa    >  2:940490/1‑81 (MQ=255)                                                                          gaaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt   >  1:940489/1‑81 (MQ=255)                                                                          gaaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt   >  1:940495/1‑81 (MQ=255)                                                                           aaggCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg  <  2:940449/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 < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑