Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,495,906 A→T 20.4% L127Q (CTA→CAA)  eutP ← conserved protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,495,9060A→T20.4% 7.3 ‑3.5 0.032380eutPconserved protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (13/12):  total (20/12)
Fisher's exact test strand distribution p-value = 2.88e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCAGTTGCTGCACGCTTTGCGGGTCGTGGCTGTTTAGCTCGAAAATCGGCTCTTCGAAACCCGTCTCAAGC  >  REL606/2495871‑2495941
                                   |                                   
ccaGTTGCTGCACGCTTTGCGGGTCGTGGCTGTTTa                                     <  1:5089411/36‑1 (MQ=255)
ccaGTTGCTGCACGCTTTGCGGGTCGTGGCTGTTTa                                     <  1:1833502/36‑1 (MQ=255)
   gTTGCTGCACGCTTTGCGGGTCGTGGCTGGTTTGCt                                  >  1:2693898/1‑36 (MQ=38)
    ttGCTGCACGCTTTGCGGGTCGTGGCTGTTTAGCTc                                 <  1:3916639/36‑1 (MQ=255)
       ctgcACGCTTTGCGGGTCGTGGCTGTTTTTCTCGaa                              >  1:4136267/1‑36 (MQ=37)
       ctgcACGCTTTGCGGGTCGTGGCTGTTTAGCTCGaa                              <  1:5356613/36‑1 (MQ=255)
         gcACGCTTTGCGGGTCGTGGCTGTTTTGCCCGgaaa                            >  1:187079/1‑36 (MQ=25)
           aCGCTTTGCGGGTCGTGGCTGTTTTGTTCGAAAATc                          >  1:5023997/1‑36 (MQ=37)
           aCGCTTTGCGGGTCGTGGCTGTTTTGATCGAAAAtt                          >  1:349575/1‑35 (MQ=37)
           aCGCTTTGCGGGTCGTGGCTGTTTAGCCCGGAAAta                          >  1:5619154/1‑35 (MQ=37)
            cGCTTTGCGGGTGGAGGCTGTTTTGCTTGAAAAACg                         >  1:1303975/1‑36 (MQ=16)
            cGCTTTGCGGGTCGTGGCTGTTTTGGTCGAAAAACg                         >  1:5254073/1‑36 (MQ=21)
            cGCTTTGCGGGTCGTGGCTGTTTAGCTCGAAAATCg                         >  1:5680049/1‑36 (MQ=255)
             gCTTTGCGGGTCGTGGCTGTTTAGGTCGAAAATCgg                        >  1:2923110/1‑36 (MQ=255)
               tttGCGGGTCGTGGCTGTTTAGCTCGAAAATCGGct                      <  1:587816/36‑1 (MQ=255)
               tttGCGGGTCGTGGCTGTTTAGCTCGAAAATCGGct                      <  1:916367/36‑1 (MQ=255)
                 tGCGGGTCGTGGCTGTTTAGCTCGAAAATCGGctct                    >  1:1580367/1‑36 (MQ=255)
                   cGGGTCGTGGCTGTTTAGCTCGAAAATCGGCTCTTc                  <  1:2304721/36‑1 (MQ=255)
                   cGGGTCGTGGCTGTTTAGCTCGAAAATCGGCTCTTc                  <  1:323207/36‑1 (MQ=255)
                     ggTCGTGGCTGTTTAGCTCGAAAATCGGCTCTTCGa                <  1:869997/36‑1 (MQ=255)
                     ggTCGTGGCTGTTTAGCTCGAAAAACGGCTCtttgg                >  1:5040290/1‑33 (MQ=37)
                        cGTGGCTGTTTAGCTCGAAAATCGGCTCTTCGAAAc             <  1:1019504/36‑1 (MQ=255)
                           ggCTGTTTAGCTCGAAAATCGGCTCTTCGAAACCCg          >  1:2837568/1‑36 (MQ=255)
                           ggCTGTTTAGCTCGAAAATCGGCTCTTCGAAACCCg          >  1:1017653/1‑36 (MQ=255)
                           ggCTGTTTAGCGCGAAAATCGGCTCTTCGAAACCCg          >  1:5500469/1‑36 (MQ=255)
                            gCTGTTTAGCTCGAAAATCGGCTCTTTGAAACCCGt         >  1:4784205/1‑36 (MQ=255)
                            gCTGTTTAGCTCGAAAATCGGCTCTTCGAAACCCGt         >  1:661457/1‑36 (MQ=255)
                             cTGTTTAGCTCGAAAATCGGCTCTTCGAAACCCGtc        >  1:3346095/1‑36 (MQ=255)
                              tGTTTAGCTCGAAAATCGGCTCTTCGAAACCCGtct       >  1:5890396/1‑36 (MQ=255)
                              tGTTTAGCTCGAAAATCGGCTCTTCGAAACCCGtct       <  1:1584093/36‑1 (MQ=255)
                               gTTTAGCTCGAAAATCGGCTCTTCGAAACCCGtctc      <  1:3220182/36‑1 (MQ=255)
                                tttAGCTCGAAAATCGGCTCTTCGAAACCCGTCTCa     <  1:2929261/36‑1 (MQ=255)
                                tttAGCTCGAAAATCGGCTCTTCGAAACCCGTCTCa     <  1:2679506/36‑1 (MQ=255)
                                tttAGCTCGAAAATCGGCTCTTCGAAACCCGTCTCa     >  1:253766/1‑36 (MQ=255)
                                ttaaGCTCGAAAATCGGCTCTTCGAAACCCGTCTCa     <  1:3071196/33‑1 (MQ=255)
                                   aGCTCGAAAATCGGCTCTTCGAAACCCGTCTCAAGc  <  1:440480/36‑1 (MQ=255)
                                   aGCTCGAAAATCGGCTCTTCGAAACCCGTCTCAAGc  <  1:2432601/36‑1 (MQ=255)
                                   aGCTCGAAAATCGGCTCTTCGAAACCCGTCTCAAGc  <  1:1860328/36‑1 (MQ=255)
                                   |                                   
CCAGTTGCTGCACGCTTTGCGGGTCGTGGCTGTTTAGCTCGAAAATCGGCTCTTCGAAACCCGTCTCAAGC  >  REL606/2495871‑2495941

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 28 ≤  ATCG  < 37 ≤  ATCG