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%	198.7 / NA	61	A325T (GCT→ACT)	dhaR	dhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/30);  total (31/30)

TTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCA  >  NC_000913/1251906‑1252166                                                                                                                                     |                                                                                                                                 tttATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            <  2:1406179/139‑1 (MQ=255) tttATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            <  2:330933/139‑1 (MQ=255)        aTGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            <  1:4207/132‑1 (MQ=255)        aTGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            >  2:4207/1‑132 (MQ=255)         tGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAg                                                                                                                    >  1:657969/1‑139 (MQ=255)                   aCCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            >  1:528531/1‑121 (MQ=255)                   aCCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            <  2:528531/121‑1 (MQ=255)                          aaCCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            <  2:33494/114‑1 (MQ=255)                          aaCCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                            >  1:33494/1‑114 (MQ=255)                          aaCCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACccg                                                                                                   <  1:1320882/139‑1 (MQ=255)                               ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCCGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCt                                                                                              >  2:985286/1‑139 (MQ=255)                               ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCat                                                                                                                             <  2:540233/108‑1 (MQ=255)                               ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCat                                                                                                                             >  1:540233/1‑108 (MQ=255)                                         cGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCAttt                                                                                    >  1:1472193/1‑139 (MQ=255)                                          gCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCAtttt                                                                                   >  1:1570805/1‑139 (MQ=255)                                              ggAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTc                                                                               >  1:55066/1‑139 (MQ=255)                                               gAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCg                                                                              <  2:657969/139‑1 (MQ=255)                                                 aCCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACccg                                                                                                   <  1:1562954/116‑1 (MQ=255)                                                 aCCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACccg                                                                                                   >  2:1562954/1‑116 (MQ=255)                                                    aGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAgg                                                                         <  1:639379/139‑1 (MQ=255)                                                                ttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                               >  1:249637/1‑105 (MQ=255)                                                                ttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                               <  2:249637/105‑1 (MQ=255)                                                                ttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGgcggcg                                                                    >  2:479678/1‑132 (MQ=255)                                                                ttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGgcggcg                                                                    <  1:479678/132‑1 (MQ=255)                                                                ttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTa                                                             >  1:1346500/1‑139 (MQ=255)                                                                        tCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcct                                                     >  1:1625326/1‑139 (MQ=255)                                                                             tgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtc                                                       <  2:1116111/132‑1 (MQ=255)                                                                             tgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtc                                                       >  1:1116111/1‑132 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCg                                                                                                           >  1:424977/1‑75 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCg                                                                                                           <  2:424977/75‑1 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                               >  1:895894/1‑87 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                               <  2:895894/87‑1 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtc                                                  <  1:236331/132‑1 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtc                                                  >  2:236331/1‑132 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCttt                                           >  1:1015366/1‑139 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCttt                                           <  2:1472193/139‑1 (MQ=255)                                                                                     aTGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTtgtg                                        >  2:136080/1‑139 (MQ=255)                                                                                      tGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTgg                                       >  2:1558221/1‑139 (MQ=255)                                                                                            aGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaaga                                 <  2:1625326/139‑1 (MQ=255)                                                                                             gTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGgcggcg                                                                    <  2:130442/103‑1 (MQ=255)                                                                                             gTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGgcggcg                                                                    >  1:130442/1‑103 (MQ=255)                                                                                                gatgaCCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCt                                                          <  1:537997/110‑1 (MQ=255)                                                                                                gatgaCCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCt                                                          >  2:537997/1‑110 (MQ=255)                                                                                                gatgaCCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGc                             <  1:1558221/139‑1 (MQ=255)                                                                                                     ccAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAggcgg                                                                      <  2:248027/93‑1 (MQ=255)                                                                                                     ccAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAggcgg                                                                      >  1:248027/1‑93 (MQ=255)                                                                                                        gTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                            <  1:880220/100‑1 (MQ=255)                                                                                                        gTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                            >  2:880220/1‑100 (MQ=255)                                                                                                        gTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCaa                     >  2:959012/1‑139 (MQ=255)                                                                                                          cAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg                   <  2:1346500/139‑1 (MQ=255)                                                                                                           aaTTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGc                  <  2:55066/139‑1 (MQ=255)                                                                                                               aGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTg              <  1:959012/139‑1 (MQ=255)                                                                                                                   aaaGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaag                                  <  2:1606850/115‑1 (MQ=255)                                                                                                                   aaaGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaag                                  >  1:1606850/1‑115 (MQ=255)                                                                                                                       tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaag                                  <  1:368669/111‑1 (MQ=255)                                                                                                                       tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaag                                  >  2:368669/1‑111 (MQ=255)                                                                                                                       tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGc                  <  2:61476/127‑1 (MQ=255)                                                                                                                       tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGc                  >  1:61476/1‑127 (MQ=255)                                                                                                                        cAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtc                                                  >  2:950453/1‑94 (MQ=255)                                                                                                                        cAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtc                                                  <  1:950453/94‑1 (MQ=255)                                                                                                                           ccATACCTTCACTCATATGCCACAGGACGATCAGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCa  >  2:276887/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                 TTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCA  >  NC_000913/1251906‑1252166

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