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

TTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCA  >  NC_000913/1251906‑1252162                                                                                                                                     |                                                                                                                             tttATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                        <  2:798549/139‑1 (MQ=255)                  aaCCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCGGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCg                                                                                                       <  1:613704/139‑1 (MQ=255)                        tAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAAcc                                                                                                 >  2:948520/1‑139 (MQ=255)                               ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCt                                                                                          >  2:917421/1‑139 (MQ=255)                               ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCt                                                                                          >  1:676753/1‑139 (MQ=255)                               ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCAGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCt                                                                                          >  1:74815/1‑139 (MQ=255)                                            aGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTgg                                                                             <  1:948520/139‑1 (MQ=255)                                              ggAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTc                                                                           >  1:1059594/1‑139 (MQ=255)                                               gAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCg                                                                          >  2:1526973/1‑139 (MQ=255)                                                      cTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAggcg                                                                   <  2:165427/139‑1 (MQ=255)                                                       tttATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAggcgg                                                                  <  1:1526973/139‑1 (MQ=255)                                                                        tCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcct                                                 >  1:1241096/1‑139 (MQ=255)                                                                             tgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcct                                            >  1:969557/1‑139 (MQ=255)                                                                                aaCAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCt                                         <  2:676753/139‑1 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                        >  1:1410349/1‑58 (MQ=255)                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCata                                                                                                                        <  2:1410349/58‑1 (MQ=255)                                                                                     aTGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTtgtg                                    <  2:1241096/139‑1 (MQ=255)                                                                                     aTGCGGAAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAgg                                                                     >  1:1598848/1‑106 (MQ=255)                                                                                     aTGCGGAAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAgg                                                                     <  2:1598848/106‑1 (MQ=255)                                                                                       gCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCGAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGga                                  >  1:925178/1‑139 (MQ=255)                                                                                       gCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGga                                  >  2:1516298/1‑139 (MQ=255)                                                                                          gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                           >  1:468117/1‑127 (MQ=255)                                                                                          gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                           <  2:468117/127‑1 (MQ=255)                                                                                                 atgaCCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                           >  2:1208896/1‑72 (MQ=255)                                                                                                 atgaCCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                           <  1:1208896/72‑1 (MQ=255)                                                                                                  tgaCCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGt                       <  2:1532743/139‑1 (MQ=255)                                                                                                       aGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                           >  2:505340/1‑66 (MQ=255)                                                                                                       aGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                           <  1:505340/66‑1 (MQ=255)                                                                                                        gTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCaa                 >  2:1587173/1‑139 (MQ=255)                                                                                                         tCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAg                <  2:1059594/139‑1 (MQ=255)                                                                                                          cAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg               <  1:1587173/139‑1 (MQ=255)                                                                                                          cAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg               <  1:917421/139‑1 (MQ=255)                                                                                                                       tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                        <  2:1431240/85‑1 (MQ=255)                                                                                                                       tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                        >  1:1431240/1‑85 (MQ=255)                                                                                                                       tCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCa  >  1:1087858/1‑139 (MQ=255)                                                                                                                           ccATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                           >  2:1535186/1‑94 (MQ=255)                                                                                                                           ccATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                           <  1:1535186/94‑1 (MQ=255)                                                                                                                              tACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtc                                                   >  2:334845/1‑83 (MQ=255)                                                                                                                              tACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtc                                                   <  1:334845/83‑1 (MQ=255)                                                                                                                                     |                                                                                                                             TTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCA  >  NC_000913/1251906‑1252162

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: ‑