Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 393,086 A→G 15.0% pseudogene (19/345 nt) yajB ← putative glycoprotein; b0404_1

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606393,0860A→G15.0% 4.1 ‑3.9 ‑0.037pseudogene (19/345 nt)yajBputative glycoprotein; b0404_1
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (13/18):  total (19/18)
Fisher's exact test strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.59e-01

CGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAG  >  REL606/393051‑393121
                                   |                                   
cgtaaaGTCCGGCTACAGCTGCGACCAGTGGCGGga                                     <  1:7862558/33‑1 (MQ=37)
cGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGga                                     <  1:1725142/36‑1 (MQ=255)
cGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGga                                     <  1:3491921/36‑1 (MQ=255)
  gAAAGTCCGGCGACAGCTGCGACCAGTGGCGGggaa                                   >  1:4986795/1‑36 (MQ=255)
   aaaGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag                                  <  1:6022126/36‑1 (MQ=255)
    aaGTCCGGCGACAGCTGCGACCAGTGGCGGgtaaga                                 >  1:513164/1‑36 (MQ=255)
     aGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaa                                <  1:4164052/36‑1 (MQ=255)
     aGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaa                                >  1:6290984/1‑36 (MQ=255)
      gTCCGGCGACAGCTGCGACCAGTGGCGGggaagaaa                               >  1:7366278/1‑36 (MQ=255)
      gTCCGGCGACAGCTGCGACCAGTGGCGGggaagaaa                               >  1:7218655/1‑36 (MQ=255)
      gTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaa                               >  1:1390412/1‑36 (MQ=255)
      gTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaa                               <  1:5332403/36‑1 (MQ=255)
       tCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaaa                              >  1:5511123/1‑35 (MQ=255)
       tCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag                              <  1:4604926/36‑1 (MQ=255)
       tCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag                              <  1:5108603/36‑1 (MQ=255)
        ccGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGt                             <  1:1041500/36‑1 (MQ=255)
         cGGCGACAGCTGCGACCAGTGGCGGGgaagaaagag                            >  1:3376985/1‑34 (MQ=37)
         cGGCGACAGCTGCGACCAGTGGCGGGGAAGAAAATg                            >  1:2831257/1‑36 (MQ=37)
          ggCGACAGCTGCGACCAGTGGCGGGGAAGAAAGGGa                           >  1:3343233/1‑36 (MQ=37)
          ggCGACAGCTGCGACCAGTGGCGGGGAAGAAAGAGa                           >  1:5311355/1‑36 (MQ=37)
           gCGACAGCTGCGACCAGTGGCGGGGAAGAAAGGGAt                          >  1:3260400/1‑36 (MQ=37)
           gCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGAt                          <  1:2680947/36‑1 (MQ=255)
           gCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGAt                          <  1:7432576/36‑1 (MQ=255)
            cGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATc                         >  1:6829108/1‑36 (MQ=255)
            cGACAGCTGCGACCAGTGGCGGGAAAGAAAATGGTc                         >  1:5098546/1‑36 (MQ=25)
             gACAGCTGCGACCAGTGGCGGGGAAGAAAGGGATcc                        >  1:4798762/1‑36 (MQ=37)
              aCAGCTGCGACCAGTGGCGGGAAAGAAAGTGATccc                       <  1:3859293/36‑1 (MQ=255)
              aCAGCTGCGACCAGTGGCGGGAAAGAAAGTGATccc                       >  1:296220/1‑36 (MQ=255)
               cAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCca                      <  1:933156/36‑1 (MQ=255)
               cAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCca                      <  1:566375/36‑1 (MQ=255)
                aGCTGCGACCAGTGGGGGGAAAGAAAGTGATCCcac                     <  1:3029580/36‑1 (MQ=255)
                aGCTCCGACCAGTGGCGGGAAAGAAAGTGATCCcac                     <  1:1066417/36‑1 (MQ=255)
                 gCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca                    <  1:1557267/36‑1 (MQ=255)
                   tGCGACCAGTGGCGGGGAAGAAAATGATCCCCCATg                  >  1:2106821/1‑36 (MQ=25)
                     cGACCAGTGGCGGGAAAGAAAGTGATCCCACATGAc                <  1:4231336/36‑1 (MQ=255)
                     cGACCAGTGGCGGGAAAGAAAGTGATCCCACATGAc                <  1:7758527/36‑1 (MQ=255)
                     cGACCAGTGGCGGGAAAGAAAGTGATCCCACATGAc                <  1:7050790/36‑1 (MQ=255)
                      gACCAGTGGCGGGAAAGAAAGTGATCCCACATGACa               >  1:6730802/1‑36 (MQ=255)
                          aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCa           <  1:339929/36‑1 (MQ=255)
                          aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCa           <  1:1318505/36‑1 (MQ=255)
                            tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGc         >  1:7772886/1‑36 (MQ=255)
                            tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGc         >  1:6432073/1‑36 (MQ=255)
                             ggCGGGAAAGAAAGTGATCCCACATGACATCCAGCg        >  1:3946057/1‑36 (MQ=255)
                             ggCGGGAAAGAAAGTGATCCCACATGACATCCAGCg        <  1:1727205/36‑1 (MQ=255)
                              gCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt       <  1:358905/36‑1 (MQ=255)
                               cGGGAAAGAAAGTGATCCCACATGACATCCAGCGTa      >  1:5834509/1‑36 (MQ=255)
                               cGGGAAAGAAAGTGATCCCACATGACATCCAGCGTa      >  1:643453/1‑36 (MQ=255)
                                gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa     >  1:6710595/1‑36 (MQ=255)
                                gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa     >  1:6461996/1‑36 (MQ=255)
                                gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa     <  1:5002028/36‑1 (MQ=255)
                                gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa     >  1:5977256/1‑36 (MQ=255)
                                 ggAAAGAAAGTGATCCCACATGACATCCAGCGTAAt    >  1:1625624/1‑36 (MQ=255)
                                 ggAAAGAAAGTGATCCCACATGACATCCAGCGTAAt    <  1:2778696/36‑1 (MQ=255)
                                   aaagaaagTGATCCCACATGACATCCAGCGTAATAg  >  1:4205106/1‑36 (MQ=255)
                                   aaagaaagTGATCCCACATGACATCCAGCGTAATAg  <  1:2980368/36‑1 (MQ=255)
                                   |                                   
CGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAG  >  REL606/393051‑393121

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