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.1% 129.4 / ‑2.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 . (1/0);  new base T (20/14);  total (21/14)
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.73e-01

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC  >  REL606/3875550‑3875694
                                                                            |                                                                      
aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                         <  1:366780/76‑1 (MQ=255)
  aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGGATTTAAACCGCTTCACCTTGAAAAAACTACAtt                                                                     <  2:1200637/76‑1 (MQ=255)
   gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                                                                    <  2:738893/76‑1 (MQ=255)
       tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                                <  2:398918/76‑1 (MQ=255)
        attTAACATACGCGTCCTGACTGTAATTTTAGAACGAATTTAAAACGCTTCACCTTGAAAAAACTACAtttttttt                                                               >  2:337239/1‑75 (MQ=255)
        attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                               >  2:602393/1‑75 (MQ=255)
         ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                              <  1:684132/76‑3 (MQ=255)
               aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                         <  1:1374473/76‑1 (MQ=255)
                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa                                                       >  1:752611/1‑76 (MQ=255)
                     gTCCTGACTGTAATTTGAGAACGAATTTAAACNGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAANTgg                                                   >  2:1303468/1‑76 (MQ=255)
                       ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                 >  1:1135375/1‑76 (MQ=255)
                       ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAc                                                 >  1:279391/1‑76 (MQ=255)
                         tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                               <  1:891578/76‑1 (MQ=255)
                           aCNGTAATTTGAGAACGAATTTAAAGCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGATagag                                             <  2:1321627/76‑1 (MQ=255)
                            cTGTAATTTGAGAANGAATTTAAACCGCATCACCNTGAAACAANTACATTTTTTTTCATCGTGGAATGGACAGAgg                                            <  2:1162491/76‑1 (MQ=255)
                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                           <  1:369936/76‑1 (MQ=255)
                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                           >  1:1114915/1‑76 (MQ=255)
                                 aTTTGAGAACGAATTGAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                       <  2:47055/76‑1 (MQ=255)
                                   ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                     <  2:1119111/76‑1 (MQ=255)
                                   ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                     >  2:159778/1‑76 (MQ=255)
                                     gagaACGAGTGTAAACNGATTCACCTTGAAGANACTACATATTTTTTCATCGTGAAATGGACAGANGATAACntnt                                   >  1:104384/1‑76 (MQ=255)
                                     gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                   >  1:1100825/1‑76 (MQ=255)
                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                >  2:382144/1‑76 (MQ=255)
                                         aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                               <  1:799771/76‑1 (MQ=255)
                                          cGAATTTAAACCGCTTCACCTTGAACAAACTCCATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                              >  2:1130445/1‑76 (MQ=255)
                                             aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                           <  2:1114915/76‑1 (MQ=255)
                                             aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                           >  2:417394/1‑76 (MQ=255)
                                             aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                           >  2:482883/1‑76 (MQ=255)
                                               ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATTGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                         >  2:1354022/1‑76 (MQ=255)
                                                tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                        >  2:938936/1‑76 (MQ=255)
                                                  aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                      <  1:894489/76‑1 (MQ=255)
                                                   aCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                     <  1:1139284/76‑1 (MQ=255)
                                                      gCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                  <  1:1366689/76‑1 (MQ=255)
                                                       cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                 <  1:403445/76‑1 (MQ=255)
                                                        ttCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc                >  1:66350/1‑76 (MQ=255)
                                                            ccTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg            >  1:1019495/1‑76 (MQ=255)
                                                             cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGANTATCCCCGCTGa           <  1:19313/76‑1 (MQ=255)
                                                               tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg         >  1:1170648/1‑76 (MQ=255)
                                                                  aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa      >  2:779654/1‑76 (MQ=255)
                                                                    aaaCTACATTTTTTTTCATCGTGAAATGGACANAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGANAgc    >  2:1300075/1‑76 (MQ=255)
                                                                     aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg   >  1:625796/1‑76 (MQ=255)
                                                                     aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg   <  1:104845/76‑1 (MQ=255)
                                                                      aCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc  >  1:985027/1‑76 (MQ=255)
                                                                            |                                                                      
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC  >  REL606/3875550‑3875694

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