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_0009131,252,0380GA96.5% 82.8 / ‑6.2 29A325T (GCT→ACT) dhaRdhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/16);  minor base T (1/0);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGC  >  NC_000913/1251905‑1252165
                                                                                                                                     |                                                                                                                               
gTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCat                                                                                                                            <  1:554632/139‑1 (MQ=255)
           cggtcgtAACCCTTAAACCGATAATAGAAACGCACGGAATGAGCTTTATTTTCTTGCTCCAGCCTGTGCCCCAGATGCGGCAGGTGTTGACCAGTCAACTGGGAAAACTGACGCATACCTTAACTCATATGCCACTGGa                                                                                                                 >  1:703915/7‑139 (MQ=255)
           cGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGa                                                                                                                 <  1:339396/139‑1 (MQ=255)
           cGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGa                                                                                                                 >  1:630699/1‑139 (MQ=255)
           cGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGa                                                                                                                 >  2:420795/1‑139 (MQ=255)
             gTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACg                                                                                                               <  1:566277/139‑1 (MQ=255)
                              gATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                              <  1:420795/139‑1 (MQ=255)
                              gATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                              <  1:709799/139‑1 (MQ=255)
                               ataataGCAACGCAGGGAACCAGCTATATATCGTTGCTCCATCCTGTGGACCAGATGCGGCAGTTGATGGCCAGTCCATTAGAAAAACTCAGCCATACCTTCTCTCATATGCCGCCGGACGAACAACACACcctccgat                                                                                             >  1:210219/1‑137 (MQ=255)
                                         cGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGa                                                                                                                 >  2:84755/1‑109 (MQ=255)
                                         cGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGa                                                                                                                 <  1:84755/109‑1 (MQ=255)
                                                 aCCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGcc                                                                           >  2:536531/1‑139 (MQ=255)
                                                         tATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGgcggcg                                                                   <  2:768546/139‑1 (MQ=255)
                                                                 tgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTAAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                           <  1:679059/139‑1 (MQ=255)
                                                                        tCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcct                                                    >  2:423529/1‑139 (MQ=255)
                                                                             tgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                              >  2:96315/1‑92 (MQ=255)
                                                                             tgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCcgccgc                                                                                              <  1:96315/92‑1 (MQ=255)
                                                                             tgGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcct                                               >  1:753106/1‑139 (MQ=255)
                                                                               gAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGa                                                                                                              >  2:129132/1‑74 (MQ=255)
                                                                               gAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGa                                                                                                              <  1:129132/74‑1 (MQ=255)
                                                                                  cAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCttt                                          >  2:320242/1‑139 (MQ=255)
                                                                                          gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                           >  2:594695/1‑114 (MQ=255)
                                                                                          gCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                           <  1:594695/114‑1 (MQ=255)
                                                                                                      cAGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGc                      <  2:585003/139‑1 (MQ=255)
                                                                                                       aGTCAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCa                     <  1:790202/139‑1 (MQ=255)
                                                                                                          cAATTAGGAAAAGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg                  <  2:680860/139‑1 (MQ=255)
                                                                                                                   aaaGTCAGCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTaa         >  2:137328/1‑139 (MQ=255)
                                                                                                                          gCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGc  <  2:630699/139‑1 (MQ=255)
                                                                                                                          gCCATACCTTCACTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGc  <  2:753106/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
GTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGC  >  NC_000913/1251905‑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: