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.8% 174.0 / ‑2.5 45intergenic (‑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/0);  new base T (24/20);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875555‑3875692
                                                                       |                                                                    
attattCAATATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                 >  2:261943/1‑76 (MQ=255)
attattCAACATACGCGTCCTGACTGTAATTTGAGAACGGATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                 <  1:1744868/76‑1 (MQ=255)
  tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                              <  2:767504/76‑1 (MQ=255)
  tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                              >  1:47089/1‑76 (MQ=255)
  tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                              >  1:941239/1‑76 (MQ=255)
    ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                            <  1:1172043/76‑3 (MQ=255)
    ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                            >  2:1773567/1‑74 (MQ=255)
     tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                           <  1:1426359/76‑4 (MQ=255)
        aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                         <  1:250079/76‑1 (MQ=255)
          aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                       >  1:481279/1‑76 (MQ=255)
            aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                     <  2:1755749/76‑1 (MQ=255)
                gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                 <  2:1863740/76‑1 (MQ=255)
                  ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                               >  1:1739482/1‑76 (MQ=255)
                     gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAACAACTACATTTTTTTTCATCGTGAAATGTACagg                                            >  1:1225218/1‑75 (MQ=255)
                      aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                           >  2:182820/1‑76 (MQ=255)
                       cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                          >  1:320813/1‑76 (MQ=255)
                        tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                         >  1:1705612/1‑76 (MQ=255)
                        tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                         >  2:826213/1‑76 (MQ=255)
                           aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                      <  1:458654/76‑1 (MQ=255)
                             tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                    >  1:945812/1‑76 (MQ=255)
                             tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                    >  1:1615555/1‑76 (MQ=255)
                              ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                   >  1:1110595/1‑76 (MQ=255)
                                 agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                                >  2:523604/1‑76 (MQ=255)
                                   aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                              <  1:1174408/76‑1 (MQ=255)
                                   aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                              <  1:1424146/76‑1 (MQ=255)
                                    aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                             <  2:1987842/76‑1 (MQ=255)
                                    aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                             <  1:519519/76‑1 (MQ=255)
                                      gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                           <  2:1755951/76‑1 (MQ=255)
                                         tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGCAATGGACAGAGGATAAATTGTTCAATAACg                        >  1:268569/1‑76 (MQ=255)
                                         tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        >  2:1331706/1‑76 (MQ=255)
                                         tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        >  1:268580/1‑76 (MQ=255)
                                         tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        <  2:418893/76‑1 (MQ=255)
                                           tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                      <  2:1533050/76‑1 (MQ=255)
                                             aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                    <  1:440554/76‑1 (MQ=255)
                                             aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                    <  2:1705612/76‑1 (MQ=255)
                                                cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                 <  1:114720/76‑1 (MQ=255)
                                                  cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc               >  1:830314/1‑76 (MQ=255)
                                                  cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc               >  1:660318/1‑76 (MQ=255)
                                                  cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc               >  1:1551073/1‑76 (MQ=255)
                                                   ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc              <  2:114848/76‑1 (MQ=255)
                                                    tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg             >  1:1514183/1‑76 (MQ=255)
                                                    tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg             >  2:1490001/1‑76 (MQ=255)
                                                    tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg             <  1:14787/76‑1 (MQ=255)
                                                     cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAGTTGTTCAATAACGGATTATCCCCGc            >  2:395971/1‑76 (MQ=255)
                                                        cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa         <  1:436459/76‑1 (MQ=255)
                                                        cTTGAAAAAACTACATTTTTCTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa         >  1:530245/1‑76 (MQ=255)
                                                         ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc        >  1:41934/1‑76 (MQ=255)
                                                          tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg       >  2:1340643/1‑76 (MQ=255)
                                                           gAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa      >  2:985857/1‑76 (MQ=255)
                                                            aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa     <  2:1991355/76‑1 (MQ=255)
                                                            aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa     <  1:1060548/76‑1 (MQ=255)
                                                              aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg   >  2:1994499/1‑76 (MQ=255)
                                                               aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  <  1:603640/76‑1 (MQ=255)
                                                               aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  >  1:1223944/1‑76 (MQ=255)
                                                                       |                                                                    
ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875555‑3875692

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