Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 751,320 A→T 13.1% N821Y (AAT→TAT)  ybgG → alpha‑mannosidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606751,3200A→T13.1% 2.1 ‑4.1 ‑0.0362461ybgGalpha‑mannosidase
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (12/19):  total (17/19)
Fisher's exact test strand distribution p-value = 1.64e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.65e-01

GAAGACCGACAAGAAGTGATTTTACGGCTGTTTAATCCGGCTGAATCAGCAACCTGTGATGCGACTGTT  >  REL606/751287‑751355
                                 |                                   
gAAGACCGACAAGAAGTGATTTTACGGCTGTTTAAt                                   <  1:6526767/36‑1 (MQ=255)
 aaGACCGACAAGAAGTGATTTTACGGCTGTTTAATc                                  <  1:1761755/36‑1 (MQ=255)
   gACCGACAAGAAGTGATTTTACGGCTGTTTTATTCg                                >  1:6731143/1‑36 (MQ=37)
    aCCGACAAGAAGTGATTTTGCGGCTGTTTAATCCgg                               <  1:2442391/36‑1 (MQ=255)
    aCCGACAAGAAGTGATTTTACGGCTGTGTTATccct                               >  1:6279923/1‑34 (MQ=37)
     ccGACAAGAAGTGATTTTACGGCTGTTTAATCCGGc                              <  1:5464091/36‑1 (MQ=255)
      cGACAAGAAGTGATTTTACGGCTGTTTAATCCGGCt                             <  1:2936737/36‑1 (MQ=255)
       gACAAGAAGTGATTTTACGGCTGTTTTATTCGGCtt                            >  1:824321/1‑35 (MQ=37)
         cAAGAAGTGATTTTACGGCTGTTTAATCCGGCTGaa                          <  1:4827952/36‑1 (MQ=255)
          aagaagTGATTTTACGGCTGTTTAATCCGGCTGAAt                         <  1:7893680/36‑1 (MQ=255)
            gaagTGATTTTACGGCTGTTTAATCCGGCTGAATCa                       <  1:1010725/36‑1 (MQ=255)
             aagTGATTTTACGGCTGTTTTATTCGGCTGTATTAg                      >  1:5709535/1‑36 (MQ=19)
             aagTGATTTTACGGCTGTTTAATCCGGCTGAATCAg                      <  1:6918988/36‑1 (MQ=255)
             aagTGATTTTACGGCTGTTTAATCCGGCTGAATCAg                      <  1:1255330/36‑1 (MQ=255)
               gTGATTTTACGGCTGTTTTATCCTGCTGAATTAGCa                    >  1:2165880/1‑36 (MQ=25)
                tGATTTTACGGCTGTTTAATCCGGCTGAATCAGCaa                   <  1:2506230/36‑1 (MQ=255)
                tGATTTTACGGCTGTTTAATCCGGCTGAATCAGCaa                   <  1:6738425/36‑1 (MQ=255)
                 gATTTTACGGCTGTTTAATCCGGCTGAATCAGCAAc                  <  1:6805249/36‑1 (MQ=255)
                 gATTTTACGGCTGTTTAAACCGGCTGAATCCGCAAc                  >  1:7735771/1‑36 (MQ=25)
                  aTTTTACGGCTGTTTAATCCTGCTTAATTATCaaac                 >  1:2536628/1‑34 (MQ=16)
                  aTTTTACGGCTGTTTAATCCGGCTGAATCAGCAAcc                 <  1:7770385/36‑1 (MQ=255)
                  aTTTTACGGCTGTTTAATCCGGCTGAATCAGCAAcc                 <  1:1922043/36‑1 (MQ=255)
                     ttACGGCTGTTTAATCCTGCTGTATCAACAACCtgt              >  1:2274410/1‑36 (MQ=25)
                      tACGGCTGTTTAATCCGGCTGAATCAGCAACCtgtg             <  1:3918024/36‑1 (MQ=255)
                        cGGCTGTTTAATCCGGCTGAATCAGCAACCTGTGAt           >  1:1027644/1‑36 (MQ=255)
                         ggCTGTTTAATCCGGCTGAATCAGCAACCTGTGATg          >  1:4058016/1‑36 (MQ=255)
                          gCTGTTTAATCCGGCTGAATCAGCAACTTGGGATGc         >  1:5909621/1‑36 (MQ=37)
                          gCTGTTTAATCCGGCTGAATCAGCAACCTGTGATGc         >  1:4980219/1‑36 (MQ=255)
                          gCTGTTTAATCCGGCTGAATCAGCAACCTGTGATGc         >  1:3863340/1‑36 (MQ=255)
                           cTGTTTAATCCGGCTGAATCAGCAACCTGTGATGCg        >  1:4501105/1‑36 (MQ=255)
                            tGTTTAATCCGGCTGAATCAGCAACCTGTGATGCGa       >  1:5325299/1‑36 (MQ=255)
                            tGTTTAATCCGGCTGAATCAGCAACCTGTGATGCGa       >  1:1707974/1‑36 (MQ=255)
                             gTTTAATCCGGCTGAATCAGCAACCTGTGATGCGat      >  1:6478469/1‑35 (MQ=255)
                             gTTTAAGCCCTCTGAGTCATGATCCTGTGATGCGAc      <  1:7753341/36‑1 (MQ=2)
                              tttAATCCGGCTGAATCAGCAACCTGTGATGCGACt     <  1:5525331/36‑1 (MQ=255)
                                tAATCCTGCTGAATCAGCAACCTGTGATGCGACTGt   <  1:5778051/36‑1 (MQ=255)
                                 aaTCCGGCTGAATCAGCAACCTGTGATGCGACTGtt  <  1:4036929/36‑1 (MQ=255)
                                 aaTCCGGCTGAATCAGCAACCTGTGATGCGACTGtt  >  1:2222861/1‑36 (MQ=255)
                                 aaTCCGGCTGAATCAGCAACCTGTGATGCGACTGtt  >  1:1508139/1‑36 (MQ=255)
                                 aaTCCGGCTGAATCAGCAACCTGTGATGCGACTGtt  <  1:948522/36‑1 (MQ=255)
                                 |                                   
GAAGACCGACAAGAAGTGATTTTACGGCTGTTTAATCCGGCTGAATCAGCAACCTGTGATGCGACTGTT  >  REL606/751287‑751355

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