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.6% 155.6 / ‑5.6 43intergenic (‑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);  major base T (29/12);  minor base C (1/0);  total (31/12)
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.95e-01

AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC  >  REL606/3875552‑3875694
                                                                          |                                                                      
aGCCTTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                      >  1:1247215/1‑76 (MQ=255)
    ttattCCACATTCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                  >  1:1223714/1‑76 (MQ=255)
       ttCAANACACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAANCTACAttttttttc                                                              >  2:1341544/1‑74 (MQ=255)
       ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                              >  2:601287/1‑74 (MQ=255)
            cATACGCGTCCTNACNGTAATTTGAGAACGAAGTNAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                          <  1:1086346/76‑1 (MQ=255)
            cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTATATTTTTTTTCATCGt                                                          >  1:1102036/1‑76 (MQ=255)
             aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                         <  1:103819/76‑1 (MQ=255)
              tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                        <  2:1166165/76‑1 (MQ=255)
               aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                       >  2:1148752/1‑76 (MQ=255)
                cgcgTCCTGACTGTACTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                                      >  1:122798/1‑76 (MQ=255)
                cgcgTCCTGACTGTAATTTGACAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                                      <  1:647802/76‑1 (MQ=255)
                   gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                   >  2:831647/1‑76 (MQ=255)
                     ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                 >  1:914957/1‑76 (MQ=255)
                         aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                             <  1:981964/76‑1 (MQ=255)
                           tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                           >  1:301710/1‑76 (MQ=255)
                            gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                          >  1:1038978/1‑76 (MQ=255)
                              aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                        >  1:407002/1‑76 (MQ=255)
                               aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                       >  1:584594/1‑76 (MQ=255)
                               aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                       >  1:842074/1‑76 (MQ=255)
                                 ttGAGAACGAATTTACACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                     >  1:127038/1‑76 (MQ=255)
                                 ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                     >  1:628861/1‑76 (MQ=255)
                                   gagaACGAATTTAAACCGCTTCACCTTNAAAAAACTACACTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                   >  1:1167924/1‑76 (MQ=255)
                                   gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                   >  2:690242/1‑76 (MQ=255)
                                            tttAAACCGCTTCACCTTGAAAAAACTNCATTTTTTGTCTTCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                          >  1:1181535/1‑76 (MQ=255)
                                             ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                         >  1:157552/1‑76 (MQ=255)
                                               aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                       <  2:834700/76‑1 (MQ=255)
                                                 aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                     >  2:407249/1‑76 (MQ=255)
                                                 aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                     <  1:941432/76‑1 (MQ=255)
                                                  ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt                    >  1:904652/1‑76 (MQ=255)
                                                   cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                   <  1:298206/76‑1 (MQ=255)
                                                    gCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                  >  2:1043063/1‑76 (MQ=255)
                                                     cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                 >  1:783345/1‑76 (MQ=255)
                                                      ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc                <  2:1106814/76‑1 (MQ=255)
                                                      ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc                >  2:863420/1‑76 (MQ=255)
                                                       tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg               >  2:550346/1‑76 (MQ=255)
                                                        cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc              >  1:837016/1‑76 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa           >  1:823936/1‑76 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa           >  2:174017/1‑76 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa           >  1:422150/1‑76 (MQ=255)
                                                           cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa           >  1:1109303/1‑76 (MQ=255)
                                                            ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGGTAAATTGTTCAATAACGGATTGTCCCCGCTGAc          <  1:1228834/76‑1 (MQ=255)
                                                            ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc          >  1:652706/1‑76 (MQ=255)
                                                             tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGCTAAATTGTTCAATAACGGATTATCCCCGCTGACg         <  2:92634/76‑1 (MQ=255)
                                                                 aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg     >  1:815078/1‑76 (MQ=255)
                                                                  aaaNTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc    >  1:1208230/1‑76 (MQ=255)
                                                                   aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg   <  2:933049/76‑1 (MQ=255)
                                                                    aCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc  >  1:500578/1‑76 (MQ=255)
                                                                          |                                                                      
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC  >  REL606/3875552‑3875694

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