Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,388,893 C→A 29.3% L45F (TTG→TTT ycjI ← murein peptide amidase A

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,388,8930C→A29.3% 10.1 ‑4.2 ‑0.036135ycjImurein peptide amidase A
Reads supporting (aligned to +/- strand):  new base (0/11):  ref base (11/14):  total (11/25)
Fisher's exact test strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.76e-01

GCTGCAGGGGCAGGAAACCAGATTAACGGCGCGCCCAATAATGAACGTCCGTAATGTTCTGTTCCGGGCGG  >  REL606/1388858‑1388928
                                   |                                   
gctgCAGGGGCAGGAAACCAGATTAACGGCGCGccc                                     <  1:7147856/36‑1 (MQ=255)
gctgCAGGGGCAGGAAACCAGATTAACGGCGCGccc                                     >  1:2171095/1‑36 (MQ=255)
  tgCAGGGGCAGGAAACCAGATTAACGGCGCGCCCaa                                   >  1:2425348/1‑36 (MQ=255)
  tgCAGGGGCAGGAAACCAGATTAACGGCGCGCCCaa                                   >  1:270741/1‑36 (MQ=255)
  tgCAGGGGCAGGAAACCAGATTAACGGCGCGCCCaa                                   >  1:5212195/1‑36 (MQ=255)
   ggaGGGTCAGGAAACCAGATTAACGGCGCGCCCaat                                  <  1:2829105/34‑1 (MQ=38)
   gCAGGGGCAGGAAACCAGATTAACGGCGCGCCCaat                                  <  1:6394781/36‑1 (MQ=255)
    cAGGGGCAGGAAACCAGATTAACGGCGCGCCCaata                                 <  1:6907401/36‑1 (MQ=255)
     aGGGGCAGGAAACCAGATTAACGGCGCGCCCaataa                                <  1:6358106/36‑1 (MQ=255)
         gCAGGAAACCAGATTAACGGCGCGCCCAATAATGaa                            <  1:210866/36‑1 (MQ=255)
          cAGGAAACCAGATTAACGGCGCGCCCAATAATGAAc                           <  1:6005131/36‑1 (MQ=255)
           aGGAAACCAGATTAACGGCGCGCCCAATAATGAACg                          >  1:890087/1‑36 (MQ=255)
             gAAACCAGATTAACGGCGCGCCCAATAATGAACGTc                        <  1:5179841/36‑1 (MQ=255)
              aaaCCAGATTAACGGCGCGCCCAATAATGAACGTcc                       >  1:7801385/1‑36 (MQ=255)
               aaCCAGATTAACGGCGCGCCCAATAATGAACGTCCg                      >  1:4577738/1‑36 (MQ=255)
                aCCAGATTAACGGCGCGCCAAATAATGAACGTCCGt                     <  1:5086196/36‑1 (MQ=255)
                  cAGAATAACGGCGCGAAAAATAAGGAACGTCCGTaa                   <  1:1797307/36‑1 (MQ=16)
                    gATTAACGGCGCGCTAAATAATGAACGTCCGTAAtg                 <  1:5225503/36‑1 (MQ=37)
                      ttAACTGCGCGCCAAATAATGAACGTCCGTAAtgtt               <  1:369955/36‑1 (MQ=37)
                       tAACGGCTCGCCCGATAATGACCGTCCGTAAtgttc              <  1:5177558/36‑1 (MQ=25)
                       tAACGGCGCGCCCAATAATGAACGTCCGTAAttttc              >  1:3843773/1‑36 (MQ=255)
                       tAACGGCGCGCCCAATAATGAACGTCCGTAAtgttc              <  1:759180/36‑1 (MQ=255)
                       tAACGGCGCGCCAAATAATGAACGTCCGTAAtgttc              <  1:7360900/36‑1 (MQ=255)
                        aaCGGCGCGCCCAATAATGAACGTCCGTAAtgttct             <  1:7185675/36‑1 (MQ=255)
                         aCGGCGCGCCAAATAATGAACGTCCGTAAtgttctg            <  1:7466776/36‑1 (MQ=255)
                          cGGCGCGCCCAATAATGAACGTCCGTAAtgttctgt           >  1:7695171/1‑36 (MQ=255)
                          cGGCGCGCCCAATAATGAACGTCCGTAAtgttctgt           >  1:656916/1‑36 (MQ=255)
                          cGGCGCGCCAAATAATGAACGTCCGTAAtgttctgt           <  1:4994817/36‑1 (MQ=255)
                          cGGCGCGCCAAATAATGAACGTCCGTAAtgttctgt           <  1:3986551/36‑1 (MQ=255)
                          cGGCGCGCAAAATAAGGAACGTCCGTAAtgttcttt           <  1:3838998/36‑3 (MQ=21)
                           ggCGCGCCAAATAATGAACGTCCGTAAtgttctgtt          <  1:4241931/36‑1 (MQ=255)
                           ggCGCGCACAATAATGAACGTCCGTAAtgttctgtt          <  1:1909397/36‑1 (MQ=255)
                             cgcgcCCAATAATGAGCGTCCGTAATGTTCTGTTcc        <  1:694439/36‑1 (MQ=255)
                              gcgcCCAATAATGAACGTCCGTAATGTTCTGTTCCg       >  1:1861578/1‑36 (MQ=255)
                              gcgcCAAATAATGAACGTCCGTAATGTTCTGTTCCg       <  1:6283944/36‑1 (MQ=255)
                               cgcCCAATAATGACCGTCCGTAATGTTCTGTTCCgg      <  1:7698741/36‑1 (MQ=255)
                               cgcCCAATAATGAACGTCCGTAATGTTCTGTTCCgg      <  1:1412505/36‑1 (MQ=255)
                                gcCCAATAATGAACGTCCGTAATGTTCTGTTCCggg     >  1:5384725/1‑36 (MQ=255)
                                 cccAATAATGAACGTCCGTAATGTTCTGTTCCGGGc    <  1:3206947/36‑1 (MQ=255)
                                   cAATAATGAACGTCCGTAATGTTCTGTTCCGGGCgg  >  1:6444656/1‑36 (MQ=255)
                                   |                                   
GCTGCAGGGGCAGGAAACCAGATTAACGGCGCGCCCAATAATGAACGTCCGTAATGTTCTGTTCCGGGCGG  >  REL606/1388858‑1388928

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG