Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,339,158 A→C 100% 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→C100.0% 188.6 55intergenic (+22/‑4)yhdG/fistRNA‑dihydrouridine synthase B/DNA‑binding protein Fis
Reads supporting (aligned to +/- strand):  new base (43/12):  ref base (0/0):  total (43/12)

AATTTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTCTGACGT  >  REL606/3339125‑3339190
                                 |                                
aaTTTTGCGTAAACAGAAATAAAGAGCTGACAGCAc                                >  1:114367/1‑36 (MQ=255)
aaTTTTGCGTAAACAGAAATAAAGAGCTGACAGCAc                                >  1:2447302/1‑36 (MQ=255)
aaTTTTGCGTAAACAGAAATAAAGAGCTGACAGCAc                                >  1:5374550/1‑36 (MQ=255)
 aTTTTGCGTAAACAGAAATAAAGAGCTGACAGCACt                               >  1:3752705/1‑36 (MQ=255)
 aTTTTGCGTAAACAGAAATAAAGAGCTGACAGCACt                               >  1:2655990/1‑36 (MQ=255)
 aTTTTGAGTAAACAGAAATAAAGAGCTGACAGCACt                               <  1:4306946/36‑1 (MQ=38)
  ttttGCGTAAACAGAAATAAAGAGCTGACAGCACTa                              >  1:1169165/1‑36 (MQ=255)
  ttttGCGTAAACAGAAATAAAGAGCTGACAGCACTa                              >  1:4186366/1‑36 (MQ=255)
  ttttGCGTAAACAGAAATAAAGAGCTGACAGCACTa                              >  1:236301/1‑36 (MQ=255)
   tttTCGTAAACAGAAATAAAGAGCTGACAGCACTAt                             <  1:5345103/36‑1 (MQ=37)
   tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                             >  1:2531544/1‑36 (MQ=255)
    ttGCGTAAACAGAAATAAAGAGCTGACAGCACTATg                            >  1:815413/1‑36 (MQ=255)
     tGCGTAAACAGAGATAAAGAGCTGACAGCACTATGt                           >  1:1351144/1‑36 (MQ=37)
     tGCGTAAACAGAAATAAAGAGCTGACAGCACTATGt                           >  1:3057886/1‑36 (MQ=255)
      gCGTAAACAGAAATAAAGAGCTGACAGCACTATGtt                          >  1:2116925/1‑36 (MQ=255)
       cGTAAACAGAAATAAAGAGCTGACAGCACTATGTTc                         <  1:1139274/36‑1 (MQ=255)
        gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                        >  1:1815978/1‑36 (MQ=255)
        gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                        >  1:4456461/1‑36 (MQ=255)
         tAAACAGAAATAAAGAGCTGACAGCCCTATGTTCGa                       <  1:5077643/36‑1 (MQ=38)
          aaaCAGAAATAAAGAGCTGACAGCACTATGTTCGaa                      >  1:1761901/1‑36 (MQ=255)
           aaCAGCAATAAAGACCTCACAACACTATGTTCGaac                     <  1:2057482/36‑1 (MQ=12)
           aaCAGAAATAAAGAGCTGACAGCACTATGTTCGaac                     <  1:3802167/36‑1 (MQ=255)
            aCAGAAATAAAGAGCTGACAGCACTATGTTCGaaca                    >  1:286259/1‑36 (MQ=255)
             cAGAAATAAAGAGCTGACAGCACTATGTTCGaacaa                   >  1:4502997/1‑36 (MQ=255)
              aGAAATAAAGAGCTGACAGCACTATGTTCGaacaac                  >  1:3566713/1‑36 (MQ=255)
               gAAATAAAGAGCTGACAGCACTATGTTCGAACAAcg                 >  1:4697735/1‑36 (MQ=255)
                aaaTAAAGAGCTGACAGCACTATGTTCGAACAAcgc                >  1:4526779/1‑36 (MQ=255)
                 aaTAAAGAGCTGACAGCACTATGTTCGAACAAcgcg               >  1:745799/1‑36 (MQ=255)
                 aaTAAAGAGCTGACAGCACTATGTGCGAACAAcgcg               >  1:3751905/1‑36 (MQ=38)
                  aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt              >  1:3788013/1‑36 (MQ=255)
                  aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt              <  1:392650/36‑1 (MQ=255)
                  aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt              <  1:907657/36‑1 (MQ=255)
                  aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt              >  1:1905125/1‑36 (MQ=255)
                   tAAAGAGCTGACAGCACTATGTTCGAACAACGCGTa             >  1:4552734/1‑36 (MQ=255)
                    aaaGAGCTGACAGCACTATGTTCGAACAACGCGTaa            >  1:3577563/1‑36 (MQ=255)
                    aaaGAGCTGAAAGCACTATGTTCGAACAACGCGTaa            >  1:175941/1‑36 (MQ=37)
                     aaGAGCTGACAGCACTATGTTCTAACAACGCGTgaa           >  1:5318317/1‑36 (MQ=37)
                     aaGAGCTGACAGCACTATGTTCGAACAACGCGTaaa           >  1:4235389/1‑36 (MQ=255)
                     aaGAGCTGACAGCACTAGGTTCGAACAACGAGCaaa           >  1:972427/1‑36 (MQ=19)
                      agagCTGACAGCACAATGTTCGAACAACGCGTAAAt          >  1:5502795/1‑36 (MQ=38)
                      agagCTGAAAGCACTATGTTCGAACAACGCGTAAAt          >  1:3591516/1‑36 (MQ=38)
                       gagCTGACAGCACTATGTTCGAACAACGCGTAAAtt         >  1:4643270/1‑36 (MQ=255)
                        agCTGACAGCACTATGTTCGAACAACGCGTAAATTc        >  1:2374036/1‑36 (MQ=255)
                        agCTGACAGCACTATGTTCGAACAACGCGTAAATTc        >  1:332685/1‑36 (MQ=255)
                        agCTGACAGCACTATGTTCGAACAACGCGTAAATTc        >  1:4040623/1‑36 (MQ=255)
                         gCTGACAGCACTATGTTCGAACAACGCGTAAATTCt       >  1:2646980/1‑36 (MQ=255)
                         gCTGACAGCACTATGTTCGAACAACGCGTAAATTCt       <  1:3953123/36‑1 (MQ=255)
                           tGACAGCACTATGTTCGAACAACGCGTAAATTCTGa     >  1:2347417/1‑36 (MQ=255)
                           tGACAGCACTATGTTCGAACAACGCGTAAATTCTGa     >  1:1173444/1‑36 (MQ=255)
                            gACAGCACTATGTTCGAACAACGCGTAAATTCTGAc    <  1:2575029/36‑1 (MQ=255)
                            gACAGCACTATGTTCGAACAACGCGTAAATTCTGAc    >  1:5258550/1‑36 (MQ=255)
                            gACAGCACTATGTTCGAACAACGCGTAAATTCTGAc    <  1:3478959/36‑1 (MQ=255)
                            gACAGCACTATGTTCGAACAACGCGTAAATTCTGAc    <  1:5655364/36‑1 (MQ=255)
                             aCAGCACTATGTTCGAACAACGCGTAAATTCTGACg   >  1:5662243/1‑36 (MQ=255)
                              cAGCACTATGTTCGAACAACGCGTAAATTCTGACGt  >  1:2540661/1‑36 (MQ=255)
                                 |                                
AATTTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTCTGACGT  >  REL606/3339125‑3339190

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG