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

AATTTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTCTGACGT  >  REL606/3339125‑3339190
                                 |                                
aaTTTTGCGTAAACAGAAATAAAGAGCTGACAGCAc                                >  1:3424855/1‑36 (MQ=255)
 aTTTTGCGTAAACAGAAATAAAGAGCTGACAGCACt                               <  2:1505929/36‑1 (MQ=255)
  ttttGCGTAAACAGAAAGAAAGAGCTGACAGCACTa                              <  2:3001750/36‑1 (MQ=37)
   tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                             <  1:3733886/36‑1 (MQ=255)
   tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                             >  2:1363365/1‑36 (MQ=255)
   tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                             >  1:89192/1‑36 (MQ=255)
   tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                             <  1:768788/36‑1 (MQ=255)
    ttGCGTAAACAGAAATAAAGAGCTGACAGCACTATg                            >  1:2050495/1‑36 (MQ=255)
     tGCGTAAACAGAAATAAAGAGCTGACAGCACTATGt                           <  1:339432/36‑1 (MQ=255)
      gCGTAAACAGAAATAAAGAGCTGACAGCACTATGtt                          <  2:2351481/36‑1 (MQ=255)
      gCGTAAACAGAAATAAAGAGCTGACAGCACTATGtt                          <  2:1788028/36‑1 (MQ=255)
      acGTAAACAGAAATAAAGAGCTGACAGCACTATGtt                          <  1:919511/35‑1 (MQ=255)
       cGTAAACAGAAATAAAGAGCTGACAGCACTATGTTc                         <  2:1849278/36‑1 (MQ=255)
       cGTAAACAGAAATAAAGAGCTGACAGCACCATGtcc                         >  1:366208/1‑34 (MQ=37)
        gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                        <  2:133926/36‑1 (MQ=255)
        gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                        >  1:2447237/1‑36 (MQ=255)
        gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                        >  2:3350326/1‑36 (MQ=255)
         tAAACAGAAATAAAGAGCTGACAGCACTATGTTCGa                       <  2:413743/36‑1 (MQ=255)
          aaaCAGAAATAAAGAGCTGACAGCACTATGTTCGaa                      <  2:152870/36‑1 (MQ=255)
          aaaCAGAAATAAAGAGCTGACAGCACTATGTTCGaa                      <  1:2075911/36‑1 (MQ=255)
          aaaCAGAAATAAAGAGCTGACAGCACTATGTTCGaa                      <  2:3246164/36‑1 (MQ=255)
          aaaCAGAAATAAAGAACTCACAGCACTATGTTCGaa                      <  1:466994/36‑1 (MQ=21)
           aaCAGAAATAAAGAGCTGACAGCACTATGTTCGaac                     >  1:1701343/1‑36 (MQ=255)
           aaCAGAAATAAAGAGCTGACAGCACTATGTTCGaac                     >  1:1178615/1‑36 (MQ=255)
            aCAGAAATAAAGAGCTGACAGCACTATGTTCGaaca                    >  1:3297216/1‑36 (MQ=255)
            aCAGAAATAAAGAGCTGACAGCACTATGTTCGaaca                    >  1:3114524/1‑36 (MQ=255)
             cAGAAATAAAGAGCTGACAGCACTATGTTCGaacaa                   <  1:1111307/36‑1 (MQ=255)
             cAGAAATAAAAAACTGACAGCACTATGTTCGaacaa                   <  1:551724/36‑1 (MQ=25)
              aGAAATAAAGAGCTGACAGCACTATGTTCGaacaac                  >  1:194947/1‑36 (MQ=255)
               gAAATAAAGAGCTGACAGCACTATGTTCGAACAAcg                 >  2:3556245/1‑36 (MQ=255)
                 aaTAAAGAGCTGACAGCACTATGTTCGAACAAcgcg               >  2:1073468/1‑36 (MQ=255)
                  aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt              <  1:1699773/36‑1 (MQ=255)
                  aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt              <  1:256298/36‑1 (MQ=255)
                   tAAAGAGCTGACAGCACTATGTTCGAACAACGCGTa             >  2:2666567/1‑36 (MQ=255)
                    aaaGTGCTGACAGCACTATGTTCGAACAACGCGTaa            <  1:713281/36‑1 (MQ=38)
                    aaaGGGCTGACAGCACTATGTTCGAACAACGCGTaa            <  1:517003/36‑1 (MQ=38)
                     aaaagCTGACAGCACTATGTTCGAACAACGCGTaaa           <  1:3322286/33‑1 (MQ=38)
                     aaGAGCTGACAGCACTATGTTCGAACAACGCGTaaa           >  2:2521065/1‑36 (MQ=255)
                       gagCTGACAGCACTATGTTCGAACAACGCGTTAAtt         >  1:2812441/1‑36 (MQ=37)
                       gagCTGACAGCACTATGTTCGAACAACGCGTAAAtt         >  2:166040/1‑36 (MQ=255)
                        agCTGACAGCACTATGTTCGAACAACGCGTCAATTc        >  2:2616639/1‑36 (MQ=38)
                        agCTGACAGCACTATGTTCGAACAACGCGTAAATTc        >  1:800911/1‑36 (MQ=255)
                          cTGACAGCACTATGTTCGTACAACGCGTTTATTCTg      >  2:2789739/1‑36 (MQ=19)
                           tGACAGCACTATGTTCGAACAACGCGTAAATTCTGa     <  1:2883435/36‑1 (MQ=255)
                            gACAGCACTATGTTCGAACAACGCGTAAATTCTGAc    >  2:1557893/1‑36 (MQ=255)
                            gACAGCACTATGTTCGAACAACGCGTAAATTCTGAc    <  2:2875182/36‑1 (MQ=255)
                             aCAGCACTATGTTCGAACAACGCGTAAATTCTGatg   >  1:1688229/1‑34 (MQ=38)
                              cAGCACTATGTTCGAACAACGCGTAAATTCTGACGt  >  1:1905883/1‑36 (MQ=255)
                              cAGCACTATGTTCGAACAACGCGTAAATTCTGACGt  >  1:618605/1‑36 (MQ=255)
                                 |                                
AATTTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTCTGACGT  >  REL606/3339125‑3339190

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG