Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,252,059	C→T	P332S (CCG→TCG)	dhaR →	DNA‑binding transcriptional dual regulator DhaR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,252,059	0	C	T	96.8%	79.3 / NA	31	P332S (CCG→TCG)	dhaR	DNA‑binding transcriptional dual regulator DhaR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/11);  minor base G (0/1);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGA  >  NC_000913/1251835‑1252303                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      tATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCACATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCCtttt                                                                                                                                                                                                                             >  2:376068/1‑251 (MQ=255)                           caTCCGCTAAAACACGTAAAAGCAACCTTTGAAACCAAGCACCATTTTATTGATGGGTTTATACCCCTAAAACCGAAAAAAGAACCGCGGGGAACCAGTTTTTTTTTGTTG‑CCCCATCCTGTGAAACAGAGGCGGCAGTTGATGCCCAGTCAATTAGGAAAA‑GTCAGCCATCCC‑TTCGCTCATATGCCACAGG‑CCGATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCtt                                                                                                                                                                                                   <  2:82675/251‑1 (MQ=255)                                       cacGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCt                                                                                                                                                                                       >  1:131601/1‑251 (MQ=255)                                                       ttGAAGGCAAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAAAACCTCAGGGAACCAGCTTTATTTTGTTT‑CGCCACCCTGGGGACCATATGGGTCAGTTGATGACCAGTCAATTAGGAATT‑GTCCGCC‑TACCCTTCGCTCATATGCCGCAGG‑AGGTTTCTAAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAggg                                                                                                                                                                        <  1:273902/250‑1 (MQ=255)                                                         aaaaGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTTGTCCCCAGGCGGCGCGCAGCAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCgg                                                                                                                                                                     >  2:12471/2‑250 (MQ=255)                                                              ccAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGc                                                                                                                                                                 >  2:393588/1‑250 (MQ=255)                                                                                    gTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTACATACTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCAAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAAt                                                                                                                                         >  1:43394/1‑251 (MQ=255)                                                                                      gATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGTGAAGAGGGGGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTc                                                                                                                                       >  2:90749/1‑251 (MQ=255)                                                                                      gATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑ATCAATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCCCTGCTA‑AGTCAGGCAATTc                                                                                                                                       >  1:478473/1‑251 (MQ=255)                                                                                      gATAACCCTTAAACCGAAAATAGAAACGCGGGGAACCAGTTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTAAATAAGGAAAA‑GTCAGCCATCCC‑TTCGCTCATATGCCACGGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTc                                                                                                                                       <  1:163502/251‑1 (MQ=255)                                                                                       atgaccctaaaacagtaaagaaaaaCCCAGGTAACCATTTTTATTTTGTGG‑CCCCATCCGGTGAAACAGAGGGGACTGTAGATGCCCAGTAAATT‑GTAAAATGTCAGCCTTCCC‑TTCCCTCATATGCCCCAGT‑CCGATTCGCATCCCCGCCGCTTGATTCATTTTGGGCGCCAGGCGGCGCGCAGTAGCTTTCCTGGCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCa                                                                                                                                      <  2:285610/229‑1 (MQ=255)                                                                                               tAAACCGATAATAGAAACGCAGGGAACAAGTGTTATTTGTTTG‑CTCCATCCGGTGAAACGGATGCGGCAGTTGATGCCCAGTCAATTAGGAAAA‑GTCAGCCACACC‑TCCGTTCATATGCCCCAGG‑ACGTTTCTCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGaa                                                                                                                               <  1:60132/250‑1 (MQ=255)                                                                                               aaaaCCAATAATAGAAACGCAGGGAACCAGCTTTTTTTTGTCG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGCCCAGCCAATTAGGAAAA‑GTCAGCCATCCC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGaaa                                                                                                                              <  2:235118/250‑1 (MQ=255)                                                                                                              aaCGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGtgctgc                                                                                                                 >  1:64489/1‑249 (MQ=255)                                                                                                                   agggaaccgtatttATTTT‑TTGACTCCAGCGTGTGGAACAGATGCGGCATTTGGTGAACAGTCAAATAGAAAAA‑TCCAGCCATACC‑GTTGATCATATGCCAAAGG‑ACGATGCGAAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCt                                                                                                          <  2:125286/240‑1 (MQ=255)                                                                                                                       aaCCAGCTTTATTTTGTTG‑CTACATCCTGTGGAACAGATGCGGCAGGTGCTGAACAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCCCAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGTTCGCATGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTATGGTGAAGAGGGAGTGGACAAGGCACTGCTACAG‑CAGGCAATTCCAAATGAAAGCGAGCAGGCTGCAgtgtcttct                                                                                                       >  2:283535/1‑242 (MQ=255)                                                                                                                                 aTTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCa                                                                                             >  1:160397/1‑250 (MQ=255)                                                                                                                                   tttGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGCGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTACTGCAGGTCCTTTTATCGCCGTCAAtt                                                                                          >  2:270029/1‑251 (MQ=255)                                                                                                                                                   tgtgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTg                                                                           >  1:473419/1‑251 (MQ=255)                                                                                                                                                    gtgGAACAGATGCGGCAGTGAATGACCAGTCAATTAGAAAAA‑GTCAGCCATCCC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTg                                                                           <  1:55471/250‑1 (MQ=255)                                                                                                                                                        aaCAGATGCGGCAGTTGATGACCAGTCAATTAGAAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGAtgct                                                                      <  1:96530/251‑1 (MQ=255)                                                                                                                                                         aCAGATGCGGCAGTTGATGACCAGTCAATTAGGTAAA‑GTCAGCCATACC‑TTCGCTCATATGCC‑CAGGGACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGCGCTGCAGGTCCTTAATTCGCCGTCATTTGTGAGTTATATGGTGAtgctg                                                                     >  2:246270/1‑251 (MQ=255)                                                                                                                                                                                 tCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATtggtg                                             >  1:65990/1‑251 (MQ=255)                                                                                                                                                                                  cAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATtggtgg                                            >  1:397191/1‑251 (MQ=255)                                                                                                                                                                                                    ccATACC‑TTCTCTCATTGGCCAAAGG‑ACGATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATg                           <  1:196058/251‑1 (MQ=255)                                                                                                                                                                                                      aTACC‑TTCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGa                          >  1:333400/1‑250 (MQ=255)                                                                                                                                                                                                      aTACC‑GTCGCTCAAAAGCCACAGG‑ATGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGaa                         <  1:203142/251‑1 (MQ=255)                                                                                                                                                                                                            ttCGCTCATATGCCACAGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGAGCAAAGAAAATgg                    >  2:348937/1‑251 (MQ=255)                                                                                                                                                                                                               ccctactcgCCACAGT‑ACTATTCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGTACAATGAAAATGGCCg                 <  2:478473/243‑1 (MQ=255)                                                                                                                                                                                                                         cacaGG‑ACGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCgg       >  2:345128/1‑251 (MQ=255)                                                                                                                                                                                                                                aCGATTCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGGCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGACGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAAAGGCAGTCTGAGTCGGCTGGa  >  1:280867/1‑250 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTG‑CTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAA‑GTCAGCCATACC‑TTCGCTCATATGCCACAGG‑ACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTA‑AGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCGCCGTCAATTGTGAGTTATATGGTGATGCTGCGCTGGCGGAAGAATTTATTGGTGGCGATCGCACGGACAATGAAAATGGCCGTCTGAGTCGGCTGGA  >  NC_000913/1251835‑1252303

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