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.5% 84.4 / 6.7 26intergenic (‑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/2);  new base T (12/11);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875553‑3875695
                                                                         |                                                                       
gCATTATTCAGCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                                                                      <  2:8371/76‑1 (MQ=255)
   ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                   >  2:433350/1‑76 (MQ=255)
     attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                                >  1:3937/1‑75 (MQ=255)
      ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                               <  1:629987/76‑3 (MQ=255)
        cAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATc                                                             <  2:402827/76‑1 (MQ=255)
          aCATACGCGTCCCGTCTCTAATTGNAGAACGAATTTAAACCGCTTCACCTTGAAACAANTACATTTTTTTTCATCg                                                            <  2:655227/76‑1 (MQ=255)
            aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                          >  1:308010/1‑76 (MQ=255)
               cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                                       <  1:520278/76‑1 (MQ=255)
                  gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    <  1:143576/76‑1 (MQ=255)
                     cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                                 <  1:182962/76‑1 (MQ=255)
                       gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAACAACTACATTTTTTTTCATCGTGAAATGGACaga                                               >  2:644559/1‑76 (MQ=255)
                           gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           >  2:595772/1‑76 (MQ=255)
                            tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                          <  1:13486/76‑1 (MQ=255)
                               tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       <  2:353580/76‑1 (MQ=255)
                                  gagaACGAATTTAAACCGCTTCACNTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                    >  2:134522/1‑76 (MQ=255)
                                   agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                                   >  2:530467/1‑76 (MQ=255)
                                   agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                                   >  1:225827/1‑76 (MQ=255)
                                       cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                               >  2:339211/1‑76 (MQ=255)
                                        gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAATGTCCAATa                              <  2:87008/76‑1 (MQ=255)
                                          aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            >  1:412346/1‑76 (MQ=255)
                                          aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           >  2:208188/1‑76 (MQ=255)
                                            ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                          >  1:646766/1‑76 (MQ=255)
                                             tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                         <  2:539509/76‑1 (MQ=255)
                                              aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        <  1:311084/76‑1 (MQ=255)
                                                      tCACCTAGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  1:121236/1‑76 (MQ=255)
                                                         ccTTGAGAAAACTACATTTTTTNTCATCGTGAANTGGACAGANGATAAATTGTTCAATAACTGATNATCCCCGCTg             <  2:686044/76‑1 (MQ=255)
                                                               aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCACTGACGaaa       <  2:397560/76‑1 (MQ=255)
                                                                 aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc     >  1:150936/1‑76 (MQ=255)
                                                                  aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg    >  2:21511/1‑76 (MQ=255)
                                                                    cTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  <  1:45330/76‑1 (MQ=255)
                                                                         |                                                                       
GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875553‑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: