Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,875,625 +T 100% intergenic (‑59/+294) glmU ← / ← atpC bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,875,6251.→T100.0% 67.4 18intergenic (‑59/+294)glmU/atpCbifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  new base (12/6):  ref base (0/0):  total (12/6)

AACGAATTTAAACCGCTTCACCTTGAAAAAACTACA.TTTTTTTCATCGTGAAATGGACAGAGGATA  >  REL606/3875590‑3875655
                                    |                              
aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                 <  1:3511380/36‑1 (MQ=255)
 aCGAATTTAAACCGCTTCACCTTGAAAAAACTACA.t                               >  1:5948736/1‑36 (MQ=255)
 aCGAATTTAAACCGCTTCACCTTGAAAAAACTACA.t                               <  1:591143/36‑1 (MQ=255)
   gAATTTAAACCGCTTCACCTTGAAAAAACTACA.ttt                             <  1:1899107/36‑1 (MQ=255)
   gAATTTAAACCGCTTCACCTTGAAAAAACTACA.ttt                             <  1:2618489/36‑1 (MQ=255)
    aaTTTAAACCGCTTCACCTTGAAAAAACTACA.tttt                            <  1:1770775/36‑1 (MQ=255)
    aaTTTAAACCGCTTCACCTTGAAAAAACTACA.tttt                            <  1:6115647/36‑1 (MQ=255)
     aTTTAAACCGCTTCACCTTGAAAAAACTACA.ttttt                           <  1:217580/36‑1 (MQ=255)
     aTTTAAACCGCTTCACCTTGAAAAAACTACA.ttttt                           >  1:506671/1‑36 (MQ=255)
     aTTTAAACCGCTTCACCTTGAAAAAACTACA.ttttt                           <  1:5050001/36‑1 (MQ=255)
      tttAAACCGCTTCACCTTGAAAAAACTACA.tttttt                          <  1:5708515/36‑1 (MQ=255)
      tttAAACCGCTTCACCTTGAAAAAACTACA.tttttt                          >  1:2647726/1‑36 (MQ=255)
       ttAAACCGCTTCACCTTGAAAAAACTACA.ttttttt                         >  1:3124000/1‑36 (MQ=255)
        tAAACCGCTTCACCTTGAAAAAACTACA.tttttttt                        >  1:5322343/1‑35 (MQ=255)
        tAAACCGCTTCACCTTGAAAAAACTACA.tttttttt                        >  1:401764/1‑35 (MQ=255)
        tAAACCGCTTCACCTTGAAAAAACTACA.tttttttt                        >  1:3753845/1‑35 (MQ=255)
        tAAACCGCTTCACCTTGAAAAAACTACA.tttttttt                        >  1:7604827/1‑35 (MQ=255)
         aaaCCGCTTCACCTTGAAAAAACTACA.ttttttttc                       >  1:4458856/1‑34 (MQ=255)
         aaaCCGCTTCACCTTGAAAAAACTACA.ttttttttc                       <  1:4690319/36‑3 (MQ=255)
         aaaCCGCTTCACCTTGAAAAAACTACA.ttttttttc                       <  1:2553574/36‑3 (MQ=255)
             cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                    <  1:7857330/36‑1 (MQ=37)
              gCTTAACCTTGAAAAAACTAAATTTTTTTTCATCGt                   >  1:6977653/1‑36 (MQ=19)
               cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                  >  1:4758057/1‑36 (MQ=37)
               cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                  >  1:7853959/1‑36 (MQ=37)
                 tCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                >  1:477533/1‑36 (MQ=37)
                  cACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa               <  1:6566118/36‑1 (MQ=37)
                  cACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa               >  1:6088581/1‑36 (MQ=37)
                  cACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa               <  1:1320005/36‑1 (MQ=37)
                   aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt              <  1:5450222/36‑1 (MQ=37)
                    ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATg             >  1:5690119/1‑36 (MQ=37)
                    ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATg             >  1:7686628/1‑36 (MQ=37)
                     cTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg            >  1:706942/1‑36 (MQ=37)
                        gAAAAAACTACATTTTTTTTCATCGTGAAATGGACa         >  1:1605474/1‑36 (MQ=37)
                        gAAAAAACTACATTTTTTTTCATCGTGAAATGGACa         <  1:926401/36‑1 (MQ=37)
                          aaaaaCTACATTTTTTTTCATCGTGAAATGGACaga       <  1:784912/36‑1 (MQ=37)
                           aaaaCTACATTTTTTTTCATCGTGAAATGGACagag      >  1:3653582/1‑36 (MQ=37)
                               cTACATTTTTTTTCATCGTGAAATGGACAGAGGATa  >  1:4846793/1‑36 (MQ=37)
                               cTACATTTTTTTTCATCGTGAAATGGACAGAGGATa  >  1:3844561/1‑36 (MQ=37)
                                    |                              
AACGAATTTAAACCGCTTCACCTTGAAAAAACTACA.TTTTTTTCATCGTGAAATGGACAGAGGATA  >  REL606/3875590‑3875655

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG