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.T95.8% 190.3 / 1.1 48intergenic (‑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/0);  new base T (27/19);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875553‑3875695
                                                                          |                                                                       
gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttt                                                                      <  1:893146/76‑1 (MQ=255)
 cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttt                                                                     >  2:482444/1‑76 (MQ=255)
  attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttt                                                                    <  2:388207/76‑1 (MQ=255)
   ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                   >  2:1366769/1‑76 (MQ=255)
   ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAtttttt                                                                   >  1:58065/1‑76 (MQ=255)
    tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                                  >  2:966160/1‑76 (MQ=255)
    tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                                  >  1:1758262/1‑76 (MQ=255)
    tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttt                                                                  <  1:1958988/76‑1 (MQ=255)
      ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttc                                                                >  1:750742/1‑74 (MQ=255)
       tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACAttttttttca                                                               >  1:1532667/1‑73 (MQ=255)
          aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                             >  1:1112762/1‑76 (MQ=255)
          aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                             <  1:1232092/76‑1 (MQ=255)
            aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                           >  2:961872/1‑76 (MQ=255)
             tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                          >  2:758774/1‑76 (MQ=255)
             tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa                                                          >  2:665699/1‑76 (MQ=255)
               cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAAt                                                       >  1:1085497/1‑76 (MQ=255)
                gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAAt                                                       <  2:1389155/76‑1 (MQ=255)
                   tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGa                                                    <  2:305124/76‑1 (MQ=255)
                     cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                                  >  2:314625/1‑76 (MQ=255)
                      tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACag                                                 >  2:843282/1‑76 (MQ=255)
                        aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                               >  1:950339/1‑76 (MQ=255)
                          tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                            <  2:1647420/76‑1 (MQ=255)
                           gcnaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           >  2:1498482/4‑76 (MQ=255)
                           gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           <  1:1051584/76‑1 (MQ=255)
                           gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           >  2:627495/1‑76 (MQ=255)
                            tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa                                          >  2:885372/1‑76 (MQ=255)
                             aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         <  1:184322/76‑1 (MQ=255)
                              aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                        <  1:140923/76‑1 (MQ=255)
                              aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                        <  2:522919/76‑1 (MQ=255)
                              aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                        >  1:1952013/1‑76 (MQ=255)
                               tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       <  2:1635745/76‑1 (MQ=255)
                               tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       >  2:1867952/1‑76 (MQ=255)
                                 tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                     <  1:398463/76‑1 (MQ=255)
                                 tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                     <  1:684712/76‑1 (MQ=255)
                                   agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGtt                                   <  1:226894/76‑1 (MQ=255)
                                          aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            >  1:903138/1‑76 (MQ=255)
                                          aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            <  1:1380911/76‑1 (MQ=255)
                                          aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            >  1:807296/1‑76 (MQ=255)
                                           tttAAACCGCTTCACCTTGAAAAAACTACCTTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           >  2:1448845/1‑76 (MQ=255)
                                            ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                          <  1:566937/76‑1 (MQ=255)
                                            ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                          >  1:1107960/1‑76 (MQ=255)
                                            ttAAACAGCTTCACCTTGAAAAAACTACATTTTTTTTCATGGGGAAATGACCGAGGAATAATTTGTTCATTaagga                          >  2:346860/1‑73 (MQ=255)
                                             tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                         >  1:690551/1‑76 (MQ=255)
                                             tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                         >  1:1326528/1‑76 (MQ=255)
                                              aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        <  1:1258034/76‑1 (MQ=255)
                                               aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                       >  2:1656692/1‑76 (MQ=255)
                                                   gCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                   <  2:535266/76‑1 (MQ=255)
                                                      tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                <  1:421340/76‑1 (MQ=255)
                                                       cACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc               >  1:134191/1‑76 (MQ=255)
                                                        aCCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt              <  2:911766/76‑1 (MQ=255)
                                                          cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa            >  1:425713/1‑76 (MQ=255)
                                                           ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           >  2:2004964/1‑76 (MQ=255)
                                                              aaaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa        <  1:187986/76‑1 (MQ=255)
                                                               aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa       >  2:563116/1‑76 (MQ=255)
                                                               aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa       >  2:1147124/1‑76 (MQ=255)
                                                                 aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc     >  2:1066563/1‑76 (MQ=255)
                                                                  aaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg    >  2:366154/1‑76 (MQ=255)
                                                                    cTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  1:892522/1‑76 (MQ=255)
                                                                          |                                                                       
GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875553‑3875695

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