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% 123.2 35intergenic (‑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 (21/13):  ref base (1/0):  total (22/13)

AACGAATTTAAACCGCTTCACCTTGAAAAAACTACA.TTTTTTTCATCGTGAAATGGACAGAGGATA  >  REL606/3875590‑3875655
                                    |                              
aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                 >  1:1041387/1‑36 (MQ=255)
aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                 >  1:1274321/1‑36 (MQ=255)
aaCGAATGTAAACCGCTTCACCTTGAAAAAAATACa                                 >  1:907302/1‑36 (MQ=37)
 aCGAATTTAAACCGCTTCACCTTGAAAAAACTACA.t                               >  1:2241192/1‑36 (MQ=255)
 aCGAATTTAAACCGCTTCACCTTGAAAAAACTACA.t                               >  1:2397956/1‑36 (MQ=255)
  cGAATTTAAACCGCTTCACCTTGAAAAAACTACA.tt                              >  1:1170350/1‑36 (MQ=255)
  cGAATTTAAACCGCTTCACCTTGAAAAAACTACA.tt                              >  1:7076816/1‑36 (MQ=255)
   gAATTTAAACCGCTTCACCTTGAAAAAACTACA.ttt                             <  1:74052/36‑1 (MQ=255)
     aTTTAAACCGCTTCACCTTGAAAAAACTACA.ttttt                           >  1:5562475/1‑36 (MQ=255)
     aTTTAAACCGCTTCACCTTGAAAAAACTACA.ttttt                           <  1:3425225/36‑1 (MQ=255)
      tttAAACCGCTTCACCTTGAAAAAACTACA.tttttt                          <  1:4404914/36‑1 (MQ=255)
      tttAAACCGCTTCACCTTGAAAAAACTACA.tttttt                          >  1:7602480/1‑36 (MQ=255)
      tttAAACCGCTTCACCTTGAAAAAACTACA.tttttt                          <  1:6056900/36‑1 (MQ=255)
         aaaCCGCTTCACCTTGAAAAAACTACA.ttttttttc                       >  1:5808169/1‑34 (MQ=255)
         aaaCCGCTTCACCTTGAAAAAACTACA.ttttttttc                       >  1:2533835/1‑34 (MQ=255)
         aaaCCGCTTCACCTTGAAAAAACTACA.ttttttttc                       >  1:7412230/1‑34 (MQ=255)
          aaCCGCTTCACCTTGAAAAAACTACA.ttttttttca                      >  1:2302902/1‑33 (MQ=255)
          aaCCGCTTCACCTTGAAAAAACTACA.ttttttttca                      >  1:292513/1‑33 (MQ=255)
          aaCCGCTTCACCTTGAAAAAACTACA.ttttttttca                      >  1:7249783/1‑33 (MQ=255)
          aaCCGCTTCACCTTGAAAAAACTACA.ttttttttca                      >  1:7641373/1‑33 (MQ=255)
          aaCCGCTTCACCTTGAAAAAACTACA.ttttttttca                      >  1:7782234/1‑33 (MQ=255)
             cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                    >  1:2843306/1‑36 (MQ=37)
             cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                    <  1:6906369/36‑1 (MQ=37)
             cGCCTCACCTTGAAAAAACTACATTTTTTTTCATCg                    <  1:709029/36‑1 (MQ=25)
               cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                  >  1:707258/1‑36 (MQ=37)
               cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                  >  1:647031/1‑36 (MQ=37)
                ttCACCTTGAAAAAACTACATTTTTTTTCATCGTgc                 <  1:4930217/36‑2 (MQ=37)
                ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                 >  1:3944741/1‑36 (MQ=37)
                 tCACCTTGAAAAAACTACATTTTTTTTCTTCGTGaa                >  1:6275395/1‑36 (MQ=25)
                 tCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                >  1:4615080/1‑36 (MQ=37)
                 tCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                >  1:3136141/1‑36 (MQ=37)
                  cACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa               >  1:946721/1‑36 (MQ=37)
                   aCCTTGAAAAAAATACATTTTTTTTCATCGTGAAAt              >  1:1593673/1‑36 (MQ=21)
                    ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATg             <  1:629074/36‑1 (MQ=37)
                    ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATg             >  1:1883164/1‑36 (MQ=37)
                     cTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg            >  1:3555208/1‑36 (MQ=37)
                     cTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg            <  1:6993565/36‑1 (MQ=37)
                     cTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg            >  1:1649272/1‑36 (MQ=37)
                      ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGa           >  1:1090933/1‑36 (MQ=37)
                       tGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc          >  1:4276323/1‑36 (MQ=37)
                       tGAAAAAACTACA.TTTTTTTCATCGTGAAATGGACa         >  1:1237718/1‑36 (MQ=255)
                       tGAAAAAAATACATTTTTTTTCATCATGAAATGGAc          >  1:3300985/1‑36 (MQ=16)
                        gAAAAAACTACATTTTTTTTCATCGTGAAATGGACa         <  1:3542639/36‑1 (MQ=37)
                        gAAAAAACTACATTTTTTTTCATCGTGAAATGGACa         >  1:7419608/1‑36 (MQ=37)
                        gAAAAAACTACATTTTTTTTCATCGTGAAATGGACa         <  1:1785076/36‑1 (MQ=37)
                         aaaaaaCTACATTTTTTTTCATCGTGAAATGGACag        <  1:6666090/36‑1 (MQ=37)
                          aaaaaCTACATTTTTTTTCATCGTGAAATGGACaga       <  1:7824367/36‑1 (MQ=37)
                          aaaaaCTACATTTTTTTTCATCGTGAAATGGACaga       >  1:1991003/1‑36 (MQ=37)
                          aaaaaCTACATTTTTTTTCATCGTGAAATGGACaga       <  1:2979167/36‑1 (MQ=37)
                          aaaaaCTACATTTTTTTTCATCGTGAAATGGACaga       >  1:6779967/1‑36 (MQ=37)
                          aaaaaCTACATTTTTTTTCATCGTGAAATGGACaga       >  1:3985760/1‑36 (MQ=37)
                           aaaaCTACATTTTTTTTCATCGTGAAATGGACagag      <  1:7546470/36‑1 (MQ=37)
                           aaaaCTACATTTTTTTTCATCGTGAAATGGACagag      >  1:1272736/1‑36 (MQ=37)
                            aaaCTACATTTTTTTTCATCGTGAAATGGACAGAgg     <  1:279763/36‑1 (MQ=37)
                              aCTACATTTTTTTTCATCGTGAAATGGACAGAGGAt   <  1:5478121/36‑1 (MQ=37)
                               cTACATTTTTTTTCATCGTGAAATGGACAGAGGATa  >  1:7298173/1‑36 (MQ=37)
                                    |                              
AACGAATTTAAACCGCTTCACCTTGAAAAAACTACA.TTTTTTTCATCGTGAAATGGACAGAGGATA  >  REL606/3875590‑3875655

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG