Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,170,490 T→A 20.9% F127L (TTT→TTA ygjR → predicted NAD(P)‑binding dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,170,4900T→A20.9% 3.9 ‑4.5 0.034381ygjRpredicted NAD(P)‑binding dehydrogenase
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (14/12):  total (14/20)
Fisher's exact test strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCGCGGGAAAATCAGGTGGTGCTGTTTGAGGCATTTAAAACCGCCTGCCTGCCGAACTTTCATTTGCTGCG  >  REL606/3170455‑3170525
                                   |                                   
gcgcGGGAAAATCAGGTGGTGCTGTTTGAGGCAttt                                     >  1:5147909/1‑36 (MQ=255)
  gcgGGAAAATCAGGTGGTGCTGTTTGAGGCATTTaa                                   >  1:4238163/1‑36 (MQ=255)
   cgGGAAAATCAGGTGGTGCTGTTTGAGGCATTTaaa                                  <  1:7771577/36‑1 (MQ=255)
   cgGGAAAATCAGGTGGTGCTGTTTGAGGCATTTaaa                                  <  1:7233674/36‑1 (MQ=255)
   cgGGAAAATCAGGTGGTGCTGTTTGAGGCATTTaaa                                  >  1:7630986/1‑36 (MQ=255)
    gggAAAATCAGGTGGTGCTGTTTGAGGCATTTaaaa                                 <  1:6240391/36‑1 (MQ=255)
    gggAAAATCAGGTGGTGCTGTTTGAGGCATTTaaaa                                 >  1:6727521/1‑36 (MQ=255)
    gggAAAATCAGGTGGTGCTGTTTGAGGCATTTaaaa                                 <  1:4151640/36‑1 (MQ=255)
     ggAAAATCAGGTGGTGCTGTTTGAGGCATTTAAAAc                                >  1:5987960/1‑36 (MQ=255)
      gAAAATCAGGTGGTGCTGTTTGAGGCATTTTAAAcc                               >  1:2006242/1‑36 (MQ=255)
      gAAAATCAGGTGGTGCTGTTTGAGGCATTTAAAAcc                               <  1:7439082/36‑1 (MQ=255)
      gAAAATCAGGTGGTGCTGTTTGAGGCATTTAAAAcc                               <  1:379564/36‑1 (MQ=255)
          aTCAGGTGGTGCTGTTTGAGGCATTTAAAACCgcct                           <  1:2139146/36‑1 (MQ=255)
           tCAGGTGGTGCTGTTTGAGGCATTTAAAACCgcctg                          >  1:4267816/1‑36 (MQ=255)
            cAGGTGGTGCTGTTTGAGGCATTTAAAACCgcctgc                         >  1:3457992/1‑36 (MQ=255)
            cAGGTGGTGCTGTTTGAGGCATTTAAAACCgcctgc                         >  1:7794141/1‑36 (MQ=255)
              ggtggtGCTGTTTGAGGCATTTAAAACCgcctgcct                       >  1:4193738/1‑36 (MQ=255)
              ggtggtGCTGTTTGAGGCATTTAAAACCgcctgcct                       >  1:5223581/1‑36 (MQ=255)
               gtggtgCTGTTTGAGGCATTTAAAACCGcctgcctg                      <  1:5933204/36‑1 (MQ=255)
                tggtgCTGTTTGAGGCATTTAAACCCGCctgcctgc                     <  1:3262325/36‑1 (MQ=255)
                tggtgCTGTTTGAGGCATTTAAAACCGCctgcctgc                     >  1:1869776/1‑36 (MQ=255)
                tggtgCTGTTTGAGGAATTAAAAACCGCctgcctgc                     <  1:1527450/36‑1 (MQ=37)
                 ggtgCTGTTTGAGGCATTTAAAACCGCCtgcctgcc                    <  1:7293320/36‑1 (MQ=255)
                 ggtgCTGTTTGAGGCATTTAAAACCGCCtgcctgcc                    <  1:1148860/36‑1 (MQ=255)
                        tttAAGGCTTATAAAACCGCCTTCCTGCCGAACttt             <  1:754465/36‑1 (MQ=19)
                         ttGAGGCATTTAAAACCGCCTGCCTGCCGAACTTTc            >  1:5081958/1‑36 (MQ=255)
                          tGAGGAATTAAAAACCGCCTGCCTGCCGAACTTTCa           <  1:7689621/36‑1 (MQ=37)
                           gAGGCATTTAAAACCGCCTGCCTGCCGAACTTTCAt          >  1:4617102/1‑36 (MQ=255)
                            agccATTAAAAACCGCCTGCCTGCCGAACTTTCAtt         <  1:1206923/33‑1 (MQ=37)
                            aGGCATTAAAACCCGCCTGCCTGCCGAACTTTCAtt         <  1:436644/36‑1 (MQ=38)
                             ggaaTTAAAACCCGCCTGCCTGCCGAACTTTCAttt        <  1:5375251/33‑1 (MQ=37)
                              gCATTTAAAACCGCCTGCCTGCCGAACTTTCATTtg       >  1:2640863/1‑36 (MQ=255)
                                aTTAAAACCCGCCTGCCTGCCGAACTTTCATTtgct     <  1:3801180/36‑1 (MQ=38)
                                aTTAAAACCCGCCTGCCTGCCGAACTTTCATTtgct     <  1:3404893/36‑1 (MQ=38)
                                aTTAAAAACCGCCTGCCTGCCGAACTTTCATTtgct     <  1:597013/36‑1 (MQ=255)
                                   tAAACCCGCCTGCCTGCCGAACTTTCATTTGCTgcg  <  1:5108601/36‑1 (MQ=255)
                                   tAAAACCGCCTGCCTGCCGAACTTTCATTTGCTgcg  >  1:5995971/1‑36 (MQ=255)
                                   |                                   
GCGCGGGAAAATCAGGTGGTGCTGTTTGAGGCATTTAAAACCGCCTGCCTGCCGAACTTTCATTTGCTGCG  >  REL606/3170455‑3170525

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG