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

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

AAGGGGACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGTCGCTGCCAGGCTT  >  REL606/3490079‑3490144
                                 |                                
aaGGGGACAGTATAAAGCGTTACGCGCCGTACGTCa                                >  1:1650378/1‑36 (MQ=255)
aaGCGGACAGTATAAAGCGTTACGCGCCGTACGTCa                                >  2:488627/1‑36 (MQ=38)
 aGGGGACAGTATAAAGCGTTACGCGCCGTACGTCAc                               >  1:1867902/1‑36 (MQ=255)
   gggACAGTATAAAGCGTTACGCGCCGTACGTCACct                             <  1:745871/36‑1 (MQ=255)
    ggACAGTATAAAGCGTTACGCGCCGTACGTCACctc                            >  2:2978429/1‑36 (MQ=255)
      aCAGTATAAAGCGTTACGCGCCGTACGTCACCTCTg                          <  1:3495867/36‑1 (MQ=255)
      aCAGTATAAAGCGTTACGCGCCGTACGTCACCTCTg                          <  2:1478613/36‑1 (MQ=255)
       cAGTATAAAGCGTTACGCGCCGTACGTCACCTCTGc                         >  1:1326779/1‑36 (MQ=255)
        aGTATAGAACGTTACGCGGCGTACGTCACCTCTGCg                        <  1:877007/36‑1 (MQ=19)
        aGTATAAAGCGTTACGCGCCGTACGTCACCTCTGCg                        >  1:1969532/1‑36 (MQ=255)
        aGTATAAAGCGTTACGCGCCGTACGTCACCTCTGCg                        <  2:2675170/36‑1 (MQ=255)
          tataAAGCGTTACGCGCCGTACGTCACCTCTGCggg                      <  1:830706/36‑1 (MQ=255)
            taAAGCGTTACGCGCCGTACGTCACCTCTGCGGGaa                    <  1:2057461/36‑1 (MQ=255)
            taAAGCGTTACGCGCCGTACGTCACCTCTGCGGGaa                    <  2:1602401/36‑1 (MQ=255)
            taAAGCGTTACGCGCCGTACGTCACCTCTGCGGGaa                    <  2:1333177/36‑1 (MQ=255)
             aaaGCGTTATGCGCCGTACGTCACCTCTGCGGGaaa                   >  1:3040800/1‑36 (MQ=37)
             aaaGCGTTACGCGCCGTACGTCACCTCTGCGGGaaa                   >  2:1858085/1‑36 (MQ=255)
             aaaCCGCTACGCGCCGTACATCACCTCTGCGGGaca                   <  2:1658337/36‑3 (MQ=16)
              aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  >  2:1809723/1‑36 (MQ=255)
              aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  <  2:2859966/36‑1 (MQ=255)
              aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  >  1:1143469/1‑36 (MQ=255)
              aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  >  1:1473750/1‑36 (MQ=255)
               aGCGTGACGCGCCGTACGTCACCTCTGCGGGAAACt                 <  1:3315203/36‑1 (MQ=38)
                gCGTTACGCGCCGTACGTCACCTCTGCGGGGAACTg                >  2:2369164/1‑36 (MQ=38)
                gCGTTACGCGCCGTACGTCACCTCTGCGGGAAACTg                <  2:3726749/36‑1 (MQ=255)
                gCGTTACGCGCCGTACGTCACCTCTGCGGGAAACTg                >  1:281499/1‑36 (MQ=255)
                 cGTTTCGCGCCGTACGTCAGCTCTGCGGGAAACTGt               <  1:2678739/36‑1 (MQ=25)
                  gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              <  2:680625/36‑1 (MQ=255)
                  gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  2:1753799/1‑36 (MQ=255)
                  gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  1:2943615/1‑36 (MQ=255)
                  gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  1:3663757/1‑36 (MQ=255)
                  gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  1:418318/1‑36 (MQ=255)
                   ttACGCGCCGTACGTCTCCTCTGCGGGAAACTGTCg             <  2:2972056/36‑1 (MQ=38)
                   ttACGCGCCGTACGTCACCTCTGCGGGAAACTGTCg             >  2:1445797/1‑36 (MQ=255)
                   ttACGCGCCGTACGTCACCTCTGCGGGAAACTGTCg             >  1:673623/1‑36 (MQ=255)
                    tACGCGCCGTACGTCACCTCTGCGGGAAACTGTCGc            >  2:650616/1‑36 (MQ=255)
                     aCGCGCCGTACGTCACCTCTGCGGGAAACTGTCGCt           >  2:9526/1‑36 (MQ=255)
                     aCGCGCCGTACGTCACCTCTGCGGGAAACTGTCGCt           <  1:1881160/36‑1 (MQ=255)
                     aCGCGCCGTACGTCACCTCTGCGCCAAACTGTCGCt           >  1:2879757/1‑36 (MQ=25)
                     aCGCGCCGTACGGCACCTCTGCGGGAAACTGTCGCt           >  2:3584773/1‑36 (MQ=255)
                     aCGCCCCGTACGTCACCTCTGCGGGAAACTGTCGCt           <  1:2439653/36‑1 (MQ=38)
                      cgcgCCGTACGTCACCTCTGCGGGAAACTGTCGCTg          <  1:2704761/36‑1 (MQ=255)
                      cgcgCCGTACGTCACCTCTCCGGGAAACTCTCGcac          >  1:977794/1‑34 (MQ=25)
                        cgcCGTACGTCACCTCTGCGGGAAACTGTCGCTGcc        <  1:1710121/36‑1 (MQ=255)
                         gcCGTACGTCGCTTCTGCGGGAAACTGTCGCTGCCa       <  1:529802/36‑1 (MQ=25)
                         gcCGTACGTCACCTCTGCGGGAAACTGTTGCTGCCa       >  2:2247685/1‑36 (MQ=37)
                         gcCGTACGTCACCTCTGCGGGAAACTGTCGCTGctt       >  1:660059/1‑34 (MQ=37)
                         gcCGTACGTCACCTCTGCGGGAAACTGTCGCTGCCa       >  2:2615451/1‑36 (MQ=255)
                          ccGTACGTCACCTCTGCGGGAAACTGTCGCTGCCAg      >  1:1360071/1‑36 (MQ=255)
                          ccGTACGTCACCTCTGCGGGAAACTGTCGCTACCAg      <  2:1082196/36‑1 (MQ=37)
                           cGTACGTCACCTCTGCGGGAAACTGTCGCTGCCAgg     <  2:3123424/36‑1 (MQ=255)
                            gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc    <  1:606757/36‑1 (MQ=255)
                            gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc    >  1:492415/1‑36 (MQ=255)
                            gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc    >  2:842592/1‑36 (MQ=255)
                             tACGTCACCTCTGCGGGAAACTGTCGCTGCCAAGCt   >  2:928302/1‑36 (MQ=37)
                             tACGCCACCCCTGCGGGAAACTGTCGCTGCCAGGCt   <  1:2247892/36‑1 (MQ=255)
                              aCGTCACCTCTGCGGGAAACTGTCGCTGCCAGGCtt  <  1:2285899/36‑1 (MQ=255)
                                 |                                
AAGGGGACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGTCGCTGCCAGGCTT  >  REL606/3490079‑3490144

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG