Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,490,112 C→T 100% A105T (GCG→ACG)  glpE ← thiosulfate sulfurtransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,490,1120C→T100.0% 176.5 49313glpEthiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  new base (24/25):  ref base (0/0):  total (24/25)

AAGGGGACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGTCGCTGCCAGGCT  >  REL606/3490079‑3490143
                                 |                               
aaGGGGACAGTATAAAGCGTTACGCGCCGTACGTCa                               >  1:3366597/1‑36 (MQ=255)
aaGGGGACAGTATAAAGCGTTACGCGCCGTACGTCa                               >  1:3445213/1‑36 (MQ=255)
aaGGGGACAGTATAAAGCGTTACGCGCCGTACGTCa                               >  1:2108955/1‑36 (MQ=255)
 aGGGGACAGTATAAAGCGTTACGCGCCGTACGTCAc                              >  1:1771310/1‑36 (MQ=255)
  ggggACATTATAAAGCGTTACGCGCCGTACGTCAcc                             <  1:4699605/36‑1 (MQ=37)
  ggggACAGTATAAAGCGTTACGCGCCGTACGTCAcc                             <  1:4881905/36‑1 (MQ=255)
   gggCCAGTCTAAAGCGTTACCCGCCTTACGTCACct                            <  1:5930424/36‑1 (MQ=16)
   gggACAGTATAAAGCGTTACGCGCCGTACGTCACct                            <  1:2449266/36‑1 (MQ=255)
   gggACAGTATAAAGCGTTACGCGCCGTACGTCACct                            <  1:3568066/36‑1 (MQ=255)
   gggACAGTATAAAGCGTTACGCGCCGGAAGTCACct                            >  1:3168853/1‑36 (MQ=25)
    ggACTGTATAAAGCGTTACGCGCCGTACGTCACctc                           <  1:4351624/36‑1 (MQ=38)
    ggACAGTATAAATCGTTGCGCGCCGTACGTCACctc                           <  1:2822021/36‑1 (MQ=37)
     gACAGTATAAAGCGTTACGCGCCGTACGTCACctct                          <  1:1606361/36‑1 (MQ=255)
     gACAGTATAAAGCGTTACGCGCCGTACGTCACctct                          <  1:4350451/36‑1 (MQ=255)
      aCAGTATAAAGCGTTACGCGCCGTACGTCACCTCTg                         <  1:4797194/36‑1 (MQ=255)
      aCAGTATAAAGCGTTACGCGCCGTACGTCACCTCTg                         >  1:5645627/1‑36 (MQ=255)
        aGTATAAAGCGTTACTCGCCGTACGTCACCTCTGCg                       <  1:932081/36‑1 (MQ=38)
        aGTATAAAGCGTTACGCGCCGTACGTCACCTCTGCg                       >  1:4428151/1‑36 (MQ=255)
           ataAAGCGTTACGCGCCGTACGTCACCTCTGCGGGa                    <  1:5298130/36‑1 (MQ=255)
            taAAGCGTTACGCGCCGTACGTCACCTCTGCGGGaa                   <  1:4917273/36‑1 (MQ=255)
             aaaGCGTTACGCGCCGTACGTCACCTCTGCGGGaaa                  <  1:5294150/36‑1 (MQ=255)
             aaaGCGTTACGCGCCGTACGTCACCTCTGCGGGaaa                  >  1:946/1‑36 (MQ=255)
              aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                 >  1:3109088/1‑36 (MQ=255)
              aaGCGTTACGCGCCGTACGTCACCTATGCGGGAAAc                 >  1:1242378/1‑36 (MQ=37)
              aaGCGTTACCCGCCGCACGTCACCTCGGCGGGGAAc                 >  1:5570017/1‑36 (MQ=16)
               aGCGTTACGCGCCGTACGTCCCCTCTGCGGGAAACt                <  1:86482/36‑1 (MQ=255)
               aGCGTTACGCGCCGTACGTCACCTCTGCGGGAAACt                <  1:3380451/36‑1 (MQ=255)
               aGCGTTACGCGCCGTACGTCACCTCTGCGGGAAACt                <  1:1457819/36‑1 (MQ=255)
                gCGTTACGCGCCGTACGTCCCCTCTGCGGGAAACTg               <  1:5670755/36‑1 (MQ=255)
                 cGTTACGCGCCGTACGTCACCTCTGCGGGAAACTGt              >  1:3008088/1‑36 (MQ=255)
                  gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc             >  1:3216799/1‑36 (MQ=255)
                  gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc             >  1:5156823/1‑36 (MQ=255)
                   ttACGCGCCGTACGTCACCTCTGCGGGAAACTGTcc            >  1:1777296/1‑35 (MQ=255)
                   ttACGCGCCGTACGTCACCTCTGCGGGAAACTGTCg            <  1:4318811/36‑1 (MQ=255)
                    tACGCGCCGTACGTCACCTCTGCGGGAAACTGTTGc           >  1:3013372/1‑36 (MQ=38)
                    tACGCGCCGTACGTCACCTCTGCGGGAAACTGTCGc           >  1:4086614/1‑36 (MQ=255)
                     aCGCGCCGTACGTCACCTCTGCGGGAAACTGTCGCt          <  1:1516313/36‑1 (MQ=255)
                      cgcgCCGTACGTCACCTCTGCGGGAAACTGTCGCTg         <  1:3902716/36‑1 (MQ=255)
                      cgcgCCGTACGTCACCTCTGCGGGAAACTGTCGCTg         <  1:1185711/36‑1 (MQ=255)
                       gcgcCGTACGTTACCGCTGCGGGAAACTGTCGCGGc        >  1:4104535/1‑36 (MQ=19)
                        cgcCGTACGTCACCTCTGCGGGAAACTGTCGCTGcc       <  1:1062356/36‑1 (MQ=255)
                         gcCGTACGTCACCTCTGCGGGAAACTGTCGCTGCCa      >  1:2198018/1‑36 (MQ=255)
                          ccGTACGTCACCTCTGCGGGAAACTGTCGCTGCCAg     >  1:4273674/1‑36 (MQ=255)
                          ccGTACGTCACCTCTGCGGGAAACTGTCGCTGCCAg     >  1:3894427/1‑36 (MQ=255)
                           cGTACGTCACCTCTGCGGGAAACTGCCGCTGCCAgg    >  1:2417365/1‑36 (MQ=37)
                            gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc   <  1:3943806/36‑1 (MQ=255)
                            gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc   <  1:1576740/36‑1 (MQ=255)
                            gTACGTCACCTCTGCGGGAAAATGTCGCTTCCAGGc   >  1:3654608/1‑36 (MQ=25)
                             tACGTCACCTCTCCGGGAAACTGTCGCTGCCAGGCt  >  1:3707331/1‑36 (MQ=37)
                                 |                               
AAGGGGACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGTCGCTGCCAGGCT  >  REL606/3490079‑3490143

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG