Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,158,576 C→T intergenic (+97/‑94) baeR → / → yegP two‑component system response regulator BaeR/YegP family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,158,5760CT86.8% 104.2 / 6.9 39intergenic (+97/‑94)baeR/yegPtwo‑component system response regulator BaeR/YegP family protein
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (12/21);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

GGAAGCCGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAAT  >  NZ_CP009273/2158452‑2158708
                                                                                                                            |                                                                                                                                    
ggAAGCCGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCAt                                                                                                                          >  2:332567/1‑137 (MQ=255)
    gCCGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGaaa                                                                                                                    <  2:392094/139‑1 (MQ=255)
     ccGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGaaaa                                                                                                                   >  1:1204829/1‑139 (MQ=255)
     ccGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGaaaa                                                                                                                   <  2:166435/139‑1 (MQ=255)
        aCGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCg                                                                                                                >  1:1061015/1‑139 (MQ=255)
          gCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCt                                                                                                              <  1:1001951/139‑1 (MQ=255)
               ccGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATAc                                                                                                         >  2:209052/1‑139 (MQ=255)
               ccGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATAc                                                                                                         >  1:56297/1‑139 (MQ=255)
                  cATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACttt                                                                                                      <  2:642639/139‑1 (MQ=255)
                         tAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCgg                                                                                                <  2:802852/138‑1 (MQ=255)
                               tAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACa                                                                                          <  2:1180641/138‑1 (MQ=255)
                                aGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAg                                                                                         <  2:56297/138‑1 (MQ=255)
                                               ttAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGt                                                                         >  2:1248389/1‑139 (MQ=255)
                                                  gCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGtttt                                                                      <  2:1123091/139‑1 (MQ=255)
                                                  gCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGtttt                                                                      <  2:979938/139‑1 (MQ=255)
                                                    cGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCg                                                                    <  2:1024382/139‑1 (MQ=255)
                                                             cGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAAt                                                           <  1:209052/139‑1 (MQ=255)
                                                                 cAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGt                                                                                               <  2:756104/99‑1 (MQ=255)
                                                                 cAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGt                                                                                               >  1:756104/1‑99 (MQ=255)
                                                                         aTCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCa                                                                                       <  2:320088/99‑1 (MQ=255)
                                                                         aTCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCa                                                                                       >  1:320088/1‑99 (MQ=255)
                                                                         aTCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCt                                                                           <  2:920244/111‑1 (MQ=255)
                                                                         aTCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCt                                                                           >  1:920244/1‑111 (MQ=255)
                                                                         aTCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGa                                               <  2:297480/139‑1 (MQ=255)
                                                                                aaCATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGAGACACAATTTCAATCAAGGAGTCATTa                                        <  1:1248389/139‑1 (MQ=255)
                                                                                   aTTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAg                                                 >  1:1016699/1‑127 (MQ=255)
                                                                                   aTTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAg                                                 <  2:1016699/127‑1 (MQ=255)
                                                                                         ccACGCAAGGTCAGCAAACTGACAGTATCTCTCCGACCTGGCAGGCATCAGAAAATCGCTCATACTTTAAGCGGTAAGCAGCGCCTTTTGAGGCTGTTTTCGATACACAATCGCATTCAAAGAGTCACTATGGCtaggt                               <  2:282684/139‑5 (MQ=255)
                                                                                            cGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCt                                                                                                              >  1:1121666/1‑57 (MQ=255)
                                                                                            cGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCt                                                                                                              <  2:1121666/57‑1 (MQ=255)
                                                                                            cGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCtggttggt                            >  1:51937/1‑139 (MQ=255)
                                                                                            cGCATGGTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCtggttggt                            >  2:622517/1‑139 (MQ=255)
                                                                                                 ggTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTggttggtt                           <  2:123334/135‑1 (MQ=255)
                                                                                                 ggTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTggttggtt                           >  1:123334/1‑135 (MQ=255)
                                                                                                  gTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAAc                      <  2:802013/139‑1 (MQ=255)
                                                                                                  gTCAGCAAACTGACATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAAc                      <  1:1084490/139‑1 (MQ=255)
                                                                                                  gTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGaa                       <  2:1061015/138‑1 (MQ=255)
                                                                                                              aCATTATCTCTCCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAagag            <  2:1204829/137‑1 (MQ=255)
                                                                                                                        tcCGTCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAAt  <  2:1196888/137‑1 (MQ=255)
                                                                                                                            |                                                                                                                                    
GGAAGCCGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAAT  >  NZ_CP009273/2158452‑2158708

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