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.T95.0% 142.1 / 1.0 41intergenic (‑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 . (2/0);  new base T (16/22);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875552‑3875695
                                                                          |                                                                       
aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                      >  2:128346/1‑76 (MQ=255)
   attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                   >  1:795037/1‑76 (MQ=255)
     tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                                 >  1:710433/1‑76 (MQ=255)
      attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                                <  2:542858/76‑2 (MQ=255)
            cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                           >  1:80506/1‑76 (MQ=255)
              tACGCGTCCTGAATGTAATTTGAGAACGAATTNAAACCGCTNCACCTTGAAAAAACTACATNTTTTTTCATCGTGa                                                         >  2:962124/1‑76 (MQ=255)
               aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                        >  1:682587/1‑76 (MQ=255)
                   gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    >  2:140873/1‑76 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                  >  2:630376/1‑76 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                  <  1:585352/76‑1 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                  <  1:574023/76‑1 (MQ=255)
                      cTGACTGTAATTAGAGAACGAATTTAAACCGCNTCGCCTTGNAAAAACTACATTTTTTTTCNTCGTGAAATGGACa                                                 >  2:966348/1‑76 (MQ=255)
                       tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                                <  1:457254/76‑1 (MQ=255)
                       tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                                <  2:757687/76‑1 (MQ=255)
                        gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                               >  2:607255/1‑76 (MQ=255)
                         aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                              <  2:420303/76‑1 (MQ=255)
                          cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                             <  1:834719/76‑1 (MQ=255)
                          cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                             >  1:463317/1‑76 (MQ=255)
                           tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            <  2:709889/76‑1 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           >  1:194466/1‑76 (MQ=255)
                             tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGTAATGGACAGAGGATa                                          <  1:1042243/76‑1 (MQ=255)
                              aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         >  2:451680/1‑76 (MQ=255)
                                tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       <  1:630376/76‑1 (MQ=255)
                                 ctaagaACGAATTNAAACCGCTTCTCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      <  1:155364/73‑1 (MQ=255)
                                    agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                                   >  2:648760/1‑76 (MQ=255)
                                      aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  2:636112/76‑1 (MQ=255)
                                      aaCGAATTTAAACCGCTGCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGATAGAGGATAAATTGTTCa                                 <  1:128346/76‑1 (MQ=255)
                                         gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              <  2:272371/76‑1 (MQ=255)
                                         gAATTTAAACCGCNTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              <  2:1063399/76‑1 (MQ=255)
                                          aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                             <  2:279388/76‑1 (MQ=255)
                                             ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                          <  2:256087/76‑1 (MQ=255)
                                                  ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                     >  1:980368/1‑76 (MQ=255)
                                                     cTTCACCTNGAAAAAACTACATTTTTTTTCATCGTGAAAAGGACAGAGGATAAATTGTTCAATAACGGATTATccc                  <  1:1019212/76‑1 (MQ=255)
                                                      ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc                 <  1:228735/76‑1 (MQ=255)
                                                      ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc                 >  1:832528/1‑76 (MQ=255)
                                                      ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc                 <  1:537652/76‑1 (MQ=255)
                                                       tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  1:111114/1‑76 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAGATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa            <  1:194659/76‑1 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa            >  1:514657/1‑76 (MQ=255)
                                                            ttGAAAACACTACATTTTTTTTCATCGTGTAATGGACAGAGGATAAATTGTTTAATAATGGATTATCCCCGCTGAc           <  1:57540/76‑1 (MQ=255)
                                                            ttGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg          >  1:729957/1‑76 (MQ=255)
                                                               aaaaatCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa        >  2:76363/1‑76 (MQ=255)
                                                                 aaaaCTACACTTTTTTTCATCGTNAAATGNACAGAGGATAAATTNTTTAATCACGGATTATCCCCGCTGACGAAAg      >  2:145169/1‑76 (MQ=255)
                                                                  aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc     <  2:200319/76‑1 (MQ=255)
                                                                     cTACATTTTTTTTCATCGTGTAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  1:8587/1‑76 (MQ=255)
                                                                          |                                                                       
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875552‑3875695

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