Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,455,519 C→G 8.7% A52P (GCC→CCC)  pyrB ← aspartate carbamoyltransferase catalytic subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,455,5190C→G 8.7% 2.5 ‑3.1 ‑0.046154pyrBaspartate carbamoyltransferase catalytic subunit
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (18/24):  total (22/24)
Fisher's exact test strand distribution p-value = 4.48e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

AAGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCA  >  REL606/4455486‑4455554
                                 |                                   
aaGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCtt                                   <  1:673197/36‑1 (MQ=255)
aaGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCtt                                   <  1:1912625/36‑1 (MQ=255)
aaGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCtt                                   <  1:4501376/36‑1 (MQ=255)
  gTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCg                                 <  1:491658/36‑1 (MQ=255)
   tttCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCga                                <  1:4479624/36‑1 (MQ=255)
   tttCGAAAGAGAGGCGGGGATGGGGAGGGGGttttt                                >  1:1944199/1‑33 (MQ=12)
    ttCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCgaa                               <  1:2214276/36‑1 (MQ=255)
    ttCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCgaa                               <  1:4970626/36‑1 (MQ=255)
      cGAAAGAGAGGCGGGGACGGGGAGAGGGTTCgaagg                             >  1:7376469/1‑35 (MQ=25)
        aaaGAGAGGCGGGTACGGGTAGAGGGTTTGGAGaaa                           >  1:95277/1‑36 (MQ=21)
         aaGGGAGGCGGGTACGGGTAGAGGCTTCGAAGAAAc                          <  1:3036549/36‑1 (MQ=255)
          agagagGCGGGGACGGGGTGAGGGTTCTAATaaaaa                         >  1:671802/1‑34 (MQ=9)
             gagGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCt                      <  1:6542413/36‑1 (MQ=255)
             gagGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCt                      >  1:1043128/1‑36 (MQ=255)
              agGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTg                     <  1:6848041/36‑1 (MQ=255)
              agGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTg                     >  1:2239164/1‑36 (MQ=255)
               ggCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTgg                    <  1:4712545/36‑1 (MQ=255)
                gCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGc                   >  1:5841777/1‑36 (MQ=255)
                  gggtacgggtaGAGGCTTCGAAGAAACAGCTGGCaa                 >  1:4179212/1‑36 (MQ=255)
                   ggtacgggtaGAGGCTTCGAAGAAACAGCTGGCAAt                <  1:2699764/36‑1 (MQ=255)
                    gtaggggtaGAGGCTTCGAAGAAACAGCTGGCAATg               <  1:3481348/36‑1 (MQ=255)
                     tacgggtaGAGGCTTCGACGAAACAGCTGGCAATGa              >  1:7075864/1‑36 (MQ=255)
                      acgggtaGAGGCTTCGAAGAAACAGCTGGCAATGAc             >  1:7861990/1‑36 (MQ=255)
                      acgggtaGAGGCTTCGAAGAAACAGCTGGCAATGAc             <  1:1065084/36‑1 (MQ=255)
                      acgggtaGAGGCTTCGAAGAAACAGCTGGCAATGAc             >  1:3336114/1‑36 (MQ=255)
                       cgggtaGAGGCTTCGAAGAACCAGCTGTCAATTACt            <  1:4863373/36‑1 (MQ=25)
                       cgggtaGAGGCTTCGAAGAAACAGCTGGCAATGACt            <  1:4754548/36‑1 (MQ=255)
                       cgggtaGAGGCTTCGAAGAAACAGCTGGCAATGACt            <  1:371816/36‑1 (MQ=255)
                       cgggtaGAGGCTTCGAAGAAACAGCTGGCAATGACt            <  1:7060784/36‑1 (MQ=255)
                        gggtaGAGGCTTCTAAGAAACAGCTGGCAATGACtt           <  1:1205307/36‑1 (MQ=255)
                        gggtaGAGGCTTCGAAGAAACAGCTGGCAATGACtt           >  1:4942812/1‑36 (MQ=255)
                        gggtaGAGGCTTCGAAGAAACAGCTGGCAATGACtt           >  1:4982794/1‑36 (MQ=255)
                        gggtaGAGGCTTCGAAGAAACAGCTGGCAATGACtt           >  1:1999253/1‑36 (MQ=255)
                         ggtaGAGGCTTCGAAGAAACAGCTGGCAATGACttt          <  1:5807701/36‑1 (MQ=255)
                         ggtaGAGGCTTCGAAGAAACAGCTGGCAATGACttt          <  1:1518799/36‑1 (MQ=255)
                          gtaGAGGCTTCGAAGAAACAGCTGGCAATGACTTtg         >  1:4598937/1‑36 (MQ=255)
                          gtaGAGGCTTCGAAGAAACAGCTGGCAATGACTTtg         >  1:7115009/1‑36 (MQ=255)
                          gtaGAGGCTTCGAAGAAACAGCTGGCAATGACTTtg         >  1:3279919/1‑36 (MQ=255)
                          gtaGAGGCTTCGAAGAAACAGCTGGCAATGACTTtg         >  1:3979152/1‑36 (MQ=255)
                           taGCGGCTTCGAAGAAACAGCTGGCAATGACTTtgt        <  1:3915210/36‑1 (MQ=255)
                           taGAGGCTTCGAAGAAACAGCTGGCAATGACTTtgt        <  1:467449/36‑1 (MQ=255)
                            agagGCTTCGAAGAAACAGCTGGCTATGACTTtgtg       >  1:7555637/1‑36 (MQ=255)
                              agGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCt     >  1:5740468/1‑36 (MQ=255)
                              agGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCt     >  1:3658421/1‑36 (MQ=255)
                               ggCTTCGAAGAACCAGCTGGCAATGACTTTGTGCtt    <  1:6456973/36‑1 (MQ=255)
                               ggCTTCGAAGAAACAGCTGGCAATGACTTTGTGCtt    >  1:4928442/1‑36 (MQ=255)
                                 cTTCTAAAAAACAGCTGGCACTGCCTTTGTGCTTCa  <  1:723599/36‑1 (MQ=16)
                                 cTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCa  <  1:6176255/36‑1 (MQ=255)
                                 cTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCa  >  1:3154317/1‑36 (MQ=255)
                                 |                                   
AAGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCA  >  REL606/4455486‑4455554

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