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.T82.0% 116.1 / 21.1 39intergenic (‑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 . (6/1);  new base T (18/14);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC  >  REL606/3875551‑3875694
                                                                           |                                                                          
tAGCATTATTCAACNTACGCGTCCTGACTGTAATTTGAGAACTAATTTAAACCGCTTCACCTTAAAAAAACTACA‑‑t                                                                          <  1:1572590/76‑1 (MQ=255)
      tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑ttttttt                                                                 <  2:1127649/76‑1 (MQ=255)
             cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAGCCGCTTCACCTTGAAAAAACTACA‑‑TTCTTTTTCATCGt                                                          <  1:96203/76‑1 (MQ=255)
               tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                        >  2:1362227/1‑76 (MQ=255)
                 cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                                      <  1:1209082/76‑1 (MQ=255)
                  gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                     >  1:340269/1‑76 (MQ=255)
                         gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                              <  1:1628711/76‑1 (MQ=255)
                         gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTTTATTTCATCGTGAAATGGACaga                                              >  2:738572/1‑76 (MQ=255)
                          aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                             <  1:1386655/76‑1 (MQ=255)
                            tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                           >  1:533972/1‑76 (MQ=255)
                            tGTAATTTGAGAACGAATCTANACNGCTTCACGTTGAAAAAACTACN‑‑TTTTTTTNCNTCGTGAAATGGACAGAGGa                                           >  2:266340/1‑76 (MQ=255)
                             gTAATTTGAGAACGAATTTAGACCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTGCATCGTGAAATGGACAGAGGAt                                          >  2:153753/1‑76 (MQ=255)
                             gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                          >  2:1391379/1‑76 (MQ=255)
                             gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTCTTTTTCATCGTGAAATGGACAGAGGAt                                          >  1:176745/1‑76 (MQ=255)
                              tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                         <  1:637491/76‑1 (MQ=255)
                              tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                         >  1:758192/1‑76 (MQ=255)
                              tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                         >  2:1068838/1‑76 (MQ=255)
                              tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                         >  1:1174348/1‑76 (MQ=255)
                              tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                         <  2:533971/76‑1 (MQ=255)
                              tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                         >  1:1372706/1‑76 (MQ=255)
                               aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                        >  2:1736384/1‑76 (MQ=255)
                               aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                        <  1:956056/76‑1 (MQ=255)
                                aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                       <  2:871982/76‑1 (MQ=255)
                                  ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                     >  1:904787/1‑76 (MQ=255)
                                   tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAGCTACA‑‑TNTTTTTTCATCGTGAAATGNACNGAGGATAAATTg                                    >  1:1492033/1‑76 (MQ=255)
                                   tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                    >  2:1324975/1‑76 (MQ=255)
                                   tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                    >  1:836671/1‑76 (MQ=255)
                                     agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                   >  2:1355861/1‑76 (MQ=255)
                                         cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATTAATTGTTCAAt                              <  1:100458/76‑1 (MQ=255)
                                          gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                             >  2:877348/1‑76 (MQ=255)
                                              ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                         <  1:529148/76‑1 (MQ=255)
                                              ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                         <  2:806807/76‑1 (MQ=255)
                                                 aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                      <  2:1772050/76‑1 (MQ=255)
                                                    cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                   >  1:1746487/1‑76 (MQ=255)
                                                    cGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                  >  1:1737959/1‑76 (MQ=255)
                                                        tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg               <  2:103165/76‑1 (MQ=255)
                                                        tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg               >  1:419638/1‑76 (MQ=255)
                                                        tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg               >  1:1706307/1‑76 (MQ=255)
                                                          aCCGTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt             >  1:1691819/1‑76 (MQ=255)
                                                                aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa       <  1:205690/76‑1 (MQ=255)
                                                                aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa       <  2:753355/76‑1 (MQ=255)
                                                                     aCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc  <  1:632749/76‑1 (MQ=255)
                                                                           |                                                                          
TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC  >  REL606/3875551‑3875694

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