Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 716,324 C→T intergenic (‑46/+188) BW25113_RS26090 ← / ← kdpE leader peptide SpeFL/two‑component system response regulator KdpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273716,3240CT88.9% 103.9 / 4.2 36intergenic (‑46/+188)BW25113_RS26090/kdpEleader peptide SpeFL/two‑component system response regulator KdpE
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (9/23);  total (10/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

TCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAA  >  NZ_CP009273/716239‑716407
                                                                                     |                                                                                   
tCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGa                                                                                 >  1:221063/1‑90 (MQ=255)
   gCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAAt                                                                              <  2:438051/90‑1 (MQ=255)
   gCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAAt                                                                              <  1:480302/90‑1 (MQ=255)
     cTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATAc                                                                            <  1:360788/90‑1 (MQ=255)
     cTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATAc                                                                            <  1:374471/90‑1 (MQ=255)
        atatGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGt                                                                         <  1:139019/90‑1 (MQ=255)
        atatGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGt                                                                         <  1:205223/90‑1 (MQ=255)
           tGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGa                                                                      <  2:565994/90‑1 (MQ=255)
           tGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGa                                                                      <  1:406293/90‑1 (MQ=255)
             gggCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGa                                                                      <  1:193591/88‑1 (MQ=255)
             gggCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGa                                                                      >  2:193591/1‑88 (MQ=255)
             gggCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAAc                                                                    >  2:85467/1‑90 (MQ=255)
             gggCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAAc                                                                    >  1:267683/1‑90 (MQ=255)
              ggCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAAt                                                                              >  2:594939/1‑79 (MQ=255)
              ggCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAAt                                                                              <  1:594939/79‑1 (MQ=255)
              ggCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACt                                                                   >  1:548766/1‑90 (MQ=255)
                cATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACtgt                                                                 >  1:48006/1‑90 (MQ=255)
                  tAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGa                                                               <  2:267683/90‑1 (MQ=255)
                  tAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCCGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGa                                                               >  1:7126/1‑90 (MQ=255)
                   aGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGaa                                                              <  2:293068/90‑1 (MQ=255)
                         gCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGAt                                                        <  1:676360/90‑1 (MQ=255)
                         gCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGAt                                                        >  1:137146/1‑90 (MQ=255)
                          cTGTTATTTTCCATTTAATTTACCTTTCCCGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATc                                                       <  2:7126/90‑1 (MQ=255)
                            gTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCaa                                                     >  2:608713/1‑90 (MQ=255)
                                       tttaatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGAt                                          <  2:702193/90‑1 (MQ=255)
                                        ttaatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATc                                         <  2:221063/90‑1 (MQ=255)
                                        ttaatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATc                                         <  2:548766/90‑1 (MQ=255)
                                                  tttCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAatt                               <  1:352670/90‑1 (MQ=255)
                                                  tttCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAatt                               <  2:137146/90‑1 (MQ=255)
                                                  tttCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAatt                               <  1:465335/90‑1 (MQ=255)
                                                       tGATGTGAATGAAGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAatttattt                          <  2:48006/90‑1 (MQ=255)
                                                              aatgaaGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAatgatg                   <  1:298827/90‑1 (MQ=255)
                                                              aatgaaGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAatgatg                   <  2:562585/90‑1 (MQ=255)
                                                              aatgaaGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAatgatg                   <  2:590594/90‑1 (MQ=255)
                                                              aatgaaGGTGAAGATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAatgatg                   <  2:678489/90‑1 (MQ=255)
                                                                          gATATTCTATCTAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTc       <  1:383463/90‑1 (MQ=255)
                                                                               tctatctaGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTaa  <  1:608713/90‑1 (MQ=255)
                                                                               tctatctaGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTaa  <  2:731970/90‑1 (MQ=255)
                                                                                     |                                                                                   
TCCGCCTTATATGGGGCATAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAA  >  NZ_CP009273/716239‑716407

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: