Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,801,076 C→A 9.4% I99I (ATC→ATA ydjQ → endonuclease of nucleotide excision repair

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,801,0760C→A 9.4% 2.4 ‑1.5 ‑0.632297ydjQendonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (16/13):  total (16/16)
Fisher's exact test strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.33e-01

GAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTT  >  REL606/1801041‑1801108
                                   |                                
gAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATc                                  <  1:5053310/36‑1 (MQ=255)
gAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATc                                  <  1:1440317/36‑1 (MQ=255)
  aaTCGGCGCTCTGCTCCTTGAAGCGCGATTAATCaa                                >  1:1329092/1‑36 (MQ=255)
   aTCGGCGCTCTGCTCCTTGAAGCGCGATTAATCaaa                               <  1:3047967/36‑1 (MQ=255)
   aTCGGCGCTCTGCTCCTTGAAGCGCGATTAATCaaa                               <  1:226770/36‑1 (MQ=255)
    tCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAg                              >  1:7129778/1‑36 (MQ=255)
    tCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAg                              <  1:2850162/36‑1 (MQ=255)
    tCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAg                              >  1:8468426/1‑36 (MQ=255)
      ggCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGaa                            <  1:6773808/36‑1 (MQ=255)
       gcactctGCTCCTTGAAGCGCGATTAATCAAAGaac                           <  1:7421611/33‑1 (MQ=255)
        cgcTCTGCTCCTTGAAGCGCGATTAATCAAAGaaca                          <  1:155685/36‑1 (MQ=255)
        cgcTCTGCTCCTTGAAGCGCGATTAATCAAAGaaca                          <  1:469627/36‑1 (MQ=255)
         gcTCTGCTCCTTGAAGCGCGATTAATCAAAGaacaa                         >  1:1962311/1‑36 (MQ=255)
         gcTCTGCTCCTTGAAGCGCGATTAATCAAAGaacaa                         <  1:5767891/36‑1 (MQ=255)
          ctctGCTCCTTGAAGCGCGATTAATCAAAGaacaac                        <  1:4767756/36‑1 (MQ=255)
           tctGCTCCTTGAAGCGCGATTAATCAAAGAacaaca                       <  1:6086392/36‑1 (MQ=255)
           tctGCTCCTTGAAGCGCGATTAATCAAAGAacaaca                       <  1:5655560/36‑1 (MQ=255)
              gCTCCTTGAAGCGCGATTAATCAAAGAACAACAGcc                    <  1:8271946/36‑1 (MQ=255)
               cTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCg                   >  1:6107584/1‑36 (MQ=255)
                  cTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTg                >  1:71081/1‑36 (MQ=255)
                    tGAAGCGCGATTAATCAAAGAACAACAGCCGCTGtt              >  1:1036330/1‑36 (MQ=255)
                      aaGCGCGATTAATAAAAGAACAACAGCCGCTGTTta            <  1:3133764/36‑1 (MQ=255)
                      aaGCGAAATAAATAAAAAAACAACAGCCGCTGTTta            <  1:4808251/36‑1 (MQ=16)
                       aGCGCGATTAATCAAAGAACAACAGCCGCTGTTtaa           >  1:7274559/1‑36 (MQ=255)
                          gaaaTTAATAAAAAAAAAANAGCCGCTGTTtaataa        <  1:1373410/33‑1 (MQ=16)
                           cgATTAATCAAAGAACAACAGCCGCTGTTTAATaaa       >  1:384746/1‑36 (MQ=255)
                            gATTAATCAAAGAACAACTGCCGCTGTTTGATAAAc      >  1:2559617/1‑36 (MQ=37)
                            gATTAATCAAAGAACAACAGCCGCTGTTTAATAAAc      >  1:4877109/1‑36 (MQ=255)
                            gATTAATCAAAGAACAACAGCCGCTGTTTAATAAAc      >  1:5976528/1‑36 (MQ=255)
                            gATTAATCAAAGAACAACAGCCGCTGTTTAATAAAc      >  1:6849544/1‑36 (MQ=255)
                            gATTAATCAAAGAACAACAGCCGCTGTTTAATAAAc      >  1:7263629/1‑36 (MQ=255)
                              ttAATCAAAGAACAACAGCAGCTGTTTAATAAACgg    <  1:4119439/36‑1 (MQ=255)
                               tAATCAAAGAACAACAGCCGATGTTTAATAAACGGt   >  1:6853547/1‑36 (MQ=255)
                                aaTCAAAGAACAACAGCCGCTGTTTAATAAACGGtt  >  1:3811276/1‑36 (MQ=255)
                                   |                                
GAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTT  >  REL606/1801041‑1801108

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG