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.T100.0% 61.0 / NA 18intergenic (‑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 . (0/0);  new base T (12/6);  total (12/6)

TTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACG  >  REL606/3875556‑3875687
                                                                      |                                                              
ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                          >  2:381764/1‑76 (MQ=255)
    tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                      >  1:406305/1‑73 (MQ=255)
    tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                      >  1:406306/1‑73 (MQ=255)
         aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                  >  2:168817/1‑76 (MQ=255)
            cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                               <  1:361280/76‑1 (MQ=255)
              cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                             >  2:360287/1‑76 (MQ=255)
               gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                            <  1:481527/76‑1 (MQ=255)
                  cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                         <  2:427856/76‑1 (MQ=255)
                  cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                         >  1:74548/1‑76 (MQ=255)
                   tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                        >  2:305758/1‑76 (MQ=255)
                         tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                  >  1:482001/1‑76 (MQ=255)
                           aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaac                                >  1:287171/1‑75 (MQ=255)
                           aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                >  1:287177/1‑76 (MQ=255)
                                agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                           >  2:506140/1‑76 (MQ=255)
                                agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                           >  2:506136/1‑76 (MQ=255)
                                      aaTTTAAACCGCTTCACCTTGAAAAATCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                     <  2:222768/76‑1 (MQ=255)
                                      aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                     <  1:137013/76‑1 (MQ=255)
                                             aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa              >  1:669970/1‑76 (MQ=255)
                                                   ncACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg        >  1:344299/2‑76 (MQ=255)
                                                       cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa    <  1:714/76‑1 (MQ=255)
                                                         tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg  >  1:482463/1‑76 (MQ=255)
                                                                      |                                                              
TTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACG  >  REL606/3875556‑3875687

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