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.T91.4% 122.1 / 5.8 35intergenic (‑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 . (3/0);  new base T (18/14);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG  >  REL606/3875550‑3875693
                                                                            |                                                                      
aTAGCATTGTTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                         >  2:1053296/1‑76 (MQ=255)
aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                         >  2:1053285/1‑76 (MQ=255)
  aGCATTATTCAACANACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtn                                                                      <  2:168154/76‑2 (MQ=255)
       tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                               <  1:283793/76‑1 (MQ=255)
        attCAACATACGCGTCCTGACTGTAATTTGAGAACGNATTTAAATCGATTCACCTTGAAAAAACTACAtttntttt                                                              >  1:150216/1‑75 (MQ=255)
        attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                              >  1:415235/1‑75 (MQ=255)
        attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                              <  1:619736/76‑2 (MQ=255)
             aCATACGCGTCCTGACTGTAATTTCAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                          >  2:261329/1‑76 (MQ=255)
              cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                         <  2:448833/76‑1 (MQ=255)
               aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                        >  2:539399/1‑76 (MQ=255)
                   gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                    <  2:982591/76‑1 (MQ=255)
                    cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                   <  1:1051083/76‑1 (MQ=255)
                         tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                              <  1:734548/76‑1 (MQ=255)
                          gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                             <  1:383474/76‑1 (MQ=255)
                          gACTGTAATTTGAGAACGAATTTAAACCGCATCACCTTGAAAAAACTACATTTTTTTTCATCNTGAAATGGACaga                                             >  2:1108503/1‑76 (MQ=255)
                           aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                            <  2:1065461/76‑1 (MQ=255)
                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                          >  1:1013833/1‑76 (MQ=255)
                               tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                        >  2:889202/1‑76 (MQ=255)
                                aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                       >  2:518852/1‑76 (MQ=255)
                                  tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                     <  2:112955/76‑1 (MQ=255)
                                   ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                    >  1:875194/1‑76 (MQ=255)
                                   ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                    <  2:848090/76‑1 (MQ=255)
                                    tGAGAACGAATTTAAACCGCTTCACCTAGAAAAAACTACATTTGNTTNCATCGTGAAATGGANAGAGGATAAATTg                                   >  2:1091995/1‑76 (MQ=255)
                                     gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATTGTGGAATGGACAGAGGGTAAATTGt                                  >  1:786627/1‑76 (MQ=255)
                                      agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAANGGACAGAGGATAAATTGtt                                 <  1:1130251/76‑1 (MQ=255)
                                       gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                               >  2:52781/1‑76 (MQ=255)
                                         aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                              >  1:520160/1‑76 (MQ=255)
                                         aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                              >  1:520166/1‑76 (MQ=255)
                                          cGAATTTAAACCGCTTCACCTTGAAAAAACTACNTTTTTTTTCATCGTGAAATGGANAGAGGATATANTGTTCAAt                             >  2:1010495/1‑76 (MQ=255)
                                               ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                        <  1:511546/76‑1 (MQ=255)
                                                tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                       >  2:477019/1‑76 (MQ=255)
                                                 aaaCGGCTTCACCTTGACAAAACTACATTTTTTTTCATCGTGAAATGGACGGAGGATAAATTGTTCAATTACGNAt                      >  2:1141382/1‑76 (MQ=255)
                                                    ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                   <  2:495749/76‑1 (MQ=255)
                                                     cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                  >  2:1067795/1‑76 (MQ=255)
                                                       cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                >  1:950744/1‑76 (MQ=255)
                                                         tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              <  1:187662/76‑1 (MQ=255)
                                                           aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt            <  2:380310/76‑1 (MQ=255)
                                                           aCCTCGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt            >  2:1193500/1‑76 (MQ=255)
                                                            ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg           >  1:547603/1‑76 (MQ=255)
                                                               tGAAAAAACTACATTTTTTNTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACCGATTATCCCCGCTGACg        >  1:1127589/1‑76 (MQ=255)
                                                                    aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc   <  1:796997/76‑1 (MQ=255)
                                                                     aaCTACATTTTTTTTCATCGTGAAATGGACAGANGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg  >  2:1022544/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 < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: