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.T88.9% 126.9 / 9.2 36intergenic (‑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/2);  new base T (18/14);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875560‑3875695
                                                                  |                                                                         
tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                              <  2:126946/76‑4 (MQ=255)
   aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                            >  1:403520/1‑76 (MQ=255)
     aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                          <  1:804449/76‑1 (MQ=255)
      tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                         <  2:424327/76‑1 (MQ=255)
       aCGCGTCCTGACTGTAATCTGAGAACGNATNTCNACTGCTTNACATTGGAAAAACTACATTTTTTTTCATCGTGaa                                                        <  1:1098435/76‑1 (MQ=255)
        cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                                       >  1:292232/1‑76 (MQ=255)
           gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    <  1:374055/76‑1 (MQ=255)
                gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                               >  2:764626/1‑76 (MQ=255)
                 aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                              >  1:698634/1‑76 (MQ=255)
                  cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                             >  2:676946/1‑76 (MQ=255)
                    gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           >  1:999709/1‑76 (MQ=255)
                      aaTTTGAGAACGAAGTCAAACCGCTTCACCGTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         <  1:21481/76‑1 (MQ=255)
                        tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       <  1:895961/76‑1 (MQ=255)
                        tttGAGAAAGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTCTTTCATCGTGAAATCGACAGAGGATAAAt                                       >  2:29250/1‑76 (MQ=255)
                         ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      >  2:1188371/1‑76 (MQ=255)
                            agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                                   <  1:416402/76‑1 (MQ=255)
                             gaACGAATTTAANCCNCTTCACCTTGAAAAAACTNCATTTTATTTCATCGTGAAATGGACAGAGGATAAATTGTTc                                  <  1:1037010/76‑1 (MQ=255)
                               aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                                <  1:991582/76‑1 (MQ=255)
                                cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                               >  2:258051/1‑76 (MQ=255)
                                 gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              >  2:902744/1‑76 (MQ=255)
                                 gAATTTAAACAGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              >  1:104448/1‑76 (MQ=255)
                                  aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                             <  1:728932/76‑1 (MQ=255)
                                   aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            >  1:351083/1‑76 (MQ=255)
                                    tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           >  1:1024422/1‑76 (MQ=255)
                                      tAAACCGCTTCANCTNGAAAAAACTACATTTTTTNTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                         <  1:1032019/76‑1 (MQ=255)
                                      tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                         <  1:1051784/76‑1 (MQ=255)
                                       aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        <  2:595213/76‑1 (MQ=255)
                                             cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                  >  2:489980/1‑76 (MQ=255)
                                               tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  1:818280/1‑76 (MQ=255)
                                                  ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg             <  1:914826/76‑1 (MQ=255)
                                                   cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa            >  1:771043/1‑76 (MQ=255)
                                                       aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa        >  2:1129379/1‑76 (MQ=255)
                                                       aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa        <  1:587794/76‑1 (MQ=255)
                                                         aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg      >  2:149862/1‑76 (MQ=255)
                                                           aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg    <  1:489058/76‑1 (MQ=255)
                                                           aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg    <  2:842889/76‑1 (MQ=255)
                                                           aaCTACATCTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCACGCTGACGAAAgcg    >  2:27233/1‑76 (MQ=255)
                                                             cTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  2:822552/1‑76 (MQ=255)
                                                                  |                                                                         
TCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875560‑3875695

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