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

CCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTT  >  NC_000913/4474270‑4474516                                                                                                                          |                                                                                                                              ccccccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGGATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              <  2:852131/139‑1 (MQ=255)     ccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTAt                                                                                                          >  2:214191/1‑139 (MQ=255)     ccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTAt                                                                                                          >  2:1515521/1‑139 (MQ=255)      cGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTAtt                                                                                                         >  1:1190542/1‑139 (MQ=255)         gCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATaaat                                                                                                      <  1:1515521/139‑2 (MQ=255)            tctACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                       <  1:408738/119‑1 (MQ=255)            tctACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                       >  2:408738/1‑119 (MQ=255)                cgcgGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               <  1:1687846/139‑1 (MQ=255)                    gCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgc                                                                                           <  1:214191/139‑1 (MQ=255)                       tttGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCAAGTGCCTGGTTGAGGTATTAAAAGAAgc                                                                                                <  2:1597462/131‑1 (MQ=255)                       tttGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCAAGTGCCTGGTTGAGGTATTAAAAGAAgc                                                                                                >  1:1597462/1‑131 (MQ=255)                                            tttAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCt                                                                   >  1:1103899/1‑139 (MQ=255)                                                 ccGTGTACTCAATGCATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCAGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          >  2:1096842/1‑111 (MQ=255)                                                 ccGTGTACTCAAGGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  1:1096842/111‑1 (MQ=255)                                                    tgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAAcggcgg                                                                     >  2:448816/1‑129 (MQ=255)                                                    tgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAAcggcgg                                                                     <  1:448816/129‑1 (MQ=255)                                                      tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctggctg                                                         >  2:396422/1‑139 (MQ=255)                                                           aaTGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTg                                                                                                                        <  1:189630/71‑1 (MQ=255)                                                           aaTGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTg                                                                                                                        >  2:189630/1‑71 (MQ=255)                                                               aaTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                <  2:1374671/139‑1 (MQ=255)                                                                aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              >  1:318588/1‑76 (MQ=255)                                                                aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              <  2:318588/76‑1 (MQ=255)                                                                     tCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGa                                                                                                               >  2:445884/1‑70 (MQ=255)                                                                     tCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGa                                                                                                               <  1:445884/70‑1 (MQ=255)                                                                     tCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                >  1:128088/1‑117 (MQ=255)                                                                     tCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                <  2:128088/117‑1 (MQ=255)                                                                     tCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAgg                                          >  2:70768/1‑139 (MQ=255)                                                                       ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctgg                                                            <  2:676173/119‑1 (MQ=255)                                                                       ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctgg                                                            >  1:676173/1‑119 (MQ=255)                                                                       ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTa                                        >  2:150659/1‑139 (MQ=255)                                                                       ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTa                                        >  2:723053/1‑139 (MQ=255)                                                                       ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTa                                        >  2:1207946/1‑139 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAgg                                          >  1:1244597/1‑135 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAgg                                          <  2:1244597/135‑1 (MQ=255)                                                                         tACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTAtt                                      >  1:935607/1‑139 (MQ=255)                                                                                 cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  2:935607/139‑1 (MQ=255)                                                                                 cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  1:1207946/139‑1 (MQ=255)                                                                                 cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  1:150659/139‑1 (MQ=255)                                                                                   aTTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATGGCCAAAACGGCGGCTGCgct                                                              >  2:230510/1‑105 (MQ=255)                                                                                   aTTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATGGCCAAAACGGCGGCTGCgct                                                              <  1:230510/105‑1 (MQ=255)                                                                                          ttGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGa                                                                                                               >  1:1402890/1‑49 (MQ=255)                                                                                          ttGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGa                                                                                                               <  2:1402890/49‑1 (MQ=255)                                                                                          ttGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          >  2:1175306/1‑70 (MQ=255)                                                                                          ttGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  1:1175306/70‑1 (MQ=255)                                                                                          ttGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                <  1:579715/112‑1 (MQ=255)                                                                                          ttGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                >  2:579715/1‑112 (MQ=255)                                                                                            gCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAaagaag                                                                                                 <  1:226059/61‑1 (MQ=255)                                                                                            gCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAaagaag                                                                                                 >  2:226059/1‑61 (MQ=255)                                                                                              cGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAgcg                 >  1:1587994/1‑139 (MQ=255)                                                                                                    tcttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              >  2:1308228/1‑40 (MQ=255)                                                                                                    tcttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              <  1:1308228/40‑1 (MQ=255)                                                                                                      ttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAgcg                 >  2:209118/1‑131 (MQ=255)                                                                                                      ttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAgcg                 <  1:209118/131‑1 (MQ=255)                                                                                                           gatgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATgg                                                    <  2:111919/91‑1 (MQ=255)                                                                                                           gatgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATgg                                                    >  1:111919/1‑91 (MQ=255)                                                                                                           gatgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                <  1:512603/95‑1 (MQ=255)                                                                                                           gatgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                >  2:512603/1‑95 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcgg                                                                                              >  1:473759/1‑47 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcgg                                                                                              <  2:473759/47‑1 (MQ=255)                                                                                                             tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGtt  <  1:70768/139‑1 (MQ=255)                                                                                                                tCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  2:1217522/108‑1 (MQ=255)                                                                                                                tCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              >  1:1217522/1‑108 (MQ=255)                                                                                                                  gTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               <  2:1027132/41‑1 (MQ=255)                                                                                                                  gTCCGGTAAGCGAGCGCCTGGTTGAGGTATTAAAAGAAgcg                                                                                               >  1:1027132/1‑41 (MQ=255)                                                                                                                          |                                                                                                                              CCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTT  >  NC_000913/4474270‑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: ‑