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.8% 186.2 / 0.5 48intergenic (‑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 . (0/2);  new base T (25/21);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG  >  REL606/3875551‑3875693
                                                                            |                                                                     
tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAt                                                                      <  1:912336/76‑1 (MQ=255)
 aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                     >  1:1527059/1‑76 (MQ=255)
     ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCANCTTGAAAAAACTACAtttttt                                                                >  2:301361/1‑76 (MQ=255)
     ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                >  1:1140858/1‑76 (MQ=255)
     ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                >  1:264020/1‑76 (MQ=255)
     ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                >  1:1683538/1‑76 (MQ=255)
     ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                <  2:1067751/76‑1 (MQ=255)
       attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                              >  2:1635461/1‑75 (MQ=255)
        ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                             <  1:216145/76‑3 (MQ=255)
            aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                          <  1:1501286/76‑1 (MQ=255)
             cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                         <  1:1256/76‑1 (MQ=255)
              aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                        <  1:136985/76‑1 (MQ=255)
              aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                        <  2:225955/76‑1 (MQ=255)
              aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                        <  1:1200584/76‑1 (MQ=255)
               tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                       >  2:668130/1‑76 (MQ=255)
                  gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                    <  2:831782/76‑1 (MQ=255)
                   cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                   >  2:181379/1‑76 (MQ=255)
                   cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                   >  1:499519/1‑76 (MQ=255)
                        tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                              <  1:396072/76‑1 (MQ=255)
                        tGACTGTAATTTGAGAACAAATTAAACCCGCTTCACCTTAAAAAAACTACTTTTTTTTTCATCGTGAAATGGACag                                              <  2:1029276/76‑1 (MQ=255)
                          aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                            >  1:1688325/1‑76 (MQ=255)
                          aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                            <  1:175268/76‑1 (MQ=255)
                             gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                         <  2:906453/76‑1 (MQ=255)
                              tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                        <  1:1316789/76‑1 (MQ=255)
                                aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                      >  1:1246924/1‑76 (MQ=255)
                                  ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                    >  1:398092/1‑76 (MQ=255)
                                    gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                  >  2:277553/1‑76 (MQ=255)
                                        aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                              >  2:745392/1‑76 (MQ=255)
                                         cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                             >  1:1224445/1‑76 (MQ=255)
                                          gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                            <  1:28429/76‑1 (MQ=255)
                                          gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                            >  1:1442145/1‑76 (MQ=255)
                                           aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                           <  1:671404/76‑1 (MQ=255)
                                              ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                        >  2:1845748/1‑76 (MQ=255)
                                              ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                        >  2:1764261/1‑76 (MQ=255)
                                              ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                        <  2:1558319/76‑1 (MQ=255)
                                               tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                       >  2:1110642/1‑76 (MQ=255)
                                               tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                       >  2:696097/1‑76 (MQ=255)
                                                aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                      <  1:840132/76‑1 (MQ=255)
                                                  aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                    >  1:1540779/1‑76 (MQ=255)
                                                  aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGAGCAGAGGATAAATTGTTCAATAACGGATTa                    <  2:71936/76‑1 (MQ=255)
                                                   ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                   <  1:525379/76‑1 (MQ=255)
                                                   ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                   <  1:345027/76‑1 (MQ=255)
                                                   ccGCTTCACCTTGAAAAAACTACATTGTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                   <  2:1306222/76‑1 (MQ=255)
                                                    cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                  >  1:731672/1‑76 (MQ=255)
                                                      cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                >  2:651941/1‑76 (MQ=255)
                                                        tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              >  2:1471886/1‑76 (MQ=255)
                                                        tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              <  2:704418/76‑1 (MQ=255)
                                                        tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              <  2:461654/76‑1 (MQ=255)
                                                        tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              >  1:235035/1‑76 (MQ=255)
                                                         cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc             >  1:1080495/1‑76 (MQ=255)
                                                          aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt            >  2:866139/1‑76 (MQ=255)
                                                           ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg           >  1:832429/1‑76 (MQ=255)
                                                               gAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa       >  1:238012/1‑76 (MQ=255)
                                                                aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa      <  1:1281745/76‑1 (MQ=255)
                                                                aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa      <  2:220373/76‑1 (MQ=255)
                                                                  aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg    >  2:880286/1‑76 (MQ=255)
                                                                    aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg  >  1:1190979/1‑76 (MQ=255)
                                                                            |                                                                     
TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG  >  REL606/3875551‑3875693

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