Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	652,063	A→G	intergenic (‑207/‑172)	citC ← / → dpiB	citrate lyase synthetase/sensory histidine kinase DpiB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	652,063	0	A	G	100.0%	103.8 / NA	27	intergenic (‑207/‑172)	citC/dpiB	citrate lyase synthetase/sensory histidine kinase DpiB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/12);  total (15/12)

ATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTA  >  NC_000913/651926‑652190                                                                                                                                          |                                                                                                                                aTAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGAt                                                                                                                               >  1:387240/1‑140 (MQ=255)              ttCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGtttt                                                                                                                  <  2:846963/140‑1 (MQ=255)                   ttGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGATATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGTTGATTTTAAGTTTTttaat                                                                                                             >  1:138949/1‑140 (MQ=255)                   ttGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTttaat                                                                                                             <  2:504534/140‑1 (MQ=255)                   ttGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTttaat                                                                                                             <  2:147306/140‑1 (MQ=255)                                ggTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACg                                                                                                >  2:148955/1‑140 (MQ=255)                                   gATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGaaa                                                                                             >  1:284107/1‑140 (MQ=255)                                    aTCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAAt                                                                                            >  2:416534/1‑140 (MQ=255)                                         ggTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTc                                                                                       <  1:148955/140‑1 (MQ=255)                                          gTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGTCTCg                                                                                      <  2:138949/140‑1 (MQ=255)                                              atCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAgg                                                                                  >  1:144593/1‑140 (MQ=255)                                              atCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAgg                                                                                  >  1:64112/1‑140 (MQ=255)                                                   tttCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTaa                                                                             <  1:352347/140‑1 (MQ=255)                                                   tttCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTaa                                                                             >  2:331627/1‑140 (MQ=255)                                                            tCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGtgatttaa                                                                    <  1:331627/140‑1 (MQ=255)                                                                       ccTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATg                                                         >  2:122826/1‑140 (MQ=255)                                                                            ggAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATaa                                                    >  1:366286/1‑140 (MQ=255)                                                                              aGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATaaaa                                                  >  1:676614/1‑140 (MQ=255)                                                                                tATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAAtt                                                >  2:224160/1‑140 (MQ=255)                                                                                  tttCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTg                                              >  2:334400/1‑140 (MQ=255)                                                                                   ttCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGc                                             <  2:239959/140‑1 (MQ=255)                                                                                     cGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCa                                           <  2:387240/140‑1 (MQ=255)                                                                                          tGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATc                                      >  1:195866/1‑140 (MQ=255)                                                                                                 ttGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCAtaatta                               <  1:416534/140‑1 (MQ=255)                                                                                                    atttatttCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAata                            <  2:284107/140‑1 (MQ=255)                                                                                                      ttatttCACATAGAGTTAGTGGTTTTTTATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATc                          <  2:366286/140‑1 (MQ=255)                                                                                                                              ttttATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTa  >  1:51887/1‑140 (MQ=255)                                                                                                                              ttttATTTATTTGATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTa  >  1:125892/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                ATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTTTAAGTGATTTAATTGATTTAATGAATAAAATTTGCCACGATCATAATTAATATCTATGTATTTTGATTCAACATTTTA  >  NC_000913/651926‑652190

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