Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,096,955 C→G 33.7% G324G (GGC→GGG baeS → sensory histidine kinase in two‑component regulatory system with BaeR

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,096,9550C→G33.7% 2.5 ‑4.5 0.025972baeSsensory histidine kinase in two‑component regulatory system with BaeR
Reads supporting (aligned to +/- strand):  new base (8/1):  ref base (5/11):  total (13/12)
Fisher's exact test strand distribution p-value = 1.12e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTAGATTTGATCCCACTGCTGGAAGTGGCGGGCGGCGC.ATTTCGCGAACGATTCGCCAGTCGTGGCCTGAA  >  REL606/2096920‑2096990
                                   |                                    
gTAGATTTGATCCCACTGCTGGAAGTGGCGggcggc                                      <  1:7024286/36‑1 (MQ=255)
    aTTTGATCCCACTGCTGGAAGTGGCGGGCGgggcgg                                  >  1:1933099/1‑34 (MQ=37)
      ttGATCCCACTGCTGGAAGTGGCGGGGGGGGC.Cttt                               >  1:1902989/1‑36 (MQ=21)
      ttGATCCCACTGCTGGAAGTGGCGGGCGGCGC.Attt                               <  1:2095204/36‑1 (MQ=255)
       tGATCCCACTGCTGGAAGTGGCGGGCGGCGC.ATTTc                              <  1:5041378/36‑1 (MQ=255)
         aTCCCACTGCTGGAAGTGGCGGGGGGGGC.Ctttttc                            >  1:5102069/1‑33 (MQ=16)
         aTCCCACTGCTGGAAGTGGCGGGCGGGGC.GCTTcgc                            >  1:2877825/1‑36 (MQ=21)
         aTCCCACTGCTGGAAGTGGCGGGCGGCGC.ATTTcgc                            <  1:1350297/36‑1 (MQ=255)
          tCCCACTGCTGGAAGTGGCGGGGGGCGG.GTTTcgcc                           >  1:1397433/1‑35 (MQ=21)
            ccACTGCTGGAAGTGGCGGGGGGGGCGATTTcgcgg                          >  1:3425909/1‑35 (MQ=16)
            ccACTGCTGGAAGTGGCGGGCGGGGG.GTTTcgcgga                         >  1:151485/1‑34 (MQ=18)
             cACTGCTGGAAGTGGCGGGGGGGGC.CTTTTGCGGAc                        >  1:3663884/1‑36 (MQ=12)
             cACTGCTGGAAGTGGCGGGCGGCGC.ATTTCGCGAAc                        <  1:2302859/36‑1 (MQ=255)
              aCTGCTGGAAGTGGCGGGCGGGGC.CTTTCGCGGACg                       >  1:3980381/1‑36 (MQ=21)
              aCTGCTGGAAGTGGCGGGCGGCGC.ATTTCGCGAACg                       <  1:1907363/36‑1 (MQ=255)
                 gctgGAAGTGGCGGCGGGCGG.CTTTTGCGCACGGtt                    >  1:2256384/1‑36 (MQ=0)
                  acggAAGGGGGGGGCGGGGC.ATTGCGCGAACGATTc                   <  1:5651984/34‑1 (MQ=16)
                    ggAAGTGGCGGGCGGCGC.ATTTCGCGAACGATTCGc                 <  1:2628302/36‑1 (MQ=255)
                      aaGTGGCGGGCGGGGC.ATTTAGCGAACGATTCGCCa               >  1:2027924/1‑36 (MQ=37)
                       aGTGGCGGGCGGCGC.ATTTCGCGAACGATTCGCCAg              <  1:401872/36‑1 (MQ=255)
                        gTGGCGGGCGGCGC.ATTTCGCGAACGATTCTCCAGt             <  1:1082485/36‑1 (MQ=255)
                        gTGGCGGGCGGCGC.ATTTCGCGAACGATTCGCCAGt             <  1:6903183/36‑1 (MQ=255)
                          ggcgggcgGCGC.ATTTCGCGAACGATTCGCCAGTCg           >  1:7026906/1‑36 (MQ=255)
                           gcgggcggCGC.ATTTCGCGAACGATTCGCCAGTCGt          <  1:7247158/36‑1 (MQ=255)
                            cgggcggCGC.ATTTCGCGAACGATTCGCCAGTCGTg         >  1:4648953/1‑36 (MQ=255)
                            cgggcggCGC.ATTTCGCGAACGATTCGCCAGTCGTg         <  1:1143297/36‑1 (MQ=255)
                             gggcggCGC.ATTTCGCGAACGATTCGCCAGTCGTgg        >  1:870145/1‑36 (MQ=255)
                                cggcgC.ATTTCGCGAACGATTCGCCAGTCGTGGCCt     >  1:4194551/1‑36 (MQ=255)
                                 ggcgC.ATTTCGCGAACGATTCGCCAGTCGTGGCCTg    <  1:7638778/36‑1 (MQ=255)
                                 ggcgC.ATTTCGCGAACGATTCGCCAGTCGTGGCCTg    >  1:2474595/1‑36 (MQ=255)
                                  gcgc.ATTTCGCGAACGATTCGCCAGTCGTGGCCTGa   >  1:5204078/1‑36 (MQ=255)
                                  gcgc.ATTTCGCAAACGATTCGCCAGTCGTGGCCTGa   <  1:6198386/36‑1 (MQ=255)
                                   cgc.ATTTCGCGAACGATTCGCCAGTCGTGGCCTGaa  <  1:3613610/36‑1 (MQ=255)
                                   |                                    
GTAGATTTGATCCCACTGCTGGAAGTGGCGGGCGGCGC.ATTTCGCGAACGATTCGCCAGTCGTGGCCTGAA  >  REL606/2096920‑2096990

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG