Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,612,333 A→G 22.0% intergenic (‑82/+240) ydfQ ← / ← ydfR predicted lysozyme/predicted protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,612,3330A→G22.0% 3.9 ‑4.1 ‑0.026intergenic (‑82/+240)ydfQ/ydfRpredicted lysozyme/predicted protein
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (8/12):  total (14/12)
Fisher's exact test strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.12e-01

TTACCTCATCAGTGATGGGCACTGGCGGGAGGAGGATGTCAATCTGATAAACACAGAGGCGACTA  >  REL606/1612301‑1612365
                                |                                
ttACCTCATCAGTGATGGGCACTGGCGggaggagga                               >  1:1294142/1‑36 (MQ=255)
 tACCTCATCAGTGCTGGGCACTGGCGGGTGGGGGAt                              >  1:4270512/1‑36 (MQ=25)
 tACCTCATCAGTGATGGGCACTGGCGGGAGGAGGAt                              <  1:3993684/36‑1 (MQ=255)
  aCCTCATCAGTGATGGGCACTGGCGGGAGGGGGATg                             >  1:3575105/1‑36 (MQ=255)
      catcaGTGATGGGCACTGGCGGGGGGGGGGTGTcta                         >  1:754611/1‑34 (MQ=25)
         caGTGATGGGCACTGGCGGGGGGGGGGTGTGAAttt                      >  1:1366251/1‑34 (MQ=16)
          aGTGATGGGCACTGGCGGGGGGGGGATGTCAATCTg                     >  1:1567340/1‑36 (MQ=38)
          aGTGATGGGCACTGGCGGGGGGAGGATGTCAATCTg                     >  1:4126933/1‑36 (MQ=255)
          aGTGATGGGCACTGGCGGGAGGAGGGTGTCAATCTg                     >  1:5613082/1‑36 (MQ=255)
             gATGGGCACTGGCGGGGGGGGGGTGTCAATCTGATa                  >  1:4456052/1‑36 (MQ=37)
             gATGGGCACTGGCGGGAGGAGGATGTCAATCTGATa                  <  1:928500/36‑1 (MQ=255)
             gATGGGCACTGGCGGGAGGAGGATGTCAATCTGATa                  <  1:2723613/36‑1 (MQ=255)
              aTGGGCACTGGCGGGGGGAGGGGGTTCATCTTATaa                 >  1:4059812/1‑36 (MQ=12)
              aTGGGCACTGGCGGGGGGAGGAAGTGAATCTGGTaa                 >  1:1527122/1‑36 (MQ=19)
                gggCACTGGCGGGAGGAGGATGTCAATCTGATAAac               <  1:838157/36‑1 (MQ=255)
                  gCACTGGCGGGAGGAGGATGTCAATCTGATAAacac             <  1:1929892/36‑1 (MQ=255)
                  gCACTGGCGGGAGGAGGATGTCAATCTGATAAacac             <  1:5332652/36‑1 (MQ=255)
                   cACTGGCGGGAGGAGGATGTCAATCTGATAAAcaca            <  1:137425/36‑1 (MQ=255)
                    aCTGGCGGGAGGAGGATGTCAATCTGATAAACACag           <  1:4274428/36‑1 (MQ=255)
                     cTGGCGGGAGGAGGATGTCAATCTGATAAACACaga          <  1:5783744/36‑1 (MQ=255)
                     cTGGCGGGAGGAGGATGTCAATCTGATAAACACaga          >  1:479025/1‑36 (MQ=255)
                      tGGCGGGAGGAGGATGTCAATCTGATAAACACagag         <  1:4093181/36‑1 (MQ=255)
                       ggCGGGAGGAGGATGTCAATCTGATAAACACAGAgg        >  1:717373/1‑36 (MQ=255)
                        gCGGGAGGAGGATGTCAATCTGATAAACACAGAGGc       >  1:3551714/1‑36 (MQ=255)
                         cGGGAGGAGGATTTCAATCTGATAAACACAGAGGCg      <  1:2610979/36‑1 (MQ=255)
                          gggAGGAGGATGTCAATCTGATAAACACAGAGGCGa     >  1:1319612/1‑36 (MQ=255)
                          gggAGGAGGATGTCAATCTGATAAACACAGAGGCGa     <  1:1304031/36‑1 (MQ=255)
                           ggaggaggaTGTCAATCTGATAAACATAGAGGCGAc    >  1:2429479/1‑36 (MQ=255)
                           ggaggaggaTGTCAATCTGATAAACACAGAGGCGAc    <  1:986748/36‑1 (MQ=255)
                            gaggaggaTGTCAATCTGATAAACACAGAGGCGACt   >  1:961530/1‑36 (MQ=255)
                             aggaggaTGTCAATCTGATAACCACAGAGGCGACTa  <  1:1731215/36‑1 (MQ=255)
                             aggaggaTGTCAATCTGATAAACACAGAGGCGACTa  >  1:352520/1‑36 (MQ=255)
                                |                                
TTACCTCATCAGTGATGGGCACTGGCGGGAGGAGGATGTCAATCTGATAAACACAGAGGCGACTA  >  REL606/1612301‑1612365

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG