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.T96.9% 111.7 / ‑2.1 32intergenic (‑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 (16/15);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875556‑3875695
                                                                      |                                                                        
ttattCAACATACGCGTCCTGACTGTAATTTCAGAACGAANNTAAACCGCTTNNCCTTGAAAAAACTACAtttttt                                                                   >  1:1182521/1‑76 (MQ=255)
ttattCAACATACGCGTCCTGACTGCAATTTTAGAACGATNNTAAACCGCTTNNCATTGAAAAAACTACAtttttt                                                                   >  1:1182512/1‑76 (MQ=255)
    tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                              >  1:641651/1‑73 (MQ=255)
    tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                              <  2:1170713/76‑4 (MQ=255)
       aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                            <  1:293429/76‑1 (MQ=255)
        cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                           >  2:1033941/1‑76 (MQ=255)
           aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                        <  1:451247/76‑1 (MQ=255)
             gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                      <  1:102169/76‑1 (MQ=255)
               gTCCTGACTNTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    >  1:31073/1‑76 (MQ=255)
               gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    >  1:1117396/1‑76 (MQ=255)
                   tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                                <  1:430169/76‑1 (MQ=255)
                    gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                               >  2:1297888/1‑76 (MQ=255)
                       tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            >  1:867926/1‑76 (MQ=255)
                        gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           >  2:424121/1‑76 (MQ=255)
                        gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           <  2:642713/76‑1 (MQ=255)
                         tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                          >  1:986008/1‑76 (MQ=255)
                         tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                          >  1:985993/1‑76 (MQ=255)
                          aaTTTGAGAACGGATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         <  1:393325/76‑1 (MQ=255)
                          aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         >  1:182350/1‑76 (MQ=255)
                             ttGAGAACTAATTAAAACCACTTCCAGTAGAAAAAACTTCATTTTTTATAAATATAAAATGGACAGAGGATAAAtt                                      <  2:533938/76‑1 (MQ=255)
                             ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      <  2:533930/76‑1 (MQ=255)
                              tGAGAACGAATTTAAACAGCTNCACCTNGAAAAAACTACATGTTTTTTNATCCTGAAATGAAACACGGATAAATTg                                     >  1:1276731/1‑76 (MQ=255)
                                  aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGNCAGAGGATAAATTGTTCa                                 <  1:1250420/76‑1 (MQ=255)
                                    cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                               <  1:734799/76‑1 (MQ=255)
                                     gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGNAATGGACAGAGGATAAATTGTTCAATa                              <  1:1314166/76‑1 (MQ=255)
                                     gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              <  1:808118/76‑1 (MQ=255)
                                       aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            <  1:612524/76‑1 (MQ=255)
                                        tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCNTGAAATGGACAGAGGATAAGTTGTTCAATAACg                           <  2:1112039/76‑1 (MQ=255)
                                        tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTGTCATCGTGAAATGGACAGAGGATAAATTGTNCAATAACg                           >  1:1155796/1‑76 (MQ=255)
                                           aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        >  2:866274/1‑76 (MQ=255)
                                                   tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  1:930813/1‑76 (MQ=255)
                                                   tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  2:465988/1‑76 (MQ=255)
                                                   tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  2:744118/1‑76 (MQ=255)
                                                      ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg             <  1:290096/76‑1 (MQ=255)
                                                        ttGAAAAAACTACATTTTTTTTCATCGTCAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           >  2:672611/1‑76 (MQ=255)
                                                              aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc     >  2:689666/1‑76 (MQ=255)
                                                               aaCTACATTTTTTTTCATCGTGAGATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg    <  1:1125060/76‑1 (MQ=255)
                                                                 cTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  1:189777/1‑76 (MQ=255)
                                                                      |                                                                        
TTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875556‑3875695

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