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%	205.3 / NA	64	G90S (GGC→AGC)	bdcR	transcriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/29);  total (35/29)

CCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCA  >  NC_000913/4474270‑4474518                                                                                                                          |                                                                                                                                ccccccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                                <  1:1057939/139‑1 (MQ=255)   ccccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGt                                                                                                              >  2:560895/1‑139 (MQ=255)    cccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTa                                                                                                             >  2:973949/1‑139 (MQ=255)    cccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTa                                                                                                             >  1:560613/1‑139 (MQ=255)            tctACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                         <  2:123572/119‑1 (MQ=255)            tctACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                         >  1:123572/1‑119 (MQ=255)                cgcgGCTTTTGGCGGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGa                                                                                                                             <  1:418731/111‑1 (MQ=255)                cgcgGCTTTTGGCGGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGa                                                                                                                             >  2:418731/1‑111 (MQ=255)                 gcgGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcgg                                                                                                <  2:716072/139‑1 (MQ=255)                   ggCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcggcg                                                                                              >  1:283010/1‑139 (MQ=255)                     cTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                            <  2:283010/139‑1 (MQ=255)                                       ggTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTAAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAAcgg                                                                          >  2:894369/1‑139 (MQ=255)                                           aTTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACggcggc                                                                      <  2:560613/139‑1 (MQ=255)                                                   gtgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctgg                                                              <  1:405488/139‑1 (MQ=255)                                                      tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcgc                                                                 <  1:122703/133‑1 (MQ=255)                                                      tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcgc                                                                 >  2:122703/1‑133 (MQ=255)                                                      tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctggctg                                                           >  1:787587/1‑139 (MQ=255)                                                              gAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                            >  2:1208870/1‑98 (MQ=255)                                                              gAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                            <  1:1208870/98‑1 (MQ=255)                                                                aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcgc                                                                 <  1:919645/123‑1 (MQ=255)                                                                aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcgc                                                                 >  2:919645/1‑123 (MQ=255)                                                                aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCtt                                                 >  2:555891/1‑139 (MQ=255)                                                                 tACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTg                                                <  1:560895/139‑1 (MQ=255)                                                                   cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGTTCGTCCGGTAAGCGa                                                                                                                             >  2:929345/1‑60 (MQ=255)                                                                   cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGa                                                                                                                             <  1:929345/60‑1 (MQ=255)                                                                   cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                                >  1:391673/1‑73 (MQ=255)                                                                   cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                                <  2:391673/73‑1 (MQ=255)                                                                   cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGaa                                              >  1:40818/1‑139 (MQ=255)                                                                       ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGat                                                                                       >  2:1211915/1‑94 (MQ=255)                                                                       ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGat                                                                                       <  1:1211915/94‑1 (MQ=255)                                                                        gTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                            <  1:1302682/88‑1 (MQ=255)                                                                        gTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                            >  2:1302682/1‑88 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATg                                                       <  1:918884/124‑1 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATg                                                       >  2:918884/1‑124 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTAtt                                        >  2:104098/1‑139 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTAtt                                        >  1:1118665/1‑139 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTAtt                                        >  1:434863/1‑139 (MQ=255)                                                                                 cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                                <  1:973949/139‑1 (MQ=255)                                                                                 cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                                <  2:453213/139‑1 (MQ=255)                                                                                   aTTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATg                              >  2:1044054/1‑139 (MQ=255)                                                                                        gCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCca                         >  1:955536/1‑139 (MQ=255)                                                                                           tGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                                 <  1:1298683/64‑1 (MQ=255)                                                                                           tGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                                 >  2:1298683/1‑64 (MQ=255)                                                                                              cGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAgcg                   >  2:349033/1‑139 (MQ=255)                                                                                              cGATATTCTTCGTGATAATCGTCCGGTAAGCGAGTGCCTGGTTAAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAgcg                   <  1:894369/139‑1 (MQ=255)                                                                                                   ttcttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGa              >  2:198327/1‑139 (MQ=255)                                                                                                      ttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGt                                                                                                              >  2:1140778/1‑40 (MQ=255)                                                                                                      ttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGt                                                                                                              <  1:1140778/40‑1 (MQ=255)                                                                                                        cGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                            >  2:147199/1‑56 (MQ=255)                                                                                                        cGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                            <  1:147199/56‑1 (MQ=255)                                                                                                        cGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                  <  1:836043/98‑1 (MQ=255)                                                                                                        cGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                  >  2:836043/1‑98 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACggcggc                                                                      >  2:312833/1‑73 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACggcggc                                                                      <  1:312833/73‑1 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGtt    <  2:434863/139‑1 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGtt    <  1:555891/139‑1 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGtt    >  2:815011/1‑139 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGtt    <  1:104098/139‑1 (MQ=255)                                                                                                               atCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCa  >  2:1188559/1‑139 (MQ=255)                                                                                                               atCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGCCATGATCCACAAGCGCGTGATATTGCCGTTCa  >  2:155235/1‑139 (MQ=255)                                                                                                                tCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAg                                             <  2:1032230/95‑1 (MQ=255)                                                                                                                tCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAg                                             >  1:1032230/1‑95 (MQ=255)                                                                                                                  gTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctgg                                                              <  2:132735/76‑1 (MQ=255)                                                                                                                  gTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctgg                                                              >  1:132735/1‑76 (MQ=255)                                                                                                                          |                                                                                                                                CCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCA  >  NC_000913/4474270‑4474518

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