Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,209,196 C→A 18.5% intergenic (‑123/‑89) pepE ← / → yjbC peptidase E/23S rRNA pseudouridine synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,209,1960C→A18.5% 4.7 ‑4.5 0.058intergenic (‑123/‑89)pepE/yjbCpeptidase E/23S rRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  new base (1/12):  ref base (21/24):  total (22/36)
Fisher's exact test strand distribution p-value = 1.11e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCAT  >  REL606/4209161‑4209228
                                   |                                
tGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTc                                  <  1:8337731/36‑1 (MQ=255)
tGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTc                                  <  1:4864862/36‑1 (MQ=255)
tGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTc                                  <  1:698544/36‑1 (MQ=255)
 gCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCa                                 >  1:5627072/1‑36 (MQ=255)
 gCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCa                                 <  1:5724745/36‑1 (MQ=255)
 gCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCa                                 >  1:6144051/1‑36 (MQ=255)
  cgtCTGGAAAGTGATTTGAGCAGGTAACGAATTCaa                                >  1:2336049/3‑36 (MQ=255)
  cTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCaa                                <  1:8125149/36‑1 (MQ=255)
  cTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCaa                                <  1:2753194/36‑1 (MQ=255)
    tCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGc                              >  1:7298891/1‑36 (MQ=255)
      tGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCa                            >  1:3579173/1‑36 (MQ=255)
       ggAAAGTGATTTGAGCAGGTAACTAATTCAAGCCaa                           <  1:3759719/36‑1 (MQ=255)
       ggAAAGTGATTTGAGCAGGTAACGAATTCAAGCCaa                           <  1:5478693/36‑1 (MQ=255)
       ggAAAGTGATTTGAGCAGGTAACGAATTCAAGCCaa                           <  1:2654363/36‑1 (MQ=255)
         aaaGTGATTTGAGCAGGTAACGAATTCAAGCCAATg                         >  1:4960983/1‑36 (MQ=255)
         aaaGTGATTTGAGCAGGTAACGAATTCAAGCCAATg                         >  1:2625425/1‑36 (MQ=255)
         aaaGGGAGTTGAGCAGGTAACGAATTCAAGCCATTg                         >  1:2659977/1‑36 (MQ=37)
          aaGTGATTTGAGCAGGGAACGAATTCAAGCCAATgg                        >  1:338409/1‑36 (MQ=255)
          aaGTGATTTGAGCAGGGAACGAATTAAATCCAATgg                        >  1:2809693/1‑36 (MQ=37)
           aGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGa                       <  1:6341659/36‑1 (MQ=255)
           aGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGa                       <  1:50962/36‑1 (MQ=255)
            gTGATTTGAGCAGGTAACGAATTCAAGCCAATGGaa                      <  1:6269843/36‑1 (MQ=255)
            gTGATTTGAGCAGGTAACGAATTCAAGCCAATGGaa                      <  1:8287507/36‑1 (MQ=255)
            gTGATTTGAGCAGGTAACGAATTCAAGCCAATGGaa                      <  1:1032894/36‑1 (MQ=255)
             tGATTTGAGCAGGTAACGAATTCAAGCCAATGGaaa                     >  1:2711624/1‑36 (MQ=255)
             tGATTTGAGCAGGTAACGAATGCAAGACAATGGaaa                     >  1:624447/1‑36 (MQ=38)
              gATTTGAGCAGGTAACGAATTCAAGCCAATGGAAAt                    <  1:6463321/36‑1 (MQ=255)
              gATTTGAGCAGGTAACGAATTCAAGCCAATGGAAAt                    <  1:1412150/36‑1 (MQ=255)
              gATTTGAGCAGGTAACGAATTCAAGCCAATGGAAAt                    <  1:743695/36‑1 (MQ=255)
              gATTTGAGAAGGTAACGAATTCAAGCCAATGGAAAt                    <  1:2773095/36‑1 (MQ=255)
               aTTTGAGCAGGTAACGAATTCAAGCCAATGGAAATa                   >  1:5009451/1‑36 (MQ=255)
               aTTTGAGCAGGTAACGAATTCAAGCCAATGGAAATa                   <  1:4007363/36‑1 (MQ=255)
                tttGAGCAGGTAACGAATTAAAGCCAATGGAAATaa                  <  1:4166894/36‑1 (MQ=255)
                tttGAGCAGGTAACGAATTAAAGCCAATGGAAATaa                  <  1:5659609/36‑1 (MQ=255)
                tttGAGAATATAACGAATTCACGCCAATGGAAATaa                  <  1:138390/36‑1 (MQ=19)
                tttGAGAAGGTAACTAATTAAAGCCAATGGAAATaa                  <  1:7454277/36‑1 (MQ=25)
                   gAGCAGGTAACGAATTAAAGCCAATGGAAATAACaa               <  1:1282897/36‑1 (MQ=255)
                   gAGCAGGCAACGAATTCATGCCAATGGAAAGAACaa               >  1:5970749/1‑36 (MQ=25)
                   gAGAAGGTAACGAATTCAAGCCAATGGAAATAACaa               <  1:2144834/36‑1 (MQ=255)
                    agaaGGTAACGAATTCAAGCCAATGGAAATAACaaa              <  1:1967189/33‑1 (MQ=255)
                    aGCAGGTAACGAATTCAAGCCAATGGAAATAACaaa              >  1:8086248/1‑36 (MQ=255)
                    aGCAGGTAACGAATTCAAGCCAATGGAAATAACaaa              >  1:1923653/1‑36 (MQ=255)
                    aGCAGGTAACGAATTAAAGCCAATGGAAATAACaaa              <  1:8463600/36‑1 (MQ=255)
                     gCAGGTAACGAATTCAAGCCAATGGAAATAACAAAt             >  1:6325811/1‑36 (MQ=255)
                      aaGGTAACGAATTAAAGCCAATGGAAATAACAAATa            <  1:3324893/35‑1 (MQ=255)
                        ggTAAAGAATTCAAGCCAATGGAAATAACAAATAGt          <  1:2270225/36‑1 (MQ=255)
                         aaaaCGAATTAAAGCCAATGGAAATAACAAATAGtt         <  1:5354375/34‑1 (MQ=38)
                          tAACGAATTCAAGCCAATGGAAATAACAAATAGTTc        >  1:2311065/1‑36 (MQ=255)
                          tAAAGAATTAAAGCCAATGGAAATAACAAATAGTTc        <  1:2796415/36‑1 (MQ=37)
                           aaagAATTCAAGCCAATGGAAATAACAAATAGTTCg       <  1:736073/33‑1 (MQ=255)
                           aaCGAATTCAAGCCAATGGAAATAACAAATAGTTCg       >  1:7244371/1‑36 (MQ=255)
                            acaaaTTAAAGCCAATGGAAATAACAAATAGTTCGc      <  1:1088055/33‑1 (MQ=37)
                            aagAATTCAAGCCAATGGAAATAACAAATAGTTCGc      <  1:7920296/34‑1 (MQ=255)
                            aaaaaTTCAAGCCAATGGAAATAACAAATAGTTCGc      <  1:7092076/33‑1 (MQ=255)
                             agAATTAAACCCAATTGAAATAACAAATAGTTCACt     <  1:1160883/35‑1 (MQ=16)
                              gAATTCAAGCCAATGGAAATAACAAATAGTTCGCTc    <  1:6224742/36‑1 (MQ=255)
                              gAATTCAAGCCAATGGAAATAACAAATAGTTCGCTc    >  1:2538215/1‑36 (MQ=255)
                              gAATTCAAGCCAATGGAAATAACAAATAGTTCGCTc    >  1:8338442/1‑36 (MQ=255)
                               aaTTCAAGCCAATGGAAATAACAAATAGTTCGCTCa   >  1:7651124/1‑36 (MQ=255)
                               aaTTAAAGCCAATGGAAATAACAAATAGTTCGCTCa   <  1:52655/36‑1 (MQ=255)
                                aTTAAAGCCAATGGAAATAACAAATAGTTCGCTCAt  <  1:2013904/36‑1 (MQ=255)
                                   |                                
TGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCAT  >  REL606/4209161‑4209228

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG