Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,481,490 T→A 15.2% L460I (TTA→ATA)  ydcR → fused predicted DNA‑binding transcriptional regulator/predicted amino transferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,481,4900T→A15.2% 2.5 ‑3.5 ‑0.0441378ydcRfused predicted DNA‑binding transcriptional regulator/predicted amino transferase
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (19/17):  total (27/17)
Fisher's exact test strand distribution p-value = 3.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.64e-01

AGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGT  >  REL606/1481455‑1481525
                                   |                                   
aGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAAt                                     <  1:7084006/36‑1 (MQ=255)
aGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAAt                                     <  1:2010574/36‑1 (MQ=255)
aGTGGGGAGAGCGTGAAGAACAGGCGGGTAAAAAAt                                     >  1:2006576/1‑36 (MQ=25)
 gTGGGGAGAGCGTGAAGACCAGGCGGTAAAACAAtt                                    <  1:5875668/36‑1 (MQ=255)
 gTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAAtt                                    >  1:4539645/1‑36 (MQ=255)
  tGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTa                                   <  1:209974/36‑1 (MQ=255)
    gggAGAGCGTGAAGAACAGGCGGTTAAACCATTaag                                 >  1:3309344/1‑34 (MQ=25)
    gggAGAGCGTGAAGAACAGGCGGTAAAACAATTAgg                                 <  1:3438988/36‑1 (MQ=255)
     ggAAAGCGTGAAGAACAGGCGGTAAAACAATTAGGc                                <  1:3057227/36‑1 (MQ=255)
       agagCGTGAAGAACAGGCGGTAAAACAATTAGGCaa                              <  1:7040415/36‑1 (MQ=255)
       agagCGTGAAGAACAGGCGGGTAAAAAAATAGGGaa                              >  1:6812636/1‑36 (MQ=12)
        gagCGTGAAGAACAGGCGGTAAAACAATTAGGCaaa                             >  1:1107791/1‑36 (MQ=255)
        gagCGTGAAGAACAGGCGGGAAAAAAATTAGGGaaa                             >  1:1512376/1‑36 (MQ=25)
         agCGTGAAGAACAGGCGGTAAAACAATTAGGCAAAc                            <  1:6233869/36‑1 (MQ=255)
         agCGTGAAGAACAGGCGGTAAAACAATTAGGCAAAc                            <  1:7443278/36‑1 (MQ=255)
          gCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACt                           <  1:4437960/36‑1 (MQ=255)
          gCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACt                           <  1:1517586/36‑1 (MQ=255)
           cGTGAAGAACAGGCGGTAAAACAAATAGGCAAACtt                          >  1:5508568/1‑36 (MQ=255)
            gTGAAGAACAGGCGGTAAAAAAAATGTGCAAACTTa                         >  1:3670897/1‑36 (MQ=19)
             tGAAGAACAGGCGGTAAAACAATTAGGCAAACTTAt                        <  1:5364899/36‑1 (MQ=255)
             tGAAGAACAGGCGGTAAAAAAATTATGCAAAATTaa                        >  1:639632/1‑35 (MQ=21)
             tGAAGAACAGGCGGTAAAAAAATTAGGCAAACtttt                        >  1:1204670/1‑34 (MQ=37)
             tGAAGAACAGGCGGTAAAAAAAATAGGCAAAATTAt                        >  1:7026021/1‑36 (MQ=25)
              gaagaaCAGGCGGTAAAAAAAATAGGGAAAATTAtt                       >  1:3729303/1‑36 (MQ=16)
              gaagaaCAGGCGGAAAAACAAGTATGCAAAATTAtt                       >  1:612886/1‑36 (MQ=16)
               aagaaCAGGCGGCAAAAAAAATAAGCAAAAATATTc                      >  1:6815482/1‑36 (MQ=9)
                   aCAGGCGGTAAAAAAAATAAACAAAACTATTTaaaa                  >  1:1784452/1‑34 (MQ=0)
                    cAGGCGGTAAAACAATTAGGCAAACTTATTCAAGaa                 >  1:1218556/1‑36 (MQ=255)
                      ggCGGTAAAACAATTAGGCAAACTTATTCAAGAACg               >  1:7502476/1‑36 (MQ=255)
                       gCGGTAAAACAATTAGGCAAACTTATTCAAGAACgg              >  1:6319758/1‑36 (MQ=255)
                        cGGTAAAACAATTAGGCAAACTTATTCAAGAACGGc             >  1:5413736/1‑36 (MQ=255)
                        cGGTAAAACAATTAGGCAAACTTATTCAAGAACGGc             <  1:4711045/36‑1 (MQ=255)
                         ggTAAAACAATTAGGCAAACTTATTCAAGAACGGCt            >  1:7921237/1‑36 (MQ=255)
                         ggTAAAACAATTAGGCAAACTTATTCAAGAACGGCt            >  1:2714910/1‑36 (MQ=255)
                         ggTAAAACAATTAGGCAAACTTATTCAAGAACGGCt            <  1:842629/36‑1 (MQ=255)
                          gTAAAACAATTAGGCAAACTTATTCAAGAACGGCTg           >  1:5367560/1‑36 (MQ=255)
                          gTAAAACAATTAGGCAAACTTATTCAAGAACGGCTg           <  1:4774269/36‑1 (MQ=255)
                           tAAAACAATTAGGCAAACTTATTCAAGAACGGCTGt          <  1:7486842/36‑1 (MQ=255)
                           tAAAACAATTAGGCAAACTTATTCAAGAACGGCTGt          >  1:3220948/1‑36 (MQ=255)
                            aaaaCAATTAGGCAAACTTATTCAAGAACGGCTGTa         >  1:2945647/1‑36 (MQ=255)
                              aaCAATTAGGCAAACTTATTCAAGAACGGCTGTAAt       <  1:982094/36‑1 (MQ=255)
                                 aaTTAGGCAAACTTATTCAAGAACGGCTGTAATaac    >  1:7340084/1‑34 (MQ=255)
                                 aaTTAGGCAAACTTATTCAAGAACGGCTGTAATAGc    >  1:3631854/1‑36 (MQ=255)
                                 aaTTAGGCAAACTTATTCAAGAACGGCTGTAATAGc    >  1:3318260/1‑36 (MQ=255)
                                 aaTTAGGCAAACTTATTCAAGAACGGCTGTAATAGc    <  1:1699378/36‑1 (MQ=255)
                                 aaTTAGGCAAACTTATTCAAGAACGGCTGTAATAGc    >  1:7997982/1‑36 (MQ=255)
                                 aaTTAGGCAAACTTATTCAAGAACGGCTGTAATAGc    <  1:2043726/36‑1 (MQ=255)
                                 aaTTAGGCAAACTTATTCAAGAACGGCTGTAATAGc    >  1:2088094/1‑36 (MQ=255)
                                 aaTTAGGAAAACTTATTCAAGAACGGCTGTAATAGc    <  1:3836972/36‑1 (MQ=255)
                                   ttAGGCAAACTTATTCAAGAACGGCTGTGATAGCGt  >  1:5022677/1‑36 (MQ=255)
                                   ttAGGCAAACTTATTCAAGAACGGCTGTAATAGCGt  <  1:1402100/36‑1 (MQ=255)
                                   |                                   
AGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGT  >  REL606/1481455‑1481525

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG