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.T96.1% 810.0 / 20.6 206intergenic (‑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 . (5/3);  new base T (87/111);  total (92/114)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 241 total aligned reads displayed.

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875550‑3875695
                                                                              |                                                                        
aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                           >  2:5997712/1‑76 (MQ=255)
  aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAACAACTACA‑‑tt                                                                       >  1:333295/1‑76 (MQ=255)
   gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttt                                                                      <  2:2684438/76‑1 (MQ=255)
     attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttt                                                                    <  2:643936/76‑1 (MQ=255)
     attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttt                                                                    >  1:1613613/1‑76 (MQ=255)
     attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttt                                                                    >  2:4039480/1‑76 (MQ=255)
     attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttt                                                                    <  1:7651549/76‑1 (MQ=255)
      ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑tttttt                                                                   >  2:3068392/1‑76 (MQ=255)
      ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑tttttt                                                                   >  2:7277105/1‑76 (MQ=255)
       tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttttt                                                                 >  2:5279702/1‑76 (MQ=255)
       tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttttt                                                                 >  2:2784149/1‑76 (MQ=255)
        attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑tttttttt                                                                >  1:4364087/1‑75 (MQ=255)
        attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑tttttttt                                                                >  2:55306/1‑75 (MQ=255)
        attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTc                                                                >  1:3755283/1‑76 (MQ=255)
         ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttttttc                                                               >  1:1215400/1‑74 (MQ=255)
          tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttttttca                                                              <  1:6039096/76‑4 (MQ=255)
             aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                            <  1:7903603/76‑1 (MQ=255)
             aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                            <  1:5892759/76‑1 (MQ=255)
               aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                          <  1:6646758/76‑1 (MQ=255)
                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                        >  1:4653034/1‑76 (MQ=255)
                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                        >  1:4200464/1‑76 (MQ=255)
                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                        <  1:929785/76‑1 (MQ=255)
                  cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                                       >  1:5363851/1‑76 (MQ=255)
                   gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                      >  1:5606021/1‑76 (MQ=255)
                   gcgTCCTGACTGTAATGTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                      <  2:7233356/76‑1 (MQ=255)
                    cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                     <  1:3922540/76‑1 (MQ=255)
                    cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                     >  1:965846/1‑76 (MQ=255)
                    cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                     >  2:8033776/1‑76 (MQ=255)
                     gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    <  1:5442357/76‑1 (MQ=255)
                     gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    <  1:1277769/76‑1 (MQ=255)
                     gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    <  1:1787925/76‑1 (MQ=255)
                       ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                  >  2:3771382/1‑76 (MQ=255)
                       ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                  <  1:5508984/76‑1 (MQ=255)
                         tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                                <  1:182518/76‑1 (MQ=255)
                           aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                              <  2:2433443/76‑1 (MQ=255)
                           aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                              <  2:2250681/76‑1 (MQ=255)
                           aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                              >  2:4150736/1‑76 (MQ=255)
                           aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                              >  2:4348381/1‑76 (MQ=255)
                            cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                             <  1:7021632/76‑1 (MQ=255)
                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            <  1:4236200/76‑1 (MQ=255)
                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            <  2:7058940/76‑1 (MQ=255)
                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            >  2:4717353/1‑76 (MQ=255)
                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            >  1:7351982/1‑76 (MQ=255)
                              gTAATTTAAAACCGAATTTAATCCGCTTCACCTTGAAAAAACTCCA‑‑TTTTTTTCATCGTGAAATGGACAGAGGATa                                          <  2:5923453/76‑1 (MQ=255)
                               tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                          >  1:4163648/1‑76 (MQ=255)
                               gcaTTTGAGGACAAATTTTAAACGGATCCGCCTGAAAAAAATTACTT‑‑TTTTTTTCATCGTGAAATGGACAGAGGAt                                           <  1:3001778/74‑1 (MQ=255)
                                aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         <  2:492702/76‑1 (MQ=255)
                                aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         <  2:6419036/76‑1 (MQ=255)
                                 aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                        >  1:2533312/1‑76 (MQ=255)
                                  tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       <  2:2533312/76‑1 (MQ=255)
                                  tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       >  1:2249276/1‑76 (MQ=255)
                                   ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      <  2:3526060/76‑1 (MQ=255)
                                    tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                     <  2:6095140/76‑1 (MQ=255)
                                    tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                     <  2:5066721/76‑1 (MQ=255)
                                       gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  1:7240191/76‑1 (MQ=255)
                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  2:6814400/76‑1 (MQ=255)
                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  1:281114/76‑1 (MQ=255)
                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  2:2347419/76‑1 (MQ=255)
                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  1:4583537/76‑1 (MQ=255)
                                         aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAGTTGTTCaa                                <  2:3304973/76‑1 (MQ=255)
                                         aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                                <  2:3168044/76‑1 (MQ=255)
                                          cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                               <  2:1168012/76‑1 (MQ=255)
                                           gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              <  1:5816061/76‑1 (MQ=255)
                                           gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              <  1:5085933/76‑1 (MQ=255)
                                            aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                             <  1:102281/76‑1 (MQ=255)
                                             aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            >  2:51057/1‑76 (MQ=255)
                                               tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           >  2:3998014/1‑76 (MQ=255)
                                               tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           >  2:185404/1‑76 (MQ=255)
                                                ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                          <  1:1562914/76‑1 (MQ=255)
                                                ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                          >  1:6986435/1‑76 (MQ=255)
                                                 tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                         >  2:2635432/1‑76 (MQ=255)
                                                 tAAACCGCTTAACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                         >  2:2968150/1‑76 (MQ=255)
                                                  aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        >  1:1497956/1‑76 (MQ=255)
                                                   aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                       <  1:7591117/76‑1 (MQ=255)
                                                   aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                       <  2:6755819/76‑1 (MQ=255)
                                                   aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                       <  1:1935471/76‑1 (MQ=255)
                                                   aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                       <  1:4837544/76‑1 (MQ=255)
                                                   aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                       <  1:3429137/76‑1 (MQ=255)
                                                    aCCTCTTCACCTTGAAAAGACTACATTTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                       <  2:3351414/76‑1 (MQ=255)
                                                    aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                      >  1:988213/1‑76 (MQ=255)
                                                    aCCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTCTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                      <  1:2561914/76‑1 (MQ=255)
                                                     ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                     <  1:413780/76‑1 (MQ=255)
                                                      cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                    >  2:7758380/1‑76 (MQ=255)
                                                        cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                  >  1:7175868/1‑76 (MQ=255)
                                                        cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                  >  1:3586027/1‑76 (MQ=255)
                                                         ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc                 >  2:452470/1‑76 (MQ=255)
                                                            aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt              >  1:4305080/1‑76 (MQ=255)
                                                            aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt              >  2:5468432/1‑76 (MQ=255)
                                                             ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg             >  1:1192589/1‑76 (MQ=255)
                                                              cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa            >  1:8050359/1‑76 (MQ=255)
                                                               ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           >  2:1795207/1‑76 (MQ=255)
                                                                tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg          >  1:2031806/1‑76 (MQ=255)
                                                                tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg          >  1:1392122/1‑76 (MQ=255)
                                                                 gAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa         >  1:2376079/1‑76 (MQ=255)
                                                                 gAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa         >  1:1918756/1‑76 (MQ=255)
                                                                  aacaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa        <  2:5339539/73‑1 (MQ=255)
                                                                   aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa       <  1:3681067/76‑1 (MQ=255)
                                                                   aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa       >  2:1026098/1‑76 (MQ=255)
                                                                     aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc     <  1:292212/76‑1 (MQ=255)
                                                                        cTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  1:6404505/1‑76 (MQ=255)
                                                                              |                                                                        
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875550‑3875695

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: