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.T86.7% 88.5 / 9.0 30intergenic (‑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 (12/14);  total (14/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 = 9.90e-01

AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875552‑3875692
                                                                          |                                                                       
aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                      >  1:853861/1‑76 (MQ=255)
 gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                                                                     >  2:155999/1‑76 (MQ=255)
   attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                  >  2:888754/1‑76 (MQ=255)
   attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                  >  2:456534/1‑76 (MQ=255)
    ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                 >  1:407774/1‑76 (MQ=255)
             aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                       <  1:279709/76‑1 (MQ=255)
             aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                       >  1:616992/1‑76 (MQ=255)
              tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                      >  1:802697/1‑76 (MQ=255)
               aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                     >  2:453668/1‑76 (MQ=255)
               aCGCGCCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                     >  2:412065/1‑76 (MQ=255)
                      cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                              <  2:575720/76‑1 (MQ=255)
                       tGACTGTAGTTTGAGAACGAATTTAAACCGCTTCACCTTGAANAAACTACATTTTTTTTCATCGTGAAATGGACag                                             >  1:962340/1‑76 (MQ=255)
                       tGACTGTAATTTGAGAACGAATTTGAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                             <  1:361257/76‑1 (MQ=255)
                           tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAggn                                         >  1:1041592/1‑75 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                        >  1:423804/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                        >  1:75493/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATNTAAACCGCTTCACCTTGAAAAAACTACATTTTTNCTCATCTTGAAATGGACAGAGGAt                                        >  2:983239/1‑76 (MQ=255)
                             tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                       <  1:736839/76‑1 (MQ=255)
                             tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                       <  2:381599/76‑1 (MQ=255)
                              aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                      <  1:974743/76‑1 (MQ=255)
                                tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                    >  1:2912/1‑76 (MQ=255)
                                  tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                  >  2:971621/1‑76 (MQ=255)
                                       aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                             <  2:532095/76‑1 (MQ=255)
                                         gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATCAATTGTTCAATa                           <  2:174466/76‑1 (MQ=255)
                                         gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                           >  1:669030/1‑76 (MQ=255)
                                         gAATTTAAACCGCTTCACCTTGAAAAAACTACANTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                           <  1:968783/76‑1 (MQ=255)
                                            tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        >  1:27096/1‑76 (MQ=255)
                                            tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        <  2:60511/76‑1 (MQ=255)
                                            tttAAACCGCTTCACCTTGAAAAAACTACATTTATTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        <  1:950918/76‑1 (MQ=255)
                                              tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        <  2:1161863/74‑1 (MQ=255)
                                              tAAACCGCTTCACCTTGAAAAAACTACATCTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                        >  1:1161864/1‑74 (MQ=255)
                                               naaCTGCTTCACCTTGAAAAANCTACATNTTTTTCCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                     <  1:994265/75‑1 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                    <  2:236226/76‑1 (MQ=255)
                                                            ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc        <  1:411811/76‑1 (MQ=255)
                                                                 aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg   <  1:412462/76‑1 (MQ=255)
                                                                  aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  <  1:568972/76‑1 (MQ=255)
                                                                          |                                                                       
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875552‑3875692

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