Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,474,391	G→A	G90S (GGC→AGC)	bdcR →	transcriptional repressor for divergent bdcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,474,391	0	G	A	100.0%	198.3 / NA	61	G90S (GGC→AGC)	bdcR	transcriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/32);  total (29/32)

TTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTT  >  NC_000913/4474255‑4474516                                                                                                                                         |                                                                                                                              ttATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGc                                                                                                                             <  1:563832/139‑1 (MQ=255)       tGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                       <  2:508161/139‑1 (MQ=255)              aaCCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                >  2:1029396/1‑139 (MQ=255)              aaCCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                <  1:382228/139‑1 (MQ=255)                   cccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTa                                                                                                           >  2:263101/1‑139 (MQ=255)                    ccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTAt                                                                                                          >  1:610999/1‑139 (MQ=255)                      gAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTa                                                                                                        >  1:358440/1‑139 (MQ=255)                           tctACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAaaga                                                                                                   >  2:551925/1‑139 (MQ=255)                                    cTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  2:610999/139‑1 (MQ=255)                                              tAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGc                                                                                >  1:420238/1‑139 (MQ=255)                                               aaaGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGcc                                                                               <  2:420238/139‑1 (MQ=255)                                                               gCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcgc                                                               <  2:157069/139‑1 (MQ=255)                                                                ccGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                       >  1:1112245/1‑82 (MQ=255)                                                                ccGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                       <  2:1112245/82‑1 (MQ=255)                                                                ccGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgct                                                              >  2:203417/1‑139 (MQ=255)                                                                ccGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgct                                                              <  2:780515/139‑1 (MQ=255)                                                                  gtgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              >  2:1042598/1‑89 (MQ=255)                                                                  gtgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              <  1:1042598/89‑1 (MQ=255)                                                                  gtgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctgg                                                            <  2:358440/139‑1 (MQ=255)                                                                   tgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  2:166773/108‑1 (MQ=255)                                                                   tgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          >  1:166773/1‑108 (MQ=255)                                                                     tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                >  2:748711/1‑132 (MQ=255)                                                                     tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                <  1:748711/132‑1 (MQ=255)                                                                          aaTGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATgg                                                    >  2:30366/1‑139 (MQ=255)                                                                           aTGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          >  2:528757/1‑100 (MQ=255)                                                                           aTGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  1:528757/100‑1 (MQ=255)                                                                             gAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTc                                                 >  2:429503/1‑139 (MQ=255)                                                                              aaTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                <  2:409958/139‑1 (MQ=255)                                                                              aaTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                <  1:203417/139‑1 (MQ=255)                                                                               aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               >  1:560703/1‑91 (MQ=255)                                                                               aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               <  2:560703/91‑1 (MQ=255)                                                                                  cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                >  2:1092372/1‑71 (MQ=255)                                                                                  cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                <  1:1092372/71‑1 (MQ=255)                                                                                  cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               >  2:368750/1‑88 (MQ=255)                                                                                  cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               <  1:368750/88‑1 (MQ=255)                                                                                      ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               >  2:1001390/1‑84 (MQ=255)                                                                                      ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               <  1:1001390/84‑1 (MQ=255)                                                                                          cGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                <  2:837588/63‑1 (MQ=255)                                                                                          cGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                >  1:837588/1‑63 (MQ=255)                                                                                             aaGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGTGCTGGCTGTATGGTTCTTGAAGGTATTCATAg                                 <  2:660662/139‑1 (MQ=255)                                                                                              aGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGt                                <  2:1029293/139‑1 (MQ=255)                                                                                                cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  1:1029396/139‑1 (MQ=255)                                                                                                cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  1:263101/139‑1 (MQ=255)                                                                                                cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  1:551925/139‑1 (MQ=255)                                                                                                  aTTTCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGGTGAGGTATTAAAAGACGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTAAGGTTCTTGAAGGTATTCATAGTCATg                            >  2:516566/1‑139 (MQ=255)                                                                                                  aTTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATg                            >  2:1114851/1‑139 (MQ=255)                                                                                                        cTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGcc                                                                                                                      >  1:53645/1‑43 (MQ=255)                                                                                                        cTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGcc                                                                                                                      <  2:53645/43‑1 (MQ=255)                                                                                                        cTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTAt                                                      >  1:368012/1‑107 (MQ=255)                                                                                                        cTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTAt                                                      <  2:368012/107‑1 (MQ=255)                                                                                                           gCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAGACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATAATCCACAAg                   >  2:656956/1‑139 (MQ=255)                                                                                                           gCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                <  1:757094/94‑1 (MQ=255)                                                                                                           gCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                >  2:757094/1‑94 (MQ=255)                                                                                                               atatTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGcgcg               <  1:30366/139‑1 (MQ=255)                                                                                                                tatTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcgc                                                               <  2:428103/90‑1 (MQ=255)                                                                                                                tatTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcgc                                                               >  1:428103/1‑90 (MQ=255)                                                                                                                       cGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGtt                                                  >  2:166992/1‑96 (MQ=255)                                                                                                                       cGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGtt                                                  <  1:166992/96‑1 (MQ=255)                                                                                                                           atgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcgg                                                                                              <  1:398540/48‑1 (MQ=255)                                                                                                                           atgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcgg                                                                                              >  2:398540/1‑48 (MQ=255)                                                                                                                            tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGtt  >  1:371264/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                              TTATCTTGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTT  >  NC_000913/4474255‑4474516

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