Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,399,853 A→G T286A (ACC→GCC)  abgR → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,399,8530AG93.8% 95.6 / ‑0.3 32T286A (ACC→GCC) abgRLysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (12/18);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

TGGGCTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGC  >  NZ_CP009273/1399740‑1399963
                                                                                                                 |                                                                                                              
tGGGCTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCAt                                                                                                       <  1:137165/123‑1 (MQ=255)
tGGGCTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCAt                                                                                                       >  2:137165/1‑123 (MQ=255)
tGGGCTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTACTCACGCAATTCCg                                                                                       >  2:597820/1‑139 (MQ=255)
tGGGCTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCg                                                                                       >  1:5779/1‑139 (MQ=255)
  ggCTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGc                                                                                     <  1:481272/139‑1 (MQ=255)
    cTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGc                                                                                     >  2:481272/1‑137 (MQ=255)
               tGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTa                                                                        >  1:275883/1‑139 (MQ=255)
                    ggACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGc                                                                   <  1:800818/139‑1 (MQ=255)
                    ggACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGc                                                                   <  1:693318/139‑1 (MQ=255)
                    ggACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGc                                                                   <  2:323843/139‑1 (MQ=255)
                      aCAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCa                                                                  >  1:62179/1‑138 (MQ=255)
                                gTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACt                                                       >  2:128925/1‑139 (MQ=255)
                                        gCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGc                                                                                     <  2:278834/101‑1 (MQ=255)
                                        gCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGc                                                                                     >  1:278834/1‑101 (MQ=255)
                                             ttAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATa                                          <  1:128925/139‑1 (MQ=255)
                                             ttAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATACCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATa                                          <  2:515828/139‑1 (MQ=255)
                                              tAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATaaaaaaa                                            >  2:866640/1‑136 (MQ=255)
                                                                aGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGcc                       <  1:485999/139‑1 (MQ=255)
                                                                   ggCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGa                    <  1:423108/139‑1 (MQ=255)
                                                                     ccTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGAct                  <  1:366301/139‑1 (MQ=255)
                                                                         ttagtTGATTCAGCGGCGTGATAGTCGCCAGACACCCCTGGCCGCGTCATTAATCACGCTATTCCTGCGTGAATGCGGCTATCGGCAAAGTTAAAAATGTATAAAAAAATAGAGTCTGTCGATCTCCGCCAGAcgct                <  2:607571/137‑1 (MQ=255)
                                                                           attTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACa             <  2:252004/138‑1 (MQ=255)
                                                                           attTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACa             <  2:446030/138‑1 (MQ=255)
                                                                                aTTCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTaca        <  1:866640/138‑1 (MQ=255)
                                                                                 ttCAGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTacaca      >  1:312905/1‑139 (MQ=255)
                                                                                   cAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCa                                                                                                 >  1:55119/1‑46 (MQ=255)
                                                                                   cAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCa                                                                                                 <  2:55119/46‑1 (MQ=255)
                                                                                    aGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGAc                              >  1:558922/1‑112 (MQ=255)
                                                                                    aGCGGCGTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGAc                              <  2:558922/112‑1 (MQ=255)
                                                                                     gcggcgTGATAGTCGCCAGACACCACTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACAcagc  <  1:135693/139‑1 (MQ=255)
                                                                                                            cacTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGAc                              >  1:768004/1‑88 (MQ=255)
                                                                                                            cacTGGCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGAc                              <  2:768004/88‑1 (MQ=255)
                                                                                                                 aCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATaga                                        <  1:796394/73‑1 (MQ=255)
                                                                                                                 aCCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATaga                                        >  2:796394/1‑73 (MQ=255)
                                                                                                                 |                                                                                                              
TGGGCTGCGATCCCTTGCACGGACAGGGGCTGGTGATGTTGCCGGTTAGCGAAATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGC  >  NZ_CP009273/1399740‑1399963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: