Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,170,704 G→T 14.2% G199C (GGT→TGT)  ygjR → predicted NAD(P)‑binding dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,170,7040G→T14.2% 2.3 ‑3.5 ‑0.044595ygjRpredicted NAD(P)‑binding dehydrogenase
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (19/18):  total (26/18)
Fisher's exact test strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.91e-01

CTTTTACTGCCTGGCGTCGGCGGTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTA  >  REL606/3170670‑3170739
                                  |                                   
cTTTTACTGCCTGGCGTCGGCGGTAGCGTTATTTgg                                    <  1:4821301/36‑1 (MQ=255)
cTTTTACTGCCTGGCGTCGGCGGTAGCGTTATTTgg                                    <  1:1656732/36‑1 (MQ=255)
 ttttACTGCCTGGCGTCGGCGGTAGGGTTATTTgtg                                   >  1:5207486/1‑34 (MQ=37)
 ttttACTGCCTGGCGTCGGCGGTAGCGTTTTTTTGt                                   >  1:6187278/1‑36 (MQ=37)
   ttACTGCCTGGCGTCGGCGGTAGCGTTATTTGGTGa                                 <  1:6692690/36‑1 (MQ=255)
   ttACTGCCTGGCGTCGGCGGTAGCGTTATTTGGTGa                                 <  1:3349982/36‑1 (MQ=255)
     aCTGCCTGGCGTCGGCGGTAGCGTTATTTGGTGAGc                               <  1:4766336/36‑1 (MQ=255)
      cTGCCTGGCGTCGGCGGTAGCGTTATTTGGTGAGcc                              <  1:1537707/36‑1 (MQ=255)
      cTGCCTGGCGTCGGCGGTAGCGGTATTTTGGGAggc                              >  1:2841529/1‑34 (MQ=21)
       tGCCTGGCGTCGGCGGTAGCGTTATTTGGTGAGCCg                             <  1:6970527/36‑1 (MQ=255)
       tGCCTGGCGTCGGCGGTAGCGTTATTTGGTGAGCCg                             <  1:3324254/36‑1 (MQ=255)
        gCCTGGCGTCGGCGGTAGCGTTTTTTTGTGAGGCGa                            >  1:5471730/1‑36 (MQ=25)
         ccTGGCGTCGGCGGTAGCGTTATTTTGTGAGCCGaa                           >  1:3443948/1‑36 (MQ=255)
         ccTGGCGTCGGCGGTAGCGTTATTTGGTGAGCCGaa                           <  1:631354/36‑1 (MQ=255)
         ccTGGCGTCGGCGGGAGCGTTTTTTGGTGAGCCGaa                           >  1:2080361/1‑36 (MQ=37)
           tGGCGTCGGCGGTAGTGGTATTTTGTGTGACGagaa                         >  1:2411768/1‑36 (MQ=9)
           tGGCGTCGGCGGTAGCGTTATTTGGTGAGCCGaaaa                         >  1:2920998/1‑36 (MQ=255)
           tGGCGTCGGCGGTAGCGTTATTTGGTGAGCCGaaaa                         <  1:1624072/36‑1 (MQ=255)
            ggCGTCGGCGGTAGCGTTTTTTGGTGAGCCGCAAAg                        >  1:3668973/1‑36 (MQ=38)
            ggCGTCGGCGGTAGCGTTATTTGGTGAGCCGAAAAg                        >  1:4523432/1‑36 (MQ=255)
             gCGTCGGCGGTAGCGTTTTTTTGTGAGCCGAAAAGc                       >  1:1940524/1‑36 (MQ=37)
              cGTCGGCGGTAGCGTTATTTGGTGAGCCGAAAAGCg                      <  1:6885855/36‑1 (MQ=255)
                tCGGCGGTAGCGTTTTTTTGTGAGCCGGAAAGCggg                    >  1:4716550/1‑34 (MQ=21)
                  ggcggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCa                  >  1:3308993/1‑36 (MQ=255)
                  ggcggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCa                  >  1:4227609/1‑36 (MQ=255)
                   gcggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAg                 <  1:5686607/36‑1 (MQ=255)
                    cggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAgg                >  1:6357656/1‑36 (MQ=255)
                    cggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAgg                <  1:5930514/36‑1 (MQ=255)
                    cggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAgg                >  1:4436812/1‑36 (MQ=255)
                    cggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAgg                <  1:4811229/36‑1 (MQ=255)
                     ggTAGCGTTATTTGGTGAGGCGAAAAGCGTGCAGGc               <  1:3070572/36‑1 (MQ=255)
                     ggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGc               >  1:7762342/1‑36 (MQ=255)
                     ggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGc               >  1:5905881/1‑36 (MQ=255)
                     ggTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGc               >  1:1354324/1‑36 (MQ=255)
                     ggTAGCGTTATTTGGTGAACCTAAAAGCGTGCAGGc               <  1:370808/36‑1 (MQ=37)
                        aGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAAc            >  1:1079104/1‑36 (MQ=255)
                         gCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAAcc           <  1:4963247/36‑1 (MQ=255)
                         gCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAAcc           <  1:5603354/36‑1 (MQ=255)
                         gCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAAcc           <  1:4864112/36‑1 (MQ=255)
                         gCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAAcc           >  1:4251090/1‑36 (MQ=255)
                         gCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAAcc           <  1:3478502/36‑1 (MQ=255)
                          cGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCg          >  1:3870259/1‑36 (MQ=255)
                          cGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCg          >  1:1137967/1‑36 (MQ=255)
                           gTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGc         >  1:5809844/1‑36 (MQ=255)
                              aTTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAg      >  1:6720237/1‑36 (MQ=255)
                               tttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGt     >  1:4712572/1‑36 (MQ=255)
                               tttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGt     >  1:3665772/1‑36 (MQ=255)
                                  ggTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTa  >  1:4130301/1‑36 (MQ=255)
                                  ggTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTa  >  1:1518832/1‑36 (MQ=255)
                                  |                                   
CTTTTACTGCCTGGCGTCGGCGGTAGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTA  >  REL606/3170670‑3170739

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG