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.0% 119.9 / ‑2.4 33intergenic (‑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 (13/19);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG  >  REL606/3875550‑3875693
                                                                            |                                                                     
aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                        >  1:819159/1‑76 (MQ=255)
 tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAt                                                                      >  2:675337/1‑76 (MQ=255)
  aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                     <  2:836344/76‑1 (MQ=255)
  aGCATTATTCAACATACGCGTCCTGACTGGAATTTGGGAACGAATTTAAACCGCTTCNCCTTGAAAAAACTACAtt                                                                     <  1:1361406/76‑1 (MQ=255)
    cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttt                                                                  >  1:1137291/1‑76 (MQ=255)
     attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                 >  1:304906/1‑76 (MQ=255)
      ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                <  1:136755/76‑1 (MQ=255)
      ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                <  1:837567/76‑1 (MQ=255)
          tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                            >  1:1141774/1‑73 (MQ=255)
               aTACTCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                        >  2:1408493/1‑76 (MQ=255)
                tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                       >  2:332103/1‑76 (MQ=255)
                 gcgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                      <  1:741363/75‑1 (MQ=255)
                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                      >  1:868674/1‑76 (MQ=255)
                         tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                              <  1:223037/76‑1 (MQ=255)
                         tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                              <  1:1375534/76‑1 (MQ=255)
                         tGACTGTAATTTGAGAACGAACTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                              <  2:213216/76‑1 (MQ=255)
                          gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                             >  1:440062/1‑76 (MQ=255)
                           aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                            <  1:1053343/76‑1 (MQ=255)
                              gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                         <  1:924142/76‑1 (MQ=255)
                              gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                         <  1:280049/76‑1 (MQ=255)
                               tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                        >  1:964585/1‑76 (MQ=255)
                                aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                       <  2:881516/76‑1 (MQ=255)
                                aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                       <  2:229288/76‑1 (MQ=255)
                                 aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                      <  1:293609/76‑1 (MQ=255)
                                          cGAATTTAAACCGCNTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaan                             <  2:199397/76‑2 (MQ=255)
                                            aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATan                           <  1:1394348/76‑2 (MQ=255)
                                            aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                           >  2:745428/1‑76 (MQ=255)
                                                  aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                     <  1:561522/76‑1 (MQ=255)
                                                  aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                     >  2:319220/1‑76 (MQ=255)
                                                   acagCTTCGCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                    <  2:1055270/73‑1 (MQ=255)
                                                   aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                    <  2:1055254/76‑1 (MQ=255)
                                                    ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                   <  1:530162/76‑1 (MQ=255)
                                                      gCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                 <  2:145590/76‑1 (MQ=255)
                                                         tCACATTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              >  2:1409572/1‑76 (MQ=255)
                                                             cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa          >  1:358448/1‑76 (MQ=255)
                                                             cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa          <  1:1029635/76‑1 (MQ=255)
                                                              ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc         >  1:288634/1‑76 (MQ=255)
                                                               tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg        >  1:155390/1‑76 (MQ=255)
                                                               tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg        <  1:495206/76‑1 (MQ=255)
                                                                 aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa      >  2:1141211/1‑76 (MQ=255)
                                                                     aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg  >  1:440189/1‑76 (MQ=255)
                                                                     aaCTACATTATTTTTCATCGTGANATGGANAGAGGATAAATTGTNCACTAATGGATTATCCACGCTGACGAAAgcg  >  2:189662/1‑76 (MQ=255)
                                                                            |                                                                     
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG  >  REL606/3875550‑3875693

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