Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,875,632 (T)7→8 intergenic (‑66/+287) glmU ← / ← atpC bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,875,6251.T97.7% 299.9 / 1.1 88intergenic (‑59/+294)glmU/atpCbifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (50/36);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 122 total aligned reads displayed.

GAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAA  >  REL606/3875589‑3875656
                                      |                               
gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                  >  2:10011937/1‑37 (MQ=255)
gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                  <  1:6621262/37‑1 (MQ=255)
 aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACAt                                >  1:11185962/1‑37 (MQ=255)
 aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACAt                                >  2:4667465/1‑37 (MQ=255)
 aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACAt                                <  2:1160946/37‑1 (MQ=255)
 aaCAAATTTAAGCCGCTTCACCTTGAAAAAACTACAt                                <  1:8363697/37‑1 (MQ=37)
  aCGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                               <  2:3135307/37‑1 (MQ=255)
   cGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                              <  2:2734703/37‑1 (MQ=255)
   cGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                              <  1:10724310/37‑1 (MQ=255)
   cGAATTAAAACCGCTTCACCTTGAAAAAACTACAttt                              <  1:4837531/37‑1 (MQ=255)
    gAATTTAAACCGCTTCACCTTGAAAAAACTACAtttt                             >  2:4021883/1‑37 (MQ=255)
    gAATTTAAACCGCTTCACCTTGAAAAAACTACAtttt                             <  1:3129304/37‑1 (MQ=255)
    gAATTTAAACCGCTTCACCTTGAAAAAACTACAtttt                             >  2:5418152/1‑37 (MQ=255)
    gAATTTAAACCGCTTCACCTTGAAAAAACTACAtttt                             <  1:7991594/37‑1 (MQ=255)
    gAATTTAAACCGCTTCACCTTGAAAAAACTACAtttt                             <  1:1626909/37‑1 (MQ=255)
    gAATTTAAACCGCTGCACCTTGAAAAAACTACAtttt                             >  1:2152732/1‑37 (MQ=255)
    gAATTAAAACCGCTTCACCTTGAAAAAACTACAtttt                             <  1:4694439/37‑1 (MQ=255)
     aaTTTAAACCGCTTCACCTTGAAAAAACTACAttttt                            >  2:1228027/1‑37 (MQ=255)
     aaTTTAAACCGCTTCACCTTGAAAAAACTACAttttt                            >  1:8225537/1‑37 (MQ=255)
     aaTTTAAACCGCTTCACCTTGAAAAAACTACAttttt                            <  1:9386725/37‑1 (MQ=255)
     aaTTTAAACCGCTTCACCTTGAAAAAACTACAttttt                            <  2:9379304/37‑1 (MQ=255)
      aTTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                           <  2:285990/37‑1 (MQ=255)
      aTTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                           <  2:8256814/37‑1 (MQ=255)
      aTTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                           >  2:948232/1‑37 (MQ=255)
      aTTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                           >  2:2991255/1‑37 (MQ=255)
      aTTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                           >  2:7608443/1‑37 (MQ=255)
      aTTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                           <  1:3999500/37‑1 (MQ=255)
       tttAAACCGCTTCACCTTGAAAAAACTACAttttttt                          >  1:9504876/1‑37 (MQ=255)
        ttAAGCCGCTTCACCTTGAAAAAACTACAtttttttt                         <  2:3146877/37‑2 (MQ=38)
        ttAAACCGCTTCACCTTGAAAAAACTACAtttttttt                         >  2:9163945/1‑36 (MQ=255)
              cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                    <  2:4986854/37‑1 (MQ=37)
              cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                    >  1:9724783/1‑37 (MQ=37)
              cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                    >  1:2384186/1‑37 (MQ=37)
              cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                    >  1:2413686/1‑37 (MQ=37)
              cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                    >  1:5273778/1‑37 (MQ=37)
              cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                    >  1:4376679/1‑37 (MQ=37)
               gCTTCACCTTGAAAAAACTACATTTTTTTTCNTCGTg                   <  2:10613718/37‑1 (MQ=37)
               gCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                   <  2:8764796/37‑1 (MQ=37)
               gCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                   >  1:2913823/1‑37 (MQ=37)
                cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                  >  1:7792418/1‑37 (MQ=37)
                cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                  >  2:2335729/1‑37 (MQ=37)
                 ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                 <  2:2915962/37‑1 (MQ=37)
                 ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                 <  2:4703711/37‑1 (MQ=37)
                 ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                 <  1:9504504/37‑1 (MQ=37)
                 ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                 <  1:7024181/37‑1 (MQ=37)
                 ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                 <  1:11140892/37‑1 (MQ=37)
                 ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                 >  2:7312756/1‑37 (MQ=37)
                  tCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                >  2:660092/1‑37 (MQ=37)
                  ccACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                >  1:5158732/2‑37 (MQ=37)
                   cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt               <  1:3663191/37‑1 (MQ=37)
                   cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt               >  2:4029145/1‑37 (MQ=37)
                   cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt               <  1:7538234/37‑1 (MQ=37)
                   cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt               >  2:7705781/1‑37 (MQ=37)
                   cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt               <  1:6616661/37‑1 (MQ=37)
                    aCCTTGAAAAAACTACGTTTTTTTTCATCGTGAAATg              <  1:534331/37‑1 (MQ=25)
                    aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg              >  1:2514226/1‑37 (MQ=37)
                    aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg              >  1:1811166/1‑37 (MQ=37)
                    aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg              >  2:10791677/1‑37 (MQ=37)
                     ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg             <  1:5438818/37‑1 (MQ=37)
                     ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg             >  2:1196125/1‑37 (MQ=37)
                     ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg             <  2:10376993/37‑1 (MQ=37)
                     ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg             >  1:4698800/1‑37 (MQ=37)
                     ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg             >  2:6061674/1‑37 (MQ=37)
                     ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg             <  1:1170283/37‑1 (MQ=37)
                     ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg             <  2:4007016/37‑1 (MQ=37)
                      cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGa            >  1:10404480/1‑37 (MQ=37)
                      cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGa            >  2:8222079/1‑37 (MQ=37)
                      cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGa            >  1:7572830/1‑37 (MQ=37)
                       ttGAAAAAACTACATTTTTTTCATCGTGAAATGGACa          >  1:8782687/1‑37 (MQ=255)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          <  2:9721241/37‑1 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          >  2:620158/1‑37 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          <  2:6543760/37‑1 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          >  2:2452090/1‑37 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          >  2:10446051/1‑37 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          >  1:4608416/1‑37 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          <  1:9681009/37‑1 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          <  1:5681578/37‑1 (MQ=37)
                        tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa          >  1:6553635/1‑37 (MQ=37)
                         gAAAAAACTACATTTTTTTTCATCGTGAAATGGACag         >  1:2243943/1‑37 (MQ=37)
                         gAAAAAACTACATTTTTTTTCATCGTGAAATGGACag         >  1:2843210/1‑37 (MQ=37)
                         gAAAAAACTACATTTTTTTTCATCGTGAAATGGACag         >  2:10763302/1‑37 (MQ=37)
                         gAAAAAACTACATTTTTTTTCATCGTGAAATGGACag         >  2:5040815/1‑37 (MQ=37)
                         gAAAAAACTACATTTTTTTTCATCGTGAAATGGACag         <  1:3002303/37‑1 (MQ=37)
                         gAAAAAACTACATTTTTTTTCATCGTGAAATGGACag         <  2:10085916/37‑1 (MQ=37)
                          aaaaaaCTACATTTTTTTTCATCGTGAAATGGACaga        <  2:6904331/37‑1 (MQ=37)
                           aaaaaCTACATTTTTTTTTATCGTGAAATGGACagag       <  1:10502643/37‑1 (MQ=25)
                           aaaaaCTACATTTTTTTTCATCGTGAAATGGACagag       <  2:6095598/37‑1 (MQ=37)
                            aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAgg      <  2:5740691/37‑1 (MQ=37)
                             aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGa     >  2:621003/1‑37 (MQ=37)
                             aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGa     <  1:422542/37‑1 (MQ=37)
                             aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGa     <  2:452024/37‑1 (MQ=37)
                             aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGa     >  2:3599848/1‑37 (MQ=37)
                              aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGAt    >  1:6674708/1‑37 (MQ=37)
                              aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGAt    >  2:8086825/1‑37 (MQ=37)
                              aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGAt    >  2:1843255/1‑37 (MQ=37)
                              aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGAt    >  1:1819449/1‑37 (MQ=37)
                              aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGAt    >  1:1557840/1‑37 (MQ=37)
                              aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGAt    >  2:9502474/1‑37 (MQ=37)
                                cTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa  <  2:6969091/37‑1 (MQ=37)
                                cTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa  <  1:8236712/37‑1 (MQ=37)
                                      |                               
GAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAA  >  REL606/3875589‑3875656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: