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.T92.9% 150.9 / 3.7 42intergenic (‑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 . (2/1);  new base T (22/17);  total (24/18)
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.96e-01

ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875555‑3875695
                                                                       |                                                                         
attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                    >  2:695731/1‑76 (MQ=255)
  tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                                 >  1:665300/1‑76 (MQ=255)
  tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAAATACAttttttt                                                                 >  1:665316/1‑76 (MQ=255)
   attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                                >  2:1066384/1‑75 (MQ=255)
   attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                                <  2:1140513/76‑2 (MQ=255)
   attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttttt                                                                <  2:1304928/76‑2 (MQ=255)
    ttCAACATACGCGTCCTGACTGTAAGTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                               <  1:596233/76‑3 (MQ=255)
        aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                            <  1:1499835/76‑1 (MQ=255)
         cATACGCGTCCTGCCTGTAATTTGAGAACGAATNTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                           >  2:1373734/1‑76 (MQ=255)
                gTCCTGACTGTAATTTGAGNACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    >  1:1513293/1‑76 (MQ=255)
                gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATgg                                                    <  1:1212249/76‑1 (MQ=255)
                   cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCTTGAAATGGACa                                                 <  2:321300/76‑1 (MQ=255)
                       cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAgg                                             >  2:1431606/1‑76 (MQ=255)
                        tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            <  2:1090871/76‑1 (MQ=255)
                        tGTAATTTGAGAACGAATTTAAACCGCNTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            >  1:1613193/1‑76 (MQ=255)
                         gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAANTGGACAGAGGAt                                           <  1:13958/76‑1 (MQ=255)
                           aaTTTGAGAACGAATTTAAACCGCTTCACCTTGATATAACTACATTTNTTCTCATCGTGAAATGGACAGAGGATaa                                         >  1:133617/1‑76 (MQ=255)
                            aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                        >  2:1037550/1‑76 (MQ=255)
                            aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                        >  2:1037559/1‑76 (MQ=255)
                              ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      >  2:38740/1‑76 (MQ=255)
                              ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      <  2:735432/76‑1 (MQ=255)
                                 agaACGAATTTAAACCGATTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGACGATACATTNtt                                   >  2:1526374/1‑76 (MQ=255)
                                  gaACGAATTTAAACCGCGTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATNGACAGAGGATAAATTGTTc                                  <  1:1512494/76‑1 (MQ=255)
                                   aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  2:1411575/76‑1 (MQ=255)
                                   aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 >  2:402945/1‑76 (MQ=255)
                                    aCGAATTTAAACCGCTTCACCTTGAAAAATCTACATCTTTTTTCATCGTGAAATGGANAGAGGATAAATTGTTCaa                                <  1:1541487/76‑1 (MQ=255)
                                      gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              <  1:864369/76‑1 (MQ=255)
                                       aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                             <  1:1084663/76‑1 (MQ=255)
                                         tttAAACCGCTTCACCTTGAAAAAACTATATTTTCTACCATCGTGAANTGCACAGAGGATAAATTGTTCAATAACg                           >  1:138429/1‑76 (MQ=255)
                                         tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           <  1:964918/76‑1 (MQ=255)
                                          ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGAGAGAGGATAAANTGTTCAATAACgc                          >  1:207181/1‑75 (MQ=255)
                                           tAAACCGCGTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTAAATAACGGa                         >  1:54665/1‑76 (MQ=255)
                                            aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAACTGGACAGAGGATAAATTGTTCAATAACGGAt                        >  2:263263/1‑76 (MQ=255)
                                            aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        >  2:263291/1‑76 (MQ=255)
                                            aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        <  1:628275/76‑1 (MQ=255)
                                               ccGCTTCACCTTGANAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTan                     <  2:217421/76‑2 (MQ=255)
                                                 gCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                   <  1:923074/76‑1 (MQ=255)
                                                    tCACATTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  1:97109/1‑76 (MQ=255)
                                                      aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGNGGATAAATTGTTCAATAACGGATTATCCCCGCt              <  2:1365268/76‑1 (MQ=255)
                                                        cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa            >  1:435408/1‑76 (MQ=255)
                                                        cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa            <  2:936681/76‑1 (MQ=255)
                                                         ttGAAAAAACTACATTTTTTTTCATGGTGAAATGGACAGANNATAAATTGTTNNACAACGGATTATCCCCGCTGAc           >  1:1449999/1‑76 (MQ=255)
                                                         ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           <  1:503811/76‑1 (MQ=255)
                                                         ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           >  1:1496761/1‑76 (MQ=255)
                                                          tGAAAAAACTACATTTTTTTTCATCGTGAAATNGATAGAGGNTAAATAGTTCAATAACGGANTATCCCCGCTGACg          >  2:1334088/1‑76 (MQ=255)
                                                            aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa        >  2:91587/1‑76 (MQ=255)
                                                             aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa       >  1:122871/1‑76 (MQ=255)
                                                               aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAGGAAAgc     >  1:1470669/1‑76 (MQ=255)
                                                                  cTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  2:1444877/1‑76 (MQ=255)
                                                                       |                                                                         
ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875555‑3875695

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