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

TGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTT  >  NC_000913/4474261‑4474516                                                                                                                                   |                                                                                                                              tGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGc                                                                                                                       <  1:110266/139‑1 (MQ=255)  ggTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGcc                                                                                                                      <  2:99992/139‑1 (MQ=255)  ggTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGcc                                                                                                                      >  2:395306/1‑139 (MQ=255)    tATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTg                                                                                                                    >  1:629565/1‑139 (MQ=255)    tATCAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGCGGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTg                                                                                                                    >  2:495874/1‑139 (MQ=255)        aaCCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                >  2:377044/1‑139 (MQ=255)        aaCCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTg                                                                                                                >  2:134626/1‑139 (MQ=255)          ccccccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAg                                                                                                              <  2:1583/139‑1 (MQ=255)              ccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTAt                                                                                                          >  1:10757/1‑139 (MQ=255)              ccGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTAt                                                                                                          >  1:115567/1‑139 (MQ=255)                gAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTa                                                                                                        >  2:502259/1‑139 (MQ=255)                              cTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  1:502259/139‑1 (MQ=255)                              cTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  1:477883/139‑1 (MQ=255)                                tttGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGAGCa                                                                                        >  1:86011/1‑139 (MQ=255)                                    gCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGata                                                                                    <  2:10757/139‑1 (MQ=255)                                                        aGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                <  2:629565/139‑1 (MQ=255)                                                        aGCCGTGTACTCAATGAATACGCCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTgcg                                                                <  1:495874/139‑1 (MQ=255)                                                          ccGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgct                                                              >  2:188496/1‑139 (MQ=255)                                                            gtgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          >  2:13943/1‑109 (MQ=255)                                                            gtgtACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          <  1:13943/109‑1 (MQ=255)                                                               tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctggctg                                                         >  1:226016/1‑139 (MQ=255)                                                               tACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCgctggctg                                                         >  2:190094/1‑139 (MQ=255)                                                                 cTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTa                                                       <  1:190094/139‑1 (MQ=255)                                                                   cAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATg                                                     <  2:321669/139‑1 (MQ=255)                                                                    aaTGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATgg                                                    >  1:43645/1‑139 (MQ=255)                                                                       gAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTc                                                 >  1:528242/1‑139 (MQ=255)                                                                       gAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTc                                                 >  2:322946/1‑139 (MQ=255)                                                                        aaTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCt                                                <  1:628228/139‑1 (MQ=255)                                                                         aTACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCtt                                               >  1:296943/1‑139 (MQ=255)                                                                           aCGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGtt                                                                                                                 >  1:433999/1‑71 (MQ=255)                                                                           aCGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGtt                                                                                                                 <  2:433999/71‑1 (MQ=255)                                                                            cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGaa                                            >  1:675727/1‑139 (MQ=255)                                                                            cGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGaa                                            <  2:43645/139‑1 (MQ=255)                                                                              tCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAgg                                          <  2:115567/139‑1 (MQ=255)                                                                                ggTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTa                                        >  1:329806/1‑139 (MQ=255)                                                                                       aaGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAg                                 <  2:226016/139‑1 (MQ=255)                                                                                         gCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTc                               >  1:433764/1‑139 (MQ=255)                                                                                          cTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  1:377044/139‑1 (MQ=255)                                                                                                 gCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              >  2:278389/1‑132 (MQ=255)                                                                                                 gCTTGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCa                              <  1:278389/132‑1 (MQ=255)                                                                                                   ttGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCcaca                     <  2:528242/139‑1 (MQ=255)                                                                                                    tGCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACaa                    <  1:466923/139‑1 (MQ=255)                                                                                                     gCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTAt                                                      <  1:527897/104‑1 (MQ=255)                                                                                                     gCCGATATTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTAt                                                      >  2:527897/1‑104 (MQ=255)                                                                                                         atatTCTTCGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGcgcg               <  1:188496/139‑1 (MQ=255)                                                                                                               ttcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGatat         >  2:518447/1‑139 (MQ=255)                                                                                                                tcGTGATGATCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATAtt        >  1:516856/1‑139 (MQ=255)                                                                                                                    gatgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCg    <  1:518447/139‑1 (MQ=255)                                                                                                                      tgatCGTCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGtt  <  2:608382/139‑1 (MQ=255)                                                                                                                            tCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAgcggcgag                                                                                          <  2:462325/45‑3 (MQ=255)                                                                                                                            tCCGGTAAGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGgcgcg                                                                                          >  1:462325/1‑45 (MQ=255)                                                                                                                                   |                                                                                                                              TGGTATTAACCCCCCGAGCCTCTACGCGGCTTTTGGCAGTAAAGCTGGGTTATTTAGCCGTGTACTCAATGAATACGTCGGTACGGAAGCTATTCCGCTTGCCGATATTCTTCGTGATGATCGTCCGGTAGGCGAGTGCCTGGTTGAGGTATTAAAAGAAGCGGCGCGCAGATATAGCCAAAACGGCGGCTGCGCTGGCTGTATGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTT  >  NC_000913/4474261‑4474516

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