Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,339,158 A→C 54.9% intergenic (+22/‑4) yhdG → / → fis tRNA‑dihydrouridine synthase B/DNA‑binding protein Fis

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,339,1580A→C54.9% 49.7 ‑1.5 ‑0.040intergenic (+22/‑4)yhdG/fistRNA‑dihydrouridine synthase B/DNA‑binding protein Fis
Reads supporting (aligned to +/- strand):  new base (10/12):  ref base (12/6):  total (22/18)
Fisher's exact test strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.18e-01

AAAATTTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTCTGACGTACT  >  REL606/3339123‑3339193
                                   |                                   
aaaaTTTTGCGTAAACAGAAATAAAGAGCTGACAGa                                     >  1:2850322/1‑36 (MQ=255)
   aTTTTGCGTAAACAGAAATAAAGAGCTGACAGAACt                                  >  1:2781006/1‑36 (MQ=255)
   aTTTTGCGTAAACAGAAATAAAGAGCTGACAGAACt                                  >  1:4595771/1‑36 (MQ=255)
    ttttGCGTAAACAGAAATAAAGAGCTGACAGCACTa                                 <  1:779238/36‑1 (MQ=255)
    ttttGCGTAAACAGAAATAAAGAGCTGACAGCACTa                                 >  1:5373430/1‑36 (MQ=255)
     tttGCTTAAACAGAAATAAAGAGCTGACAGAACTAt                                <  1:6218415/36‑1 (MQ=255)
     tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                                <  1:308439/36‑1 (MQ=255)
     tttGCGTAAACAGAAATAAAGAGCTGACAGAACTAt                                >  1:784901/1‑36 (MQ=255)
     tttGCGTAAACAGAAATAAAGAGCTGACAGAACTAt                                <  1:6049872/36‑1 (MQ=255)
      ttGCGTAAACAGAAATAAAGAGCTGACAGAACTATg                               >  1:6684109/1‑36 (MQ=255)
         cGTAAACAGAAATAAAGAGCTGACAGAACTATGTTc                            <  1:4760907/36‑1 (MQ=255)
          gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                           <  1:1974598/36‑1 (MQ=255)
          gTAAACAGAAATAAAGAGCTGACAGAACTATGTTCg                           <  1:5287336/36‑1 (MQ=255)
             aaCAGAAATAAAGAGCTGACAGCACTATGTTCGaac                        >  1:3420122/1‑36 (MQ=255)
             aaCAGAAATAAAGAGCTGACAGCACTATGTTCGaac                        >  1:5806668/1‑36 (MQ=255)
               cAGAAATAAAGAGCTGACAGCACTATGTTCGaacaa                      >  1:4287249/1‑36 (MQ=255)
                 gAAATAAAGAGCTGACAGCACTATGTTCGAACAAcg                    >  1:3431811/1‑36 (MQ=255)
                 gAAATAAAGAGCTGACAGCACTATGTGAGAACAAcg                    >  1:4268125/1‑36 (MQ=37)
                 gAAATAAAGAGCTGACAGAACTATGTTCGAACAAcg                    >  1:3511907/1‑36 (MQ=255)
                  aaaTAAAGAGCTGACAGAACTATGTTCGAACAAcgc                   >  1:6861586/1‑36 (MQ=255)
                   aaTAAAGAGCTGACAGCACTATGTTCGAACAAcgcg                  <  1:875434/36‑1 (MQ=255)
                    aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt                 <  1:3984671/36‑1 (MQ=255)
                    aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt                 >  1:3999363/1‑36 (MQ=255)
                    aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt                 >  1:5309010/1‑36 (MQ=255)
                    aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt                 <  1:3169738/36‑1 (MQ=255)
                      aaaGAGCTGACAGAACTATGTTCGAACAACGCGTaa               >  1:3329014/1‑36 (MQ=255)
                       aaGAGCTGACAGCAATATGTTCGAACAACGCGTaaa              <  1:51076/36‑1 (MQ=37)
                       aaGAGCTGAAAGCATTATTTTCTAACAACGCGTaaa              <  1:4380529/36‑1 (MQ=16)
                        agagCTGACAGAACTATGTTCGAACAACGCGTAAAt             >  1:381336/1‑36 (MQ=255)
                         gagCTGACAGCACTATGTTCGAACAACGCGTAAAtt            <  1:2615538/36‑1 (MQ=255)
                         gagCTGACAGCACTATGTTCGAACAACGCGTAAAtt            <  1:2612656/36‑1 (MQ=255)
                          agCTGACAGCACTATGTTCGAACAACGCGTAAATTc           >  1:4167079/1‑36 (MQ=255)
                           gCTGACAGCACTATGTTCGAACAACGCGTAAATTCt          <  1:7708761/36‑1 (MQ=255)
                           gCTGACAGAACTATGTTCGAACAACGCGTAAATTCt          <  1:3255369/36‑1 (MQ=255)
                           gCTGACAGAACTATGTTCGAACAACGCGTAAATTCt          >  1:3127878/1‑36 (MQ=255)
                             tGACAGCACTATGTTCGAACAACGCGTAAATTCTGa        <  1:1559236/36‑1 (MQ=255)
                             tGACAGAACTATGTTCGAACAACGCGTAAATTCTGa        >  1:7368843/1‑36 (MQ=255)
                              gACAGCACTATGTTCGAACAACGCGTAAATTCTGAc       >  1:29892/1‑36 (MQ=255)
                                cAGAACTATGTTCGAACAACGCGCAAATTCTGACGt     >  1:7623075/1‑36 (MQ=255)
                                  gAACTATGTTCGAACAACGCGTAAATTCTGACGTAc   <  1:4565107/36‑1 (MQ=255)
                                  gAACTATGTTCGAACAACGCGTAAATTCTGACGTAc   >  1:5346561/1‑36 (MQ=255)
                                   aaCTATGTTCGAACAACGCGTAAATTCTGACGTACt  >  1:6003846/1‑36 (MQ=255)
                                   aaCTATGTTCGAACAACGCGTAAATTCTGACGTACt  >  1:6699860/1‑36 (MQ=255)
                                   |                                   
AAAATTTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTCTGACGTACT  >  REL606/3339123‑3339193

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG