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.T93.9% 118.6 / 0.8 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 . (2/0);  new base T (17/14);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875552‑3875692
                                                                          |                                                                    
aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                   >  1:865179/1‑76 (MQ=255)
aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                   <  1:52584/76‑1 (MQ=255)
aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                   <  1:448144/76‑1 (MQ=255)
 gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                                                                  <  2:959329/76‑1 (MQ=255)
 gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCACTTCACCTTGAAAAAACTACAttt                                                                  >  1:329270/1‑76 (MQ=255)
        tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                           >  1:1200608/1‑73 (MQ=255)
            cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                        >  1:499301/1‑76 (MQ=255)
             aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                       <  1:315718/76‑1 (MQ=255)
              tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                      >  2:2740/1‑76 (MQ=255)
              tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                      >  1:548989/1‑76 (MQ=255)
                  cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                  >  1:683159/1‑76 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                               <  2:72772/76‑1 (MQ=255)
                       tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                             <  1:596280/76‑1 (MQ=255)
                       tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATATTTTTTCATAGTGAAATGGACag                                             >  1:1113449/1‑76 (MQ=255)
                        gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                            <  1:846173/76‑1 (MQ=255)
                          cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                          <  1:827871/76‑1 (MQ=255)
                          cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                          <  1:855729/76‑1 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGANAAAACTAGATGTTTTATCATCGTGAAATGGACAGNGGAt                                        >  1:151561/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                        >  2:938170/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                        >  2:56214/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                        >  2:938161/1‑76 (MQ=255)
                              aaTTTGAGACCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                      >  2:1115588/1‑76 (MQ=255)
                                 ttGAGAATGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                   >  1:1267052/1‑76 (MQ=255)
                                    agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTCAAATGGACAGAGGATAAATTGtt                                >  1:1122226/1‑76 (MQ=255)
                                      aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                              <  1:773545/76‑1 (MQ=255)
                                           aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                         <  1:1059853/76‑1 (MQ=255)
                                             ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                       >  1:782797/1‑76 (MQ=255)
                                             ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                       >  2:1266119/1‑76 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCTATAACGGAtt                    <  1:1198979/76‑1 (MQ=255)
                                                aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                    <  2:529601/76‑1 (MQ=255)
                                                  ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                  <  1:283300/76‑1 (MQ=255)
                                                       tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg             >  1:519249/1‑76 (MQ=255)
                                                         aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt           <  1:423557/76‑1 (MQ=255)
                                                          ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg          <  1:1253795/76‑1 (MQ=255)
                                                          ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg          <  1:379554/76‑1 (MQ=255)
                                                                aaaaaCTACATTTTTTTTCATCGTGAANTGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa    >  1:1297893/1‑76 (MQ=255)
                                                                aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa    >  2:274706/1‑76 (MQ=255)
                                                                  aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  >  2:280857/1‑76 (MQ=255)
                                                                  aaaCTACATTTTTTTTCATCGTCAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  >  2:145901/1‑76 (MQ=255)
                                                                          |                                                                    
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875552‑3875692

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