Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,252,038	G→A	A325T (GCT→ACT)	dhaR →	dhaKLM operon transcription activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,252,038	0	G	A	100.0%	106.3 / NA	35	A325T (GCT→ACT)	dhaR	dhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

CGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGC  >  NC_000913/1251916‑1252165                                                                                                                           |                                                                                                                                cGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGa                                                                                                                 <  1:107544/139‑1 (MQ=255)                aaCCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACccg                                                                                                  <  1:22644/139‑1 (MQ=255)                     ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCt                                                                                             >  1:1057970/1‑139 (MQ=255)                        ataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGa                                                                                          <  2:1057970/139‑1 (MQ=255)                               cGCAGTGAACCAGCTTTATTTTGTTTCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCAttt                                                                                   <  1:534528/139‑1 (MQ=255)                                                    tgttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAg                                                              <  2:103263/139‑1 (MQ=255)                                                       tgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                           <  2:1038586/139‑1 (MQ=255)                                                           ccATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCttt                                                       >  1:163674/1‑139 (MQ=255)                                                                      aaCAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCTATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCt                                            <  1:878141/139‑1 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                           >  1:636343/1‑58 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                           <  2:636343/58‑1 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                              >  2:1105979/1‑135 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                              <  1:1105979/135‑1 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCttt                                          >  2:578823/1‑139 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCttt                                          >  2:123505/1‑139 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCttt                                          >  1:927322/1‑139 (MQ=255)                                                                        cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAACTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCttt                                          >  2:779119/1‑139 (MQ=255)                                                                           aTGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTtgtg                                       >  1:410840/1‑139 (MQ=255)                                                                            tGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTgg                                      <  2:163674/139‑1 (MQ=255)                                                                                gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAgg                                                                                                                  <  2:877466/59‑1 (MQ=255)                                                                                gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAgg                                                                                                                  >  1:877466/1‑59 (MQ=255)                                                                                gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaa                                  <  1:578823/139‑1 (MQ=255)                                                                                              gTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCa                     <  1:976553/138‑1 (MQ=255)                                                                                              gTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCa                     >  2:976553/1‑138 (MQ=255)                                                                                                cAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg                  <  2:927322/139‑1 (MQ=255)                                                                                                 aaTTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGgcgcgc                                                                >  1:709277/1‑92 (MQ=255)                                                                                                 aaTTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGgcgcgc                                                                <  2:709277/92‑1 (MQ=255)                                                                                                         aaaGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTaa         >  1:594311/1‑139 (MQ=255)                                                                                                            gTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTc      <  1:779119/139‑1 (MQ=255)                                                                                                             tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTgg                                      >  1:515720/1‑106 (MQ=255)                                                                                                             tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTgg                                      <  2:515720/106‑1 (MQ=255)                                                                                                                gCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGc  >  1:533548/1‑139 (MQ=255)                                                                                                                gCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTCGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGc  <  1:39045/139‑1 (MQ=255)                                                                                                                 ccATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                           <  1:714507/81‑1 (MQ=255)                                                                                                                 ccATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGAAGTAg                                                           >  2:714507/1‑81 (MQ=255)                                                                                                                           |                                                                                                                                CGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGC  >  NC_000913/1251916‑1252165

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